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'in  which  the  flow  was  obliquely  interrogated  over  its  complete  cross  section  by 
a  square-enveloped  ultrasonic  beam.  The  cell  included  screens  for  minimiz¬ 
ing  eddies  at  the  transducer  ports,  gradual  circular  to  rectangular  pipe  tranBi-| 
tions,  flushing  means,  and  a  bleed  port.  Flow  related  parameters  were  mea¬ 
sured  by  me«Bi  of  a  novel  approaeh-  utilizing  a  fully  coherent  electronic  system! 
which  operated  synchronously  from  a  5 -MHz  crystal  controlled  oscillator.  TheJ 
transmitted  waveform  was  generated  by  means  of  synchronously  gating  the 
output  of  the  5 -MHz  oscillator  for  a  50-microsecond  interval  at  a  2-kHz  repe¬ 
tition  rate.  Upstream  and  downstream  common-path  transmissions  were 
consecutive  and  spaced  250  microseconds  apart.  Independent  measurements 
were  made  of  the  phase  difference  between  the  upstream  and  downstream 
5-MHz  spectral  components  of  the  received  waveforms  and  the  phase  sum  of 
the  2-kHz  modulation  components  of  these  waveforms.  These  two  independent 
measurements  were  combined  to  yield  the  Mach  number  v/c  and  v,  where 
v  =  flow  velocity  and  c  =  sound  speed.^High  signal-to-noise  ratio  was  achieved| 
by  means  of  narrowband  receivers  anafully  coherent  detection,  thereby  pro¬ 
viding  distinct  advantages  in  precision  and  stability  over  conventional  flow¬ 
metering  techniques  which  utilize  broadband  pulses.  One  cell  configuration 
contained  an  ultrasonic  reflectometer  which  generated  a  pair  of  echoes.  The 
echo  amplitudes  were  processed  to  yield  the  fuel's  characteristic  impedance 
pc,  where  p  =  density.  The  product  of  the  v/c  and  pc  terms  in  pv,  which 
equals  the  mass  flow  rate  M.  In  a  second  cell,  presently  the  preferred  con¬ 
figuration,  p  was  measured  downstream  of  the  velocixneter  region  in  a  serial 
branch  by  a  standard  capacitance-type  di  nsitometer.  M  then  was  determined 
as  the  pv  product.  The  flowmeter  was  programmed  to  perform  such  com¬ 
putations  digitally.  In  six  independent  J]  »-4  calibration  tests  at  50,  75  and 
100°F,  and  using  the  second  cell  and  the  corresponding  formulation  id  =  p  v, 


the  maximum  average  deviation  from  li 
the  mean  of  the  average  deviations  was 
reproducibility,  the  maximum  average 
while  the  mean  of  the  average  deviation^ 
results  demonstrated  that  the  calibrate* 
flow  velocimeter  and  a  capacitance-typ< 


arity  was  0.  2%  of  full  scale,  while 
13%  of  full  scale.  With  respect  to 
eviation  was  0.  2s%  of  full  scale, 
was  0.  2%  of  full  scale.  These  test 
p  v  system  consisting  of  an  ultrasonic 
densitometer  could  measure  M  to 
nd  JP-5  and  the  conditions  of  interest. 


within  _+  1%  of  range  for  the  fuels  JP-4 
With  appropriate  modification,  the  system  could  be  used  to  measure  dynamic 
volumetric  and  mass  flow  rates.  Respjonse  time,  presently  several  seconds, 
could  be  reduced  to  the  range  10  to  10Q  milliseconds. 
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INTRODUCTION 


This  introductory  section  includes  paragraphs  on  the  need  for  a  mass 
flowmeter,  results  of  the  previous  program,  approaches  to  solve  spe¬ 
cific  problems,  and  a  summary  of  work  accomplished  in  the  present 
program. 

NEED  FOR  MASS  FLOWMETER 


Before  proposing  that  the  present  program  be  conducted,  we  felt  it  was 
necessary  to  first  determine  whether  or  not  an  industry-wide  need  ex¬ 
isted  for  a  diagnostic  meter  to  measure  the  mass  flow  rate  M  for  aviation 
fuels.  Therefore,  we  supplemented  information  available  in  the  literature 
by  contacting  instrumentation  and  control  engineers  at  AiResearch,  Avco, 
General  Electric  and  Pratt  &  Whitney  (Table  1).  Of  ten  engineers  inter¬ 
viewed,  eight  expressed  a  definite  need  for  M  measurements  for  engine 
diagnostics  and/or  control.  One  of  the  two  objectors  presented  the  view 
that  since  engine  efficiencies  and  fuel  heating  values  vary  by  several 
percent,  thrust  cannot  be  controlled  accurately  enough  from  M  alone. 
However,  it  was  acknowledged  that  (a)  M  might  be  the  correct  parameter 
to  base  control  on,  in  future  electronic-computer-controlled  engines;  and 
(b)  M,  the  time  rate  of  change  in  M,  could  be  an  important  control  param¬ 
eter  in  fuel  scheduling,  to  avoid  undesirable  transients.  (M  =  fuel  density 
times  acceleration.  )  Another  objector  felt  that  rpm,  pressure  and  tem¬ 
perature  were  more  important  than  M,  for  control.  However,  this  indi¬ 
vidual  felt  that  a  diagnostic  need  might  become  apparent  during  the  de¬ 
velopment  and  testing  of  the  T700  engine.  Some  other  individuals  placed 
emphasis  either  on  diagnostics  or  on  control,  depending  on  their  own 
primary  responsibilities.  Several  indicated  the  need  for  M  data  in  both 
disciplines.  A  response  time  of  about  10  ms  apparently  would  satisfy 
all  those  surveyed.  An  engineer  concerned  with  large  engines  such  as 
GE's  CF6,  which  develops  a  thrust  of  40,  000  lb  and  which  burns  fuel  at 
the  rate  of  nearly  20,000  lb/hr  on  takeoff,  pointed  out  that  a  demonstra¬ 
tion  of  our  proposed  flowmeter  in  the  2000  lb/hr  range  would  be  valuable, 
and  of  course  would  be  applicable  to  many  "smaller"  engines  such  as 
GE's  J85  or  other  helicopter  engines. 

Discussions  showed  that  military  applications  are  generally  more  severe 
than  commercial.  The  military  situation  is  more  likely  to  include  fuel 
mixtures  or  fuel  substitutions,  and  contaminants.  Engineers  are  appar¬ 
ently  reluctant  to  run  tests  with  contaminated  fuel  flowing  through  new 
turbine  flowmeters,  because  of  the  bearing  wear  and  maintenance  prob¬ 
lems.  Instead,  older,  less  accurate  flowmeters  are  often  used.  Clearly, 
an  accurate  ultrasonic  mass  flowmeter  would  be  useful  in  such  tests. 
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Table  1  also  shows  that  there  is  some  interest  in  M  ranging  below  50  lb/hr, 
and  considerable  interest  up  to  5000  lb/hr.  Again,  for  larger  engines,  M 
ranges  over  25,000  lb/hr,  with  some  requirements  as  high  as  100,000 
lb/hr.  In  most  engines,  the  range  of  M  does  not  exceed  about  25:1.  For 
example,  the  GE  T700  engine  burns  about  100  to  2500  lb/hr  o  1  JP-4  or 


JP-5. 


Regarding  fuel  temperature,  ground  testing  diagnostics,  if  carried  out 
in  Alaska  or  desert  sun,  approach  the  limits  of  the  present  temperature 
span,  -65  F  and  +160  F.  High  altitude  (~30,000  ft)  generates  the  mini¬ 
mum  temperature  limit.  The  maximum  temperature  requirement  depends 
partly  on  the  engine  and  application,  on  the  meter  location,  and  on  heat-up 
after  engine  turnoff,  in  the  absence  of  cooling  by  fuel  flow.  Thus,  a  250  F 
requirement  may  be  applicable  in  some  instances. 

In  most  diagnostic  testing,  the  fuel  is  known.  Also,  a  manual  check  on 
density  and  temperature  may  not  be  objectionable.  Besides  the  diagnostic 
need  in  ground  testing,  Edwards  Air  Force  Base  and  other  sources  have 
indicated  a  need  in  flight  testing.  Here,  in  the  development  .md  evaluation 
of  new  engines,  flight  test  data  is  telemetered  to  the  ground. 

Fuel  pressure  is  often  below  100  psi,  but  can  be  >,  1000  psi,  depending  on 
meter  location,  and  so  a  2000-psi  capability  is  desired. 

The  accuracy  requirement  is  generally  felt  to  be  about  0.  5%;  but  it  was 
pointed  out  that  for  some  control  applications,  this  might  be  relaxed  to 
2%  in  the  interests  of  reliability,  simplicity,  small  size  and  light  weight. 

Lastly,  we  should  briefly  indicate  the  diagnostic  need  for  an  adequate 
mass  flowmeter,  based  on  our  own  experiences  in  conducting  calibration 
and  contamination  tests  at  Avco-Lycoming.  In  calibration  tests  the  avail¬ 
able  turbine  flowmeters  individually  covered  M  ranges  of  only  10:1.  They 
had  to  be  calibrated  on  the  particular  fuel  used  in  the  test.  That  is,  on 
switching  from  JP-4  to  JP-5,  they  had  to  be  recalibrated.  Even  on  a 
given  fuel,  changes  in  fuel  temperature  (density)  required  the  turbine 
flowmeter  reading  to  be  multiplied  by  a  density  ratio  correction  factor. 

In  contaminant  tests,  lifetimes  of  turbine  flowmeters  are  usually  less 
than  about  10  hours. 

It  seemed  reasonable  that  the  basic  ultrasonic  approach  of  the  previous 
program  could  be  extended  to  linearly  respond  to  an  M  range  of  50  to 
5000  lb/hr,  could  operate  from  -65  to  +160  F  and  at  pressures  to  2000 
psi,  and  could  respond  in  10  ms.  We  therefore  proposed  a  program  which 
generally  retained  the  positive  features  of  the  earlier  approach,  but  did 
not  retain  the  negative  elements. 


15 


RESULTS  OF  PREVIOUS  PROGRAM 


Reviewing  the  previous  program,1  five  major  results  may  be  identified: 

•  A  new  ultrasonic  transmission  flow  velocimeter  technique  for  con¬ 
tinuously  measuring  v/c^  was  demonstrated.  This  technique  was 
believed  to  inherently  offer  the  best  accuracy  consistent  with  fast 
response,  20  ms,  of  any  ultrasonic  approach  demonstrated  so  far. 
Operation  up  to  5000  lb/hr  was  demonstrated. 

•  A  new  ultrasonic  reflection  densitometer  technique  for  measuring 
p  c  was  demonstrated  on  liquids  ranging  in  density  from  0.  7  to  1.  6 
g/cm3. 

•  A  new  weighted-area-averaging  method  was  demonstrated  on  water, 
Avgas-100,  JP-4,  JP-5  and  a  contaminated  solvent,  for  interrogating 
100%  of  the  flowing  liquid's  square  cross  section  using  a  rectangularly 
collimated  plane  wave.  This  method  potentially  offers  linearity  and 
accuracy  essentially  independent  of  Reynolds  number  or  flow  profile. 

•  Using  a  cell  design  containing  nonintrusive  transducers  and  no  rotat¬ 
ing  parts,  endurance  of  over  100  hours  was  demonstrated  using  a 
test  fluid  flowing  at  ~1900  lb/hr  and  contaminated  essentially  in  con¬ 
formance  with  military  specifications  for  flow  measuring  instruments. 

•  A  new  ultrasonic  mass  flowmeter  was  built,  consisting  of  a  flow  veloc¬ 
imeter  (v/c^),  a  densitometer  (pc)  and  a  time  intervalometer  (T~l/c). 
The  combination  of  these  outputs  has  been  shown  to  yield  an  output 
voltage  approximately  proportional  to  the  mass  flow  rate  M  over 
specified  ranges. 

It  was  concluded  that  a  refinement  of  the  new  technology  demonstrated  in 
the  previous  program  should  result  in  an  ultrasonic  mass  flowmeter 
capable  of  meeting  current  U.  S.  Army  requirements  for  diagnostics, 
and,  ultimately,  engine  control. 

Our  approach  to  overcoming  the  remaining  problems  is  summarized  next. 

APPROACHES  TO  SOLVE  SPECIFIC  PROBLEMS 

Three  key  problems  were  to  be  solved  in  the  present  program: 

(1)  Improve  the  accuracy  to  4^  1%  at  2000  lb/hr,  5%  at  minimum  flow 
(50  lb/hr). 
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(2)  Eliminate  the  effects  due  to  temperature -dependence  of  impedance 
of  transducers  used  in  the  flow  velocimeter. 

(3)  Avoid  nonlinearity  at  high  flow  velocities. 

We  approached  these  problems  as  follows: 

(1 )  Accuracy 

We  wanted  to  retain  the  advantages  of  continuous  wave  as  demonstrated 
in  the  previous  program,  and  yet  avoid  the  stringent  requirements  on 
impedance  matching  associated  with  phase  coding  with  single  transmitter/ 
receiver  transducers  on  each  side  of  the  acoustic  path.  This  led  us  to 
select  the  synchronous  radio  frequency  burst  method  described  later.  The 
phase  difference  of  the  5 -MHz  carrier  is  measured  after  heterodyning. 
Additionally,  the  system  measures  the  sum  of  upstream  and  downstream 
transit  times  at  2  kHz,  which  is  the  pulse  repetition  frequency.  This 
eliminates  the  previous  requirement  for  a  separate  time  inter valometer. 
Further,  this  data  is  processed  to  yield  the  Mach  No.  ,  v/c,  which,  when 
multiplied  by  the  p  c  term,  yields  a  product  proportional  to  M,  the  mass 
flow  rate.  M  can  also  be  determined  as  the  product  of  v  and  p . 

In  the  general  pc  subsystem  under  study,  where  assumptions  about  spe¬ 
cific  fuel  type  or  specific  p-c-T  (density -sound  speed-temperature) 
correlations  cannot  be  made,  we  chose  to  measure  the  echo  amplitudes 
of  a  nearly- symmetrical  pc  reflectometer  probe  using  a  manufacturer- 
modified  version  of  their  standard  commercial  instrument.  The  instru¬ 
ment,  called  an  Echo  Comparator,  resolves  echo  amplitudes  to  0.  01  dB, 
or  3.  1%,  which  is  adequate  for  the  present  purposes.  After  some  modi¬ 
fication,  its  thermal  stability  appeared  to  be  reasonable.  The  pc  probe, 
however,  appeared  to  be  considerably  more  sensitive  to  temperature, 
with  drift  attributed  to  coupling  or  adherence  changes  in  the  bond  layers 
between  the  transducer  and  the  buffer  rods.  Several  bonding  methods 
were  evaluated,  but  an  adequate  solution  has  not  yet  been  demonstrated. 

As  of  this  writing,  it  appears  that  some  additional  development  work  is 
still  required  to  reduce  pc  errors  to  a  fraction  of  1%  over  the  tempera¬ 
ture  and  time-scale  ranges  anticipated  for  the  system.  Specific  ap¬ 
proaches  for  accomplishing  this  are  discussed  in  Appendix  C,  based  on 
either  reflection  or  transmission  coefficient  principles,  and  utilizing 
existing  electronic  equipment,  or  a  new  adaptation  of  the  narrowband 
concepts  used  in  the  flow  velocimeter.  On  the  other  hand,  a  satisfac¬ 
tory  densitometer  alternative  has  been  identified  for  JP-4  and  JP-5. 

This  alternative  is  based  on  the  p-dependence  of  the  fuel's  dielectric 
constant,  and  the  required  equipment  is  commercially  available. 
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(2)  Transducer  Improvements 


By  choosing  an  rf  burst  mode  referred  to  in  paragraph  (1)  wherein  the 
pulses  are  interlaced,  each  transducer  still  functions  as  transmitter  and 
receiver,  but  now  these  two  functions  are  separated  in  time  by  250ps. 
Therefore,  the  stringent  requirements  on  transducer  impedance  remain¬ 
ing  constant  could  be  relaxed  considerably,  without  sacrificing  electrical 
isolation  between  transmitter  and  receiver. 

Additionally,  partly  as  a  result  of  other  related  transducer  activities,  we 
departed  from  a  conventional  tungsten-loaded-epoxy  transducer  backing 
and  chose  a  graphite  backing  instead.  The  acoustic  impedance  of  graphite 
is  relatively  insensitive  to  temperature,  and  the  differential  thermal  ex¬ 
pansion  between  the  ceramic  transducer  element  and  graphite  is  relatively 
small.  Therefore,  a  more  stable,  reliable  thin  epoxy  bond  between  the 
two  can  be  achieved.  Such  bonds  have  successfully  withstood  thermal 
cycling  tests  from  liquid  nitrogen  temperature,  about  -200  C,  up  to  boil¬ 
ing  water  temperature,  +  100  C. 

(3)  Linearity 

Analysis  of  test  data  obtained  in  the  first  contract  led  us  to  conclude  that 
the  observed  departures  from  linearity  were  due  mainly  to  eddies  gener¬ 
ated  in  the  45  transducer  recesses,  and  due  secondarily  to  eddies  gen¬ 
erated  at  the  abrupt  inlet  transition  of  the  cell. 

After  comparing  several  alternatives,  we  chose  to  install  a  corrosion- 
resistant  40-mesh  gcreen  over  each  transducer  recess,  to  maintain  a 
substantially  continuous  wall,  as  far  as  the  flow  was  concerned.  Next, 
we  designed  a  circular-to- rectangular  inlet  transition  pipe. 

Test  and  calibration  data  obtained  in  June  and  December  1974  showed  that 
the  nonlinearity  previously  observed  at  high  flow  rates  was  eliminated. 

SUMMARY  OF  WORK  ACCOMPLISHED 


1.  Panametrics  furnished  the  personnel,  material,  equipment,  and  facil¬ 
ities  to  perform  experimental,  developmental  and  research  work  to  design, 
fabricate,  and  test  an  ultrasonic  mass  flowmeter  for  use  on  Army  aircraft 
gas  turbine  engines.  The  redesigned  flowmeter  was  based  on  the  flow¬ 
meter  design  previously  designed  by  us  under  Contract  No.  DAAJ02-71-C- 
0061. 1 

2.  We  performed  work  directed  toward  attaining  the  following  objectives 
as  related  to  the  ultrasonic  mass  flowmeter: 
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a.  The  flowmeter  was  designed  to  operate  over  the  temperature  range 
of  -65°F  to  +1 60°F  (-548C  to  +71°C). 

b.  The  flowmeter  was  designed  to  measure  fuel  flow  over  the  range  of 
50  pounds  per  hour  to  over  4,  000  pounds  per  hour.  M  was  to  be  de¬ 
termined  as  the  product  of  flow  velocity  and  density- related  terms. 

c.  The  transducer  cell  was  to  have  no  moving  parts  associated  with  the 
M  measurement  per  se.  However,  manually  operated  valves  are 
included  as  part  of  a  contaminant-flushing  branch. 

d.  The  transducer  cell  was  to  contain  no  significant  pressure-dropping 
restrictions  or  obstructions  in  the  flow  path. 

e.  The  transducer  cell  was  to  be  sized  and  configured  so  as  to  be 
adaptable  for  mounting  on  Army  aircraft  gas  turbine  engines. 

3.  The  work  accomplished  may  be  further  understood  by  considering  the 
following  three  tasks: 

a.  Task  I  -  Design  Optimization: 

(1)  In  this  initial  task  we  optimized  the  previous  ultrasonic  mass 
flowmeter  to  combine  the  parametric  measurements  of  a  veloc- 
imeter,  a  densitometer,  and  a  time  intervalometer,  such  that  the 
output  is  in  direct  proportion  to  mass  flow  rate.  It  turned  out  that 
the  time  intervalometer  function  could  be  accomplished  within  the 
velocimeter,  without  the  need  for  a  separate  instrument. 

(2)  In  this  effort  we  addressed  the  problems  associated  with 
transducer  ports  and  the  inlet  transition.  We  optimized  the 
design  of  the  transducer  cell  so  as  to  eliminate  the  undesirable 
effects  associated  with  stationary  eddies  and  nonlaminar  flow. 
Specifically,  we  used  screens  along  the  flow  channel  walls,  a 
gradual  inlet  transition,  and  a  mixer  just  upstream  of  the  v  mea¬ 
suring  region.  (The  need  for  the  mixer  is  not  yet  fully  established.  ) 

(3)  We  optimized  the  velocimeter  design  so  as  to  eliminate  the 
undesirable  effects  associated  with  temperature  variations  by 
eliminating  the  coded  cw  approach  in  favor  of  an  interlaced 
coherent  burst  approach.  Despite  encouraging  results  in  our 
attempts  to  eliminate  thermal  effects  in  the  pc  system,  we 
finally  interrupted  this  activity  in  favor  of  a  standard  capac¬ 
itance  type  densitometer.  This  densitometer,  as  described  in 
a  1974  publication  by  Stuart,  provides  adequate  p  accuracy  for 
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JP-4  and  JP-5  over  the  full  temperature  range  of  interest.  This 
densitometer,  or  even  other  standard  densitometer  alternatives, 
can  be  installed  downstream  of  the  v-measuring  region  in  a  serial 
or  parallel  branch. 

b.  Task  II  -  Fabrication: 

(1)  We  fabricated  one  mass  flowmeter  system  based  on  the  opti¬ 
mized  design  developed  under  Task  I.  The  mass  flowmeter  sys¬ 
tem  included  a  velocimeter,  a  densitometer,  a  signal  processor 
containing  programmed  calculator  chips  to  automatically  and 
digitally  compute  mass  flow  rate  and  other  terms  of  interest. 

(2)  The  transducer  portion  of  the  flowmeter  system,  which  was 
on  the  fuel  line  containing  the  fuel  being  metered,  was  designed 
and  fabricated  to  withstand  the  thermo-mechanical  environment 
of  an  Army  aircraft  gas  turbine  engine  installation.  This  was 
verified  by  thermal  and  vibration  tests. 

c.  Task  HI  -  Tests: 

(1)  We  conducted  verification  tests  on  the  velocimeter  and 
densitometer,  including  long  term  stability  and  drift  tests. 

(2)  We  conducted  preliminary  calibration  tests  and  endurance 
tests,  using  water  flowing  at  rates  up  to  4000  lb/hr  at  room 
temperature. 

(3)  We  tested  the  p  c  densitometer  probe  from  cryogenic  levels 
up  to  boiling  water  temperature  under  no-flow  conditions. 

(4)  We  conducted  calibration  tests,  but  not  contaminant  tests. 

The  calibration  tests  used  JP-4  and  diesel  fuel  I  VVF  800.  The 
tests  were  conducted  at  room  temperature  and  at  higher  and 
lower  temperatures  (i.  e.  ,  50,  75,  100,  115  F)  over  the  mass 
flow  rate  range  of  50  to  6000  lb/hr. 

(5)  Vibration  tests  were  conducted  on  the  pc  probe,  the  trans¬ 
ducers,  and  a  test  block  containing  a  screen  soldered  along 
the  flow  channel  sidewall,  generally  in  accordance  with  MIL- 
F-5009D.  No  significant  detrimental  effects  were  observed. 
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THEORY  OF  OPERATION 


BASIC  EQUATIONS 

The  flow  measuring  system  used  in  this  effort  differs  from  conventional 
ultrasonic  flowmetering  systems.  Many  conventional  systems  use  a  dual 
"sing -around"  technique,  in  which  a  series  of  pulses  is  transmitted  in 
first  one  direction  (e.  g.  ,  upstream)  and  then  in  the  opposite  direction 
(downstream).  The  repetition  rate  of  the  pulses  in  a  given  direction  is 
predominantly  determined  by  the  transit  time  across  the  fluid  medium: 
When  a  pulse  is  received,  this  causes  the  next  pulse  to  be  transmitted  in 
the  same  direction.  The  repetition  rate  thus  becomes  essentially  the 
reciprocal  of  the  transit  time.  Neglecting  second-order  effects,  the 
equations  appropriate  to  a  conventional  system  associated  with  a  flow 
cell  similar  to  that  shown  in  Fig.  1,  are 

fj  =  ~  sin  9  (c  +  v  cos  9  )  (1) 

f2  =  =  “jj  sin  9  (c  -  v  cos  9  )  (2) 


The  frequency  difference  is  measured,  giving 


„  sin  9  cos  9  v  .  „  „ 

2v  - g -  =  —  sin  2 0 


(3) 


where  v  =  flow  velocity,  c  =  sound  speed,  D  -  pathlength  across  the  fluid 
perpendicular  to  the  pipe  axis,  and  9  =  angle  of  transmission  relative  to 
the  pipe  axis.  The  main  difficulty  with  this  type  of  system  is  that  the  fre¬ 
quency  difference,  fj  -  f2*  is  very  small  indeed.  As  can  be  seen  by  add¬ 
ing  Eqs.  (1)  and  (2)  and  noting  that  v  <<  c,  we  see  that 


If,  for  example,  D  =  1",  9  =  45°,  and  c  =  5000  ft/sec,  we  have  from  Eqs. 
(1)  or  (2),  fj  =  f2  =  21  kHz,  and  from  Eq.  (3),  for  v  =  1  ft/sec,  fj  -  f2  = 
12  Hz.  If  the  system  is  to  have  a  resolution  of,  say,  0.  01  ft/sec,  from 
Eq.  (4)  we  see  that  the  resolution  of  the  frequency  difference  must  be 
~0.  12  Hz,  with  the  basic  sing-around  frequency  on  the  order  of  20  kHz, 
which  means  that  each  of  the  sing-around  oscillators  must  be  stable  to 
within  a  few  parts  per  million. 
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Figure  1.  Basic  Geometry. 
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In  terms  of  time  resolution,  it  can  easily  be  shown  (assuming  a  45°  path) 
that  each  pulsed  oscillator  must  have  a  basic  resolving  time  of  6  T  £ 

(D/c  )  5  v.  If,  as  before,  5  v  =  0.  01  ft/sec,  D  =  1"  and  c  =  5  x  lO^ft/sec, 
we  find  that  6T  $  3.  3  x  10“^  sec  =  33  picoseconds.  This  iB,  of  course, 
a  very  heavy  burden  upon  the  electronics.  Another  way  of  stating  the 
problem  is  to  note  that  the  required  average  phase  stability  of  each  in¬ 
dividual  oscillator  (for  6  v  =  0.  01  ft/sec,  D  =  1",  c  =  5000  ft/sec)  would 
be  ~2  x  10” ^  x  360°/2tt  =  1.  1  x  10"^  electrical  degrees. 

In  addition  to  the  above  problem  in  precision,  one  is  also  confronted  with 
a  situation  where  the  higher  the  required  precision,  the  longer  must  be 
the  measurement  time.  In  the  above  example,  a  measurement  time  on 
the  order  of  8.  3  seconds  would  be  required.  If  the  pipe  diameter  is  in¬ 
creased,  the  measurement  time  increases,  as  can  be  seen  from  Eq.  (3). 
For  example,  for  a  10"  pipe,  83  seconds  would  be  required  to  resolve  the 
flow  velocity  to  within  .01  ft/ sec. 

Many  of  the  above  problems  have  been  eliminated  in  the  present  measuring 
technique.  Rather  than  measure  one  quantity  to  extremely  high  accuracy, 
as  in  conventional  ultrasonic  flowmeters,  the  present  system  utilizes  two 
measurements,  (Tj  -  T2)  and  (Tj  +  T2),  each  of  which  need  be  accurate 
to  only  a  relative  precision  on  the  order  of  the  overall  required  instru¬ 
mental  precision.  (Measurement  of  the  latter  term,  Tj  +  may  be 
thought  of  as  accomplishing  the  objective  associated  with  the  separate 
time  intervalometer  of  the  previous  contract.  ) 


It  was  shown  above  that  the  time  differential  between  the  upstream  and 
downstream  transits  must  be  resolved  to  a  small  fraction  of  a  nanosecond. 
This  can  most  easily  and  reliably  be  accomplished  to  high  accuracy  by 
means  of  a  phase  measurement  between  two  high-frequency  sine  waves 
of  the  same  fixed  frequency. 

The  mathematical  theory  of  operation  of  the  system  is  straightforward 
and  is  given  below. 


The  phase  shifts,  0^  and  0  2*  °*  the  (fixed  frequency)  upstream  and  down¬ 
stream  waves,  respectively,  are  given  below: 


0.  =  2  irf  T. 
r  1  0  1 


0 


0 

1 

(5) 

sin  9 

c  -  v  cos  9 

2  it  f  D 

0 

1 

(6) 

sin  9 

c  +  v  cos  6 

23 

2  irf  T, 
o  2 


where  fQ  =  5  MHz  (the  high-frequency  component).  T,  j  phase  difference 
(01  -  <t>2)  is 
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The  frequency  fQ  is  chosen  so  that  when  v  =  vmax  the  phase  difference 
approaches  but  does  not  exceed  180  . 

The  high-frequency  (fQ)  carrier  wave  is  actually  modulated  by  a  rectangular 
pulse  having  a  width  of  50  microseconds  and  a  repetition  rate  of  2  kHz. 

The  phase  shifts  of  the  fundamental  frequencies  (2  kHz)  of  the  modulation 
waveforms  are  similarly  given  by 
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The  sum  of  these  phase  shifts  (ip  ^  +  ip^)  is 
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We  next  divide  Eq.  (7)  by  Eq.  (10): 
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The  output  of  the  reflectometer  densitometer  is  proportional  to  (pc), 
where  p  =  fluid  density.  Multiplying  Eq.  (11)  by  this  factor  yields 


pc=  fe)  <pvco8e,i(i)“ 


M  cos  0 


where  M  is  the  mass  flow  rate  of  fluid  down  the  pipe.  Note  that  scale 
factors  have  been  omitted  and  M  is  normalized  per  unit  area. 
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Also,  if  we  divide  Eq.  (11)  by  (0  }  +  ip^),  we  obtain 
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Noting  that,  even  under  the  most  extreme  flow  conditions,  (v/c)  <<<  1, 
we  can  simplify  (13): 
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Thus,  in  addition  to  mass  flow  rate,  we  can  also  obtain  an  output  propor¬ 
tional  to  flow  velocity  v.  This  v,  in  turn,  can  be  multiplied  by  p  ,  the 
output  of  a  densitometer  of  more  conventional  design,  which  measures 
p  alone,  not  p  c.  The  product  of  v  and  p  is  M. 

Let  us  next  investigate  the  required  precision  to  which  the  various  quan¬ 
tities  must  be  measured.  Three  "independent"  quantities  are  measured. 
Let  us  consider  the  case  when  these  three  are  (0  j  “  02)»  +  0 2 

In  the  case  of  the  M  determination 
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where  the  magnitudes  of  the  various  relative  errors  are  added  to  take  into 
account  the  worst  case.  Thus,  if  we  wish  to  obtain  a  0.  6%  worst-case 
error  in  M,  the  relative  errors  of  each  of  the  measured  quantities  should 
be  kept  below  ~0.  2%. 

In  the  case  of  the  velocity  determination  (Eq.  (14)), 
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It  should  be  noted  that  we  have  no  frequency  stability  problems  here, 
since  the  fundamental  frequency,  fQ,  is  fixed  and  is  derived  from  a  crys¬ 
tal  oscillator.  The  fundamental  modulating  frequency,  f  ,  is  derived  by 
dividing  down  from  f  =  5  MHz  by  exactly  2500  to  give  exactly  2  kHz. 
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Since  these  frequencies  are  fixed,  narrowband  receivers  can  be  used, 
thereby  providing  extremely  high  signal-to-noise  ratio,  while  operating 
the  transmitter  at  only  a  few  volts  peak-to-peak. 

From  the  above,  it  can  be  seen  that,  by  means  of  independently  (and  un¬ 
ambiguously)  measuring  two  quantities  (A0  and  +  4>2)>  instead  of  the 
single  quantity  (Af)  of  conventional  systems,  the  required  measurement 
precision  can  be  significantly  relaxed,  and  also  the  response  time  of  the 
instrument  can,  if  need  be,  be  made  very  short  without  greatly  sacrific¬ 
ing  measurement  precision. 

The  entire  system  is  digitally  synchronized  to  the  5 -MHz  crystal -controlled 
clock  oscillator,  including  the  digital  phase  measuring  units.  As  will  be 
shown,  the  Tj  -  T£  and  Tj  +  T2  outputs  are  actually  a  series  of  pulses 
which  are  accumulated  over  the  period  of  time  required  for  the  electronic 
computations  of  v/c,  v,  and  M  to  be  made.  Upon  completion  of  a  given  set 
of  computations,  the  stored  pulse  counts  are  fed  to  the  calculator  and  a  new 
set  of  computations  is  begun,  while  at  the  same  time  new  (Tj  -  and 
(Tj  +  T2)  data  is  being  collected  for  the  next  computation.  The  (pc)  data 
is  fed  to  the  calculator  chip  in  BCD  format.  The  digital  computations  are 
presently  being  made  with  the  Mostek  5012  and  5013  calculator  chips. 
Several  seconds  are  required  to  compute  M.  This  time  can  be  decreased 
substantially,  if  necessary,  to  ~1  millisecond  by  using  high-speed  (and 
higher  cost)  digital  computing  elements  or,  with  some  sacrifice  in  com¬ 
puting  accuracy,  by  using  high-quality  analog  components. 

DETAILED  DESCRIPTION  OF  SYSTEM  ELEMENTS 


As  can  be  seen  from  the  foregoing,  the  complete  ultrasonic  mass  flow¬ 
meter  system  consists  of  the  following  subsystems: 

(1)  A  flow  sensor  in  which  the  flow-related  parameters  are  measured. 

(2)  A  densitometer  which  measures  the  density-dependent  acoustic 
impedance,  or,  alternatively,  a  densitometer  which  measures  p  . 

(3)  A  computing  system  which  combines  the  outputs  of  (1)  and  (2)  to 
provide  mass  flow  rate  and  flow  velocity,  as  well  as  intermediate 
quantities  such  as  v/c^  and  v/c. 

The  flow  cell  used  in  most  of  this  program  may  be  denoted  the  (v/c)  times 
p  c  cell,  to  contrast  it  with  the  p  times  v  cell  used  at  the  end  of  the  pro¬ 
gram,  and  which  used  a  dielectric -constant-type  densitometer.  The  (v/c) 
times  (pc)  or  (v/c)  (pc)  cell  included  a  body,  an  inlet,  an  outlet,  two 
transducers  for  the  flow  velocity  measurement,  and  a  densitometer  probe 
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for  the  pc  measurement  (Figs.  2-4).  It  was  fabricated  principally  of 
stainless  steel  304,  for  corrosion  resistance,  and  to  permit  inlet  and 
outlet  transitions  to  be  welded  to  it.  The  rectilinear  channels,  1/2"  x  1" 
for  flow,  and  1/2"  x  1/2"  for  ultrasound,  were  produced  by  electric  dis¬ 
charge  machining.  Transducers  communicate  over  a  3"  path  at  45  to 
the  flow  axis,  essentially  uniformly  interrogating  a  full  cross  section  to 
provide  area-averaging  independent  of  profile.  (Appendix  B  contains 
construction  details,  shows  alternative  designs,  and  includes  background 
information  on  the  need  for  accommodating  various  flow  profiles.  ) 

The  v/c  times  pc  cell  design  principally  differs  from  its  predecessor  of 
the  previous  contract  as  follows.  (1)  Cell  material  is  stainless  steel. 

(2)  Screens  are  soldered  along  sidewalls  to  minimize  eddy  effects  near 
transducer  recesses.  (3)  The  inlet  gradually  transitions  from  circular 
to  rectangular,  to  further  minimize  eddies.  (4)  Sight  ports  are  included, 
with  line-of-sight  being  orthogonal  to  line -of -sound.  (5)  Transducers  are 
backed  by  graphite,  for  thermal  immunity.  (6)  The  pc  densitometer  probe, 
instead  of  being  stepped,  is  nearly  symmetrical  (suggestion  of  N.  E. 
Pedersen)  and  may  be  separated  from  the  "flow  velocimeter"  portion  of 
the  cell.  (7)  Flushing  and  gas  bleed  ports  are  provided  to  allow  entrapped 
air  to  be  released,  and  to  allow  prompt  clearing  out  of  the  transducer 
cavities  and  screens,  should  they  eventually  fill  with  contaminant  particles. 
(8)  Space  is  available  for  installing  a  Kenics  static  mixer  upstream  of  the 
v  measuring  region.  (This  mixer  was  used  in  the  December  tests.  ) 

The  flow  parameter  electronic  measurement  system  provides  two  separate 
outputs.  The  first  is  the  phase  difference  between  the  5-MHz  components 
of  the  upstream  and  the  downstream  waveforms,  while  the  second  output 
is  the  sum  of  the  upstream  and  the  downstream  transit  times.  These 
transit  times  are  determined  by  independently  measuring  the  phase  shifts 
of  the  low-frequency  (2-kHz)  components  of  the  modulation  waveform 
associated  with  the  upstream  and  the  downstream  waves.  These  wave¬ 
forms  are  identical  but  are  shifted  in  time  with  respect  to  each  other. 

The  basic  waveform  is  a  5-MHz  cw  wave,  modulated  by  a  50-ps  rectan¬ 
gular  pulse  having  a  repetition  rate  of  2  kHz  (500  microseconds  between 
pulses).  The  upstream  and  downstream  transmissions  are  shifted  by 
250  (j.  s  but  are  otherwise  identical.  The  transit  time  across  the  fluid 
medium  is  approximately  50  microseconds.  Figure  5  is  a  diagram  show¬ 
ing  the  modulation  waveforms  and  the  timing  sequence  for  upstream  and 
downstream  transmission  and  reception. 

The  transmitted  waveforms  can  be  described  mathematically  as  a  cw 
wave  (the  5  MHz)  modulated  by  a  Fourier  series  representing  the  repet¬ 
itive  pulse: 
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J  igure  3.  View  Through  Transducer  Port.  Screen  is  Soldered  Along  Sidewall  of  Flow 
Channel  to  Reduce  Eddies  in  Measuring  Region.  Screen  is  40-Mesh  Nickel, 
Having  Wire  Diameter  of  0.  1  3  mm  (0.  005  in.  ). 


1974.  Cell  Consists  of  Stainless  Steel  Body,  Velocimeter  Transducer  Fair 
Housed  pc  Densitometer,  Sight  Ports,  Screens,  Inlet  and  Outlet  Transition 
Pipes,  Flushing  and  Vent  System.  Additionally,  in  Tests  Conducted  During 
December,  a  Static  Mixer  Was  Installed  Just  Upstream  of  the  Measuring 
Region,  and  a  Capacitance-Type  Densitometer  Was  Installed  Downstream. 
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Noting  that  the  modulation  waveform  is  positive  or  zero,  depending  upon 
whether  transmission  is  occurring  or  is  not  occurring,  it  is  obvious  that 
Aq  >  0.  This  is  true  because  the  time  integral  of  the  modulation  wave¬ 
form  is  positive,  and  AQ  is  the  only  term  in  the  series  which  does  not 
integrate  to  zero. 

From  the  above  we  see  that  there  will  be  a  5 -MHz  component  in  the  re¬ 
ceived  spectrum  which,  with  appropriate  filtering,  can  be  recovered. 
Also,  if  the  received  wavj  is  amplitude  detected  and  the  resultant  is 
passed  through  a  narrowband  2-kHz  filter,  the  fundamental  frequency  of 
the  modulation  waveform  can  be  recovered. 

If  the  5 -MHz  components  of  the  upstream  and  downstream  waves  are  each 
recovered  as  described  above,  and  are  compared  in  phase,  an  accurate 
measure  of  the  transit  time  difference  between  upstream  and  downstream 
waves  is  obtained.  Also,  the  phases  of  the  2-kHz  fundamentals  of  the  up¬ 
stream  and  downstream  modulations  are  independently  and  unambiguously 
measured  relative  to  a  2-kHz  reference  wave.  These  phases  are  then 
added  to  give  the  sum  of  the  upstream  and  downstream  transit  times. 

Thus,  the  flow  parameter  measurement  system  provides  the  following 
two  outputs:  (T2  -  Tj)  and  (T2  +  Tj).  A  block  diagram  of  this  system 
is  given  in  Fig.  6.  Referring  to  this  figure,  the  system's  operation  will 
now  be  described. 

A  5-MHz  crystal-controlled  oscillator  provides  the  overall  system  syn¬ 
chronization,  as  well  as  the  carrier  wave  to  be  transmitted.  The  output 
of  this  oscillator  is  divided  by  2500  and  serves  to  trigger  two  50  micro¬ 
second  gates,  which  then  modulate  the  5-MHz  cw  waveform.  The  gates 
are  sequentially  operated  250  microseconds  apart,  the  overall  duty  cycle 
being  500  microseconds.  The  switching  logic  operates  analog  solid-state 
switches  which  serve  to  alternately  connect  the  upstream  and  downstream 
transducers  to  their  respective  transmitter  or  receiver  as  required. 

Each  of  the  receivers  employs  a  common  broadband  input  stage,  at  the 
output  of  which  the  signal  is  split  into  two  channels:  the  high-frequency 
and  low-frequency  channels. 

The  operation  of  each  high-frequency  receiver  channel  is  as  follows. 

After  a  stage  of  isolation,  the  broadband  received  waveform  is  fed  to  a 
narrowband  5-MHz  crystal  filter  of  the  bandpass  type.  The  filter 
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Figure  6.  System  Block  Diagram. 


characteristics  are  such  as  to  reject  each  of  the  2-kHz  modulation  spectral 
lines  by  more  than  70  dB.  Thus,  the  filter  output  is  a  5-MHz  cw  wave.  A 
coherent  "local  oscillator"  frequency  of  5.  002*MHz  is  derived  by  means  of 
mixing  the  5-MHz  oscillator  output  with  the  (5-MHz/ 2500)  =  2-kHz  signal, 
and  (by  means  of  a  crystal  filter)  passing  only  the  5.  002-MHz  component. 
This  is  then  mixed  with  the  5-MHz  crystal  filter  output  to  produce  a  2-kHz 
i.  f.  wave  which,  after  bandpassing  in  an  active  analog  filter,  is  converted 
to  a  square  wave  by  means  of  a  zero  cross  detector.  The  phase  shift  of 
this  square  wave  is  equal  to  the  phase  shift  of  the  5-MHz  component  of  the 
waveform  transmitted  across  the  fluid  path. 

The  operation  of  the  2-kHz  portion  of  each  receiver  is  straightforward. 
After  a  stage  of  isolation,  the  received  waveform  is  amplitude  detected 
and  passed  into  a  narrowband  2-kHz  active  analog  bandpass  filter.  The 
filter  output  is  then  fed  to  a  zero  cross  detector,  which  produces  a  square 
wave  at  its  output.  The  phase  shift  of  this  square  wave  corresponds  to 
the  shift  in  phase  of  the  2-kHz  modulation  component  of  the  waveform 
transmitted  across  the  fluid  path. 

The  difference  in  transit  time,  Tj  -  T^,  is,  of  course,  very  Brnall  and 
must  be  resolved  to  a  fraction  of  a  nanosecond.  This  is  done  by  means 
of  measuring  the  time  of  coincidence  of  the  two  5-MHz-derived  2-kHz 
square  waves.  Let  the  time  of  flight  upstream  and  downstream  be  repre¬ 
sented  by  T  +  At  and  T  -  At,  respectively.  A  180  phase  shift  between 
these  square  waves  results  for  At  =  100  nanoseconds.  The  actual  Tj  -  T£ 
measurement  consists  of  counting  the  number  of  5-MHz  pulses  which  occur 
within  the  coincidence  time.  Thus,  each  half  cycle  is  divided  into  2500/2  = 
1250  parts.  Thus,  the  phase  difference  can  be  measured  to  ~0. 14  degree, 
which  results  in  an  instrumental  time  differential  measurement  precision 
of  0.08  nanosecond,'1'  or  80  picoseconds.  Actually,  since,  during  the  mea¬ 
surement  period,  small  flow  fluctuations  result  in  slight  transit  time 
variations,  interpolation  automatically  takes  place,  and  the  resolution  is 
really  better  than  0.  08  nanosecond. 

The  quantity  Tj  +  T£  is  measured  by  counting  the  number  of  5-MHz 
pulses  during  the  coincidence  between  each  modulation-derived  2-kHz 
square  wave  and  a  reference  2-kHz  square  wave,  and  adding  each  count. 
(This  accomplishes  the  time  intervalometer  function.  )  These  pulse  se¬ 
quences,  as  well  as  the  pulse  sequences  from  the  T,  -  T£  measurement, 
are  fed  to  counters  in  the  flow  calculator,  which  will  be  described  in  the 
following  section. 

*As  mentioned  later  in  this  report,  a  measurement  technique  was  devised 
which  cancels  differences  in  propagation  delay  and  phase  shift  variations 
of  the  two  5-MHz  crystal  filters  by  means  of  systematically  reversing  the 
upstream  and  downstream  transducer  connections. 
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Figures  7-9  illustrate  the  flowmeter  electronic  system. 
FLOW  CALCULATOR 


The  purpose  of  this  subsystem  is  to  compute  flow  velocity  v  and  mass  flow 
rate  M,  and  to  provide  a  digital  display  of  these.  Other  quantities  which 
can  be  read  out  are  v/c^,  v/c,  p  c,  and  Tj  +  T2»  where  c  =  sound  speed, 
p  =  fluid  density,  and  Tj  +  =  sum  of  upstream  and  downstream  transit 

times. 

The  calculator  section  makes  the  above  computations  from  previously 
accumulated  data  and,  as  the  calculation  proceeds,  simultaneously  col¬ 
lects  updated  data  for  the  next  computation.  Thus,  both  the  computational 
and  the  storage  functions  operate  continuously,  rather  than  consecutively. 

A  block  diagram  of  the  calculator  section  is  shown  in  Fig.  10.  Four  inputs 
are  provided.  Two  of  these  (the  0  j  -  $  2  anc*  ^1  +  ^2  inPu*8)  accept  the 
pulse  sequences  from  the  flow  parameter  measurement  section.  These 
two  sets  of  pulse  trains  are  counted  by  means  of  a  pair  of  10®  BCD  coun¬ 
ters.  The  data  necessary  for  the  pc  calculation,  namely,  KA  -  LB  and 
KA  +  LB,  are  analog  voltages  and  are  fed  to  two  analog -to- digital  con¬ 
verters.  (As  explained  in  Ref.  1  and  Appendix  C,  A  and  B  are  echo 
amplitudes  generated  with  the  p  c  probe,  and  K  and  L  are  null-balancing 
constants.  )  Upon  completion  of  a  computation,  the  contents  of  the  A/D 
converters  are  fed  to  separate  registers,  and  the  counters  are  reset  (and 
begin  new  counts).  The  register  outputs  are  fed  via  a  data  buss  to  the 
calculator,  which  utilizes  the  Mostek  5012  and  5013  calculator  chips.  The 
control  function  is  provided  by  a  programmable  read-only  memory.  The 
outputs  of  the  calculator  are  selectively  read  out  (bv  means  of  a  multi¬ 
position  switch)  on  a  12-digit  LED  panel  display. 

The  overall  cycle  time  of  the  calculator  is  presently  several  seconds.  As 
mentioned  earlier,  this  time  could  be  shortened  to  about  one  second  by 
means  of  more  efficient  programming,  or  to  less  than  100  ms  through  the 
use  of  one  of  the  new  microcomputer  assemblies.  At  some  reduction  in 
accuracy,  the  computations  could  also  be  performed  by  means  of  integrated 
circuit  analog  computing  elements. 

The  outputs  of  the  calculator  are  presently  unsealed.  We  can  think  of  the 
Avco/ Lycoming  tests  as  the  means  of  experimentally  determining  the 
scale  factor. 

The  flow  calculator  can  also  be  programmed  to  compute  M  from  products 
proportional  or  related  algebraically  to  vc  (for  a  given  fuel)  or  from  v  and 
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Figure  10.  Flow  Calculator. 

Note:  The  KA-LB  and  KA+LB  inputs  are  not  used  in  flow  velocity  mea 
surements. 
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independent  density  data  (M  =  /  v).  Density  is  measured  directly  by  de¬ 
vices  such  as  ITT-Barton's  resonant  vane  densitometer  or  Simmonds 
Precision's  dielectric  constant  densitometer.  Despite  their  intrusive¬ 
ness,  such  devices  appear  to  have  sufficient  resistance  against  clogging, 
and  sufficiently  low  pressure  drop,  as  to  be  potentially  usable  with  the 
ultrasonic  flow  velocimeter  in  some  instances.  To  avoid  the  0.  1  to  <  10 
psi  pressure  drop  that  may  be  associated  with  intrusive  densitometers 
installed  in  series,  one  can  install  the  densitometer  in  a  parallel  branch 
through  which  the  flow  velocity  is  <  1  ft/s,  thereby  assuring  that  the  pres¬ 
sure  drop  will  be  negligibly  small  (Fig.  11).  In  addition  to  the  above  reson¬ 
ance  and  dielectric  types,  nuclear  densitometers  utilizing  one  or  two 
energy  levels  may  also  be  applicable  in  some  cases.  These,  too,  are 
commercially  available  (e.  g.  ,  Tyco).  For  completeness,  these  alter¬ 
natives  to  the  (v/c)  (pc)  formulation  are  indicated  here. 

THERMAL  EFFECTS 

Appendix  E  analyzes  the  effects  of  temperature  on  the  acoustic  properties 
of  the  fuels,  on  the  cell  dimensions  and  integrity,  on  the  transducers  used 
in  the  velocimeter,  on  the  transducers  and  coupling  interfaces  used  in  the 
p  c  densitometer,  and  on  the  fused  silica  buffer  rods  used  in  the  p  c  probe. 
Appendix  E  includes  an  analysis  of  refraction  and  "acoustic  short  circuit" 
(a  term  denoting  the  early  arrival  of  metal-borne  waves  which  interfere 
with  the  measurement  of  fuel-borne  waves).  These  analyses  support  the 
present  choice  of  retaining  the  oblique  fixed-angle  interrogation  path 
formed  by  45  channel  recesses  (but  with  the  addition  of  screens  soldered 
along  the  flow  channel  sidewalls).  One  potential  alternative  to  the  lixed 
incidence  angle  screened  channel  was  a  solid  wedge  installed  in  the  45 
channel  recesses.  However,  its  use  was  discounted  in  the  present  cell, 
because  of  the  refraction  complications. 

Some  of  the  analysis  of  Appendix  E  applies  not  only  to  the  "wetted  trans¬ 
ducer"  concept  used  in  the  present  program,  but  also  to  nonwetted  clamp- 
on  or  externally  mounted  transducers. 

In  a  few  instances,  Appendix  E  includes  test  results  to  indicate  perfor¬ 
mance  of  specific  flow  cell,  transducer  or  p  c  probe  components  at  the 
temperatures  of  liquid  nitrogen  and  boiling  water,  essentially  -200  C  and 
+100  C.  These  particular  temperatures,  while  beyond  the  contractual 
limits  of  -54  C  to  +71  C,  were  selected  for  some  tests,  because  of  the 
simplicity  of  attaining  them  experimentally  and  because  they  illustrate 
performance  beyond  minimum  requirements. 
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Figure  11.  Branch  Concept  for  Densitometry,  (a)  Density  Cell  in 
Series,  Downstream,  (b)  Density  Cell  in  Parallel 
Branch,  Downstream,  (c)  Density  Cell  Branch  in 
Parallel  With  Flow  Velocity  Cell. 
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TEST  RESULTS 


IN-HOUSE  TESTS  OF  FLOWMETER 


In  order  to  determine  that  the  overall  system  was  operating  properly,  it 
was  necessary  to  conduct  a  series  of  s emiquantitative  tests  which  were  to 
be  conducted  prior  to  the  calibration  tests  at  Avco/ Lycoming.  These  pre¬ 
liminary  tests  consisted  of  determination  of  (1)  repeatability,  (2)  zero  drift, 
and  (3)  long-term  stability. 

The  lack  of  an  in-house  calibrated  flow  loop  precluded  a  highly  quantitative 
preliminary  test  program.  However,  it  was  determined  that,  in  the  short 
term,  nearly  identical  flow  conditions  could  be  reproduced  by  consecutively 
valving  off  the  flow  to  zero  and  then  fully  opening  a  gate  valve.  Of  course, 
creation  of  the  zero  flow  condition  is  a  trivial  matter.  Since  we  could  not 
depend  upon  precise  long-term  repeatability  of  a  particular  flow  rate,  the 
long-term  stability  measurement  was  really  a  measurement  of  zero  drift 
over  periods  of  ~l6  hours. 

These  tests  were  most  useful  in  that  they  revealed  two  mechanical  prob¬ 
lems.  The  first  related  to  the  fact  that  the  transducers  were  sealed  to  the 
faces  of  the  mechanical  housing  by  means  of  O-rings.  Changes  in  pressure 
gave  rise  to  slight  shifts  in  the  angular  positions  of  the  transducers,  which, 
in  turn,  caused  slight  changes  in  the  high-frequency  phase  difference.  The 
second  problem  was  more  difficult  to  overcome  in  a  truly  positive  way:  A 
system  of  vent  tubes  (described  elsewhere  in  this  report)  was  installed, 
which  made  it  possible  to  bleed  both  transducer  cavities,  and  also  to  re¬ 
verse  flush  the  cavities.  This  arrangement  made  it  possible  not  only  to 
eliminate  any  air  bubbleb  from  the  cavities,  but  also  to  clean  the  cavities 
and  screens  by  means  of  alternate  forward  and  reverse  flushing. 

After  the  annoying  problems  described  above  were  taken  care  of,  the  re¬ 
peatability  and  short-term  zero  drift  of  the  system  were  consistently  on 
the  order  of  0.  1%  to  0.  2%  of  full  scale.  However,  the  long-term  zero 
drift  of  the  system  was  on  the  order  of  1%,  which  is,  of  course,  unaccept¬ 
able.  Because  of  this  latter  situation,  a  measurement  procedure  was 
developed  which  virtually  negates  zero  drifts  occurring  over  periods 
greater  than  one  minute.*  The  observed  zero  drift  is  due  to  very  slight 
and  nonidentical  thermally  induced  phase  differences  in  the  RF  sections 
of  the  two  receivers.  By  means  of  making  two  consecutive  flow  measure¬ 
ments,  the  latter  with  the  transducer  connections  reversed,  then  adding 


*The  technique  can  be  electronically  implemented  and,  when  this  is  done, 
will  eliminate  zero  drift  for  periods  greater  than  a  few  seconds. 
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the  two  measurements  and  dividing  by  2,  the  effect  is  negated.  The  reason 
is  that  a  relative  phase  lead  for  one  set  of  connections  remains  a  relative 
phase  lead  for  reversed  transducer  connections.  The  flow  induced  relative 
phase,  however,  reverses  when  the  two  transducer  connections  are  re¬ 
versed.  The  procedure  is  mathematically  illustrated  below: 

Suppose  we  have  an  upstream  phase  of  a  j  and  a  downstream  phase  of  <?£> 
and  that  the  phase  shifts  in  receivers  1  and  2  are  (3  j  and  (3  £>  respectively. 
Since  the  system  measures  only  the  magnitude  of  the  phase  shift,  the  first 
measurement,  M^,  is 

=  \(ax  +  pj)  -  (*2  +  P2)|  (18) 

Next,  reversing  the  transducer  connections,  we  perform  a  second  mea¬ 
surement,  M2: 

M2  =  \(az  +  Px)  -  +  P2)|  (19) 

Now,  as  long  as  |a-2  "  ai  I  >  IP  2  '  P 1  I  *  and  a8Buming  flow  in  the  posi¬ 
tive  direction  (e.  g.  ,  >  a  2>» 

M  =  Mj  +  M2  =  («x  -  a2)  +  (pj  -  P2)  +  (ofj  -  o2)  +  (P2  -  pj)  (20) 

J  M  =  («x  -  az)  (21) 

Thus,  as  long  as  the  magnitude  of  the  zero  drift  is  smaller  than  the  ab¬ 
solute  value  of  the  flow  measurement  in  the  absence  of  zero  drift  (~1%  of 
full  scale),  the  zero  drift  is  eliminated  by  this  procedure. 

The  above  technique  was  extensively  and  successfully  tested  in  the 
Panametrics  laboratory  prior  to  the  calibration  tests  at  A vco/ Lycoming. 

CALIBRATION  TESTS  AT  AVCO/ LYCOMING,  JUNE  1974 

The  purpose  of  these  tests  was  to  determine  the  linearity  and  reproducibil¬ 
ity  of  the  volumetric  flowmeter,  and  of  the  overall  mass  flowmeter.  As 
previously  mentioned,  for  purposes  of  analysis  and  diagnostics,  a  multi¬ 
plicity  of  output  readings  was  provided.  These  readings  are  listed  below: 

2 

(a)  v/c  .  This  quantity  is  proportional  to  the  phase  difference  be¬ 
tween  the  5 -MHz  components  of  the  upstream  and  downstream 
waves. 
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(b)  (Tj  +  T2).  This  is  the  sum  of  the  upstream  and  downstream 
transit  times  and  is  derived  from  the  sum  of  the  phases  of  the 
2-kHz  components  of  the  upstream  and  downstream  transmitted 
waveforms.  Note  that  (Tj  +  T2)  is  proportional  to  the  reciprocal 
of  the  sound  speed. 

(c)  v/c.  The  ratio  of  fluid  velocity  to  sound  speed  is  obtained  by 
dividing  (v/c2)  by  (Tj  +  T2). 

(d)  p  c.  The  product  of  fluid  density  and  sound  speed  is  the  acoustic 
impedance,  and  is  the  output  quantity  of  the  pc  densitometer. 

(e)  (v/c)  (pc)  =  M.  This  is  the  mass  flow  rate,  which  is  the 
quantity  of  greatest  interest  in  the  program.* 

2 

(f)  v.  The  fluid  velocity  is  computed  by  means  of  dividing  (v/c  ) 
by  (Tj  +  T2)2. 

As  previously  discussed,  the  various  outputs  were  recorded  in  the  normal 
and  reversed  transducer  connection  conditions,  and  the  average  value  of 
each  reading  was  computed  as  a  part  of  the  data  analysis. 

Flow  measurements  were  made  using  a  Cox  calibration  test  stand  as  the 
M  standard  at  the  Avco/ Lycoming  flow  calibration  laboratory.  The  two 
fuels  used  were  JP-4  and  marine  diesel  fuel,  I  VVF  800.  Measurements 
were  made  on  each  fuel  at  room  temperature  (75°F)  and  at  elevated  tem¬ 
perature  (~  1 1 5  F). 

The  equipment  was  set  up  and  checked  out  on  the  morning  of  June  19. 

JP-4  measurements  were  made  on  the  afternoon  of  June  19  and  on  June  20, 
and  diesel  fuel  measurements  were  made  on  June  21.  In  all,  four  JP-4 
runs  and  two  diesel  fuel  runs  were  made.  Unfortunately,  a  fluid  leakage 
problem  and  a  cracked  fused  silica  buffer  rod  in  the  (p  c)  unit  prevented 
measurement  of  (pc)  in  the  first  three  JP-4  runs. 

The  raw  data  for  each  of  the  six  runs  is  given  in  Appendix  A.  Note  that 
the  "a"  and  "b"  of  each  point  represent  the  normal  and  reversed  trans¬ 
ducer  connection  configurations,  while  "av"  denotes  the  average  of  the 
"a"  and  "b"  measurements.  Note  also  that  the  column  labeled  (v/c)  (pc) 


*Due  to  a  minor  programming  problem  involving  scale  factor,  this  quan¬ 
tity  was  not  directly  read  out,  but  was  computed  subsequent  to  the  data 
runs  by  means  of  directly  multiplying  outputs  (c)  and  (d)  above.  Note  that 
no  loss  of  accuracy  was  incurred  because  of  this. 
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was  subsequently  computed  as  the  product  of  (v/c)  and  (pc),  and  that  M 
is  the  mass  flow  rate  measured  by  Avco. 

Analysis  of  the  data  proceeded  as  follows:  First,  the  data  for  a  given  col¬ 
umn  were  plotted  as  a  function  of  M.  The  best  straight-line  fit  (through 
the  origin)  was  drawn  for  each  plot,  and  the  slope  was  computed.  Then 
the  percentage  of  deviation  of  each  data  point  was  computed  (1)  as  a  per¬ 
centage  of  reading  and  (2)  as  a  percentage  of  full  scale  (which  we  have 
taken  to  be  4000  lb/hr).  These  are  listed  with  subscripts  "lin"  and  "max", 
respectively.  The  results  of  these  computations  are  given  in  Tables  2 
through  7,  in  which  the  average  deviation  was  computed  as  %  of  full  scale 
for  each  quantity  in  each  run.  Any  points  having  deviations  greater  than 
the  statistical  limit  of  error  were  next  excluded  and  the  average  deviation 
recomputed.* 

The  experimental  data,  in  both  tabular  and  graphical  form,  are  presented 
for  all  six  runs  in  Appendix  A.  At  this  point  two  experimentally  observed 
facts  should  be  mentioned: 

First,  it  was  found  that  the  positive  peak  of  the  2-kHz  reference  wave 
(used  to  measure  Tj  and  T2)  was  running  saturated.  Since  this  wave  is 
fed  to  a  zero  cross  detector,  one  might  suppose  that  this  would  be  of  no 
consequence.  However,  since  capacitive  coupling  into  the  zero  cross 
detector  was  used,  a  change  in  amplitude  of  this  reference  wave  would 
result  in  a  slight  shift  of  the  zero  point  as  seen  by  the  zero  cross  detector. 
This  effect  is  additive,  since  we  are  measuring  (Tj  + 

Secondly,  it  was  noted  that  on  two  of  the  data  runs  (#3  and  #5)  an  abnormal 
departure  from  linearity  was  observed  for  mass  flow  rates  below  ~500 
lb/hr,  e.  g.  ,  for  laminar  flow  at  ~10%  of  full  scale,  t  Observation  of  the 
flowing  fuel  through  the  sight  port  revealed  visible,  wavy  striations  in  the 
flow.  The  results  of  this  inhomogeneity  can  be  seen  in  the  graphs  of 
Appendix  A. 


*Such  deviations  were  rare  and  were  generally  found  to  be  isolated  to  only 
one  of  the  several  measured  quantities  for  a  given  M.  Thus,  it  is  reason¬ 
able  to  conclude  that  these  are  reading  errors. 

This  problem  has  been  corrected  since  the  return  of  the  equipment  to 
the  Panametrics  laboratory. 

^Use  of  a  Kenics  "Static  Mixer"  apparently  eliminated  this  nonlinearity 
in  the  December  tests.  Alternatively,  one  could  use  a  flow  cell  dimen¬ 
sioned  for  a  particular  flow  range  to  minimize  departures  from  linearity. 
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It  should  also  be  noted  that  the  range  gates  of  the  two  receivers  were  ad¬ 
justed  between  Run  1  and  Run  2.  This  precludes  legitimate  comparison 
of  the  absolute  slopes  between  Run  1  and  the  subsequent  runs.  However, 
the  internal  consistency  of  the  various  runs  is  not  affected. 

Referring  to  Table  8,  we  see  that  the  average  deviation  (%  of  full  scale) 
of  the  various  measured  quantities  has  a  tendency  to  increase  from  v/c^ 
to  v/c  to  v.  The  reason  for  this  is  attributed  to  the  slight  asymmetric 
saturation  of  the  2-kHz  reference  wave,  as  noted  above.  During  precal¬ 
ibration  tests  at  Panametrics  (when  this  wave  was  not  saturated)  it  was 
noted  that  the  readout  of  (Tj  +  T£)  remained  stable  to  within  better  than 
0.  1%  under  constant  temperature  conditions.  Therefore,  we  would  expect 
that  the  basic  linearities  for  v/c^  would  prevail  for  v/c  and  v  for  an  un¬ 
saturated  reference  wave. 

Note  also  in  Table  8  that  the  deviations  from  linearity  appear  to  be  two 
to  three  times  higher  in  the  case  of  diesel  fuel  than  in  the  case  of  JP-4. 

We  do  not  have  an  explanation  for  this,  but  we  regard  the  basic  linearity 
of  the  v/c^  data  of  the  JP-4  runs  to  represent  the  inherent  precision  of 
the  velocity  dependent  portion  of  the  flowmeter. 

Since  the  density  vs.  temperature  characteristic  of  the  JP-4  and  diesel 
fuel  samples  used  in  these  tests  was  provided  by  Avco,  and  since  one  can 

directly  compare  the  slopes  of  the  j  (pc)  =  ^meaB  vs*  M  curves,  it 

is  possible  to  make  a  number  of  run-to-run  comparisons.  Specifically, 
we  will  take  into  account  temperature -dependent  density  change  and  com¬ 
pare  the  modified  slopes  of  fluid  velocity  vs.  M  for  Runs  2,  3,  4,  5  and  6. 
We  will  also  directly  compare  the  slopes  of  Mmeas  vs.  M  for  Runs  4  and 
5-6.  In  these  comparisons,  we  will  arbitrarily  use  Run  4  as  the  reference 
run. 

The  following  data  are  needed  for  the  comparisons: 

(  p(-80°F)  =  0.  836  gm/cc 
JP-4  J 

*p(  240°F)  =  0.692 

(p(  76°F)  =  0.  830  gm/cc 
Marine  Diesel  I  VVF  800  | 

(  p  (113°F)  =  0.823 
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TABLE  8.  AVERAGE  DEVIATIONS  (%  OF  FULL  SCALE) 


Data  taken  at  Avco/ Lycoming ,  June  19-21,  1974. 

All  data  plotted  vs.  Ivl  (+_  0.  2%)  as  measured  by  Avco. 


Using  the  above  data,  we  find  that  the  temperature  dependences  of  the 
densities  are 


I  J_£_ 

p  0  T 


(  JF-4:  -  0.  06%  Deg  F 

I  Diesels  -  0.  23%  Deg  F 


Applying  the  above  velocity  vs.  M  slopes  as  computed  in  Appendix  A,  we 
have  the  following  direct  comparisons  among  the  runs: 


Run  #  Slope  deviation  (%) 


2 

3 

5 

6 


+4.  5 
-1.  0 
+0.  2 
+1.  1 


Run  4  used  as  standard 


Next,  we  directly  compare  the  slope  of  the  combined  5  and  6  with  Run  4. 
Even  though  the  fuels  are  different,  the  M  curves  should  agree.  Inspec¬ 
tion  of  the  curves  of  (v/c)  (p  c)  vs.  M  for  Runs  4,  5  and  6  (Appendix  A) 
for  M  =  4000  lb/hr  shows: 

Run  4  (JP-4):  (~)  (pc)  =  148 

Runs  5  and  6  (Diesel):  (pc)  =  155 

The  ratio  of  the  above  values  is  0.  955.  Thus,  the  two  runs  agree  to 
within  4.  5%. 

It  is  interesting  to  inspect  the  curves  of  Runs  5  and  6  as  follows:  First, 
referring  to  v/c^  vs.  M,  note  the  large  difference  in  slopes  between  the 
elevated  temperature  (5)  curve  and  the  room  temperature  (6)  curve. 

Next,  inspecting  the  v/c  vs.  M  curves,  we  see  that  the  slopes  are  much 
closer.  Next,  inspecting  the  (v/c)  (pc)  vs.  M  curve,  we  see  that  the 
slopes  coincide  much  more  closely.  This  shows  that  the  system  is  com¬ 
pensating  for  the  temperature  difference  between  the  5  and  6  runs.  Also, 
note  that  the  slopes  essentially  coincide  in  the  case  of  v  vs.  M.  The 
reason  for  this  is  that  the  fractional  sound  speed  variation  of  the  diesel 
fuel  is  -1.  73  x  10"^  deg  F"*,  while  the  density  dependence  on  temperature 
is  only  -2.  29  x  10'4  deg  F"*.  This  fact,  however,  does  not  suggest  that 
the  (pc)  unit  was  not  performing  an  important  correction.  This  correc¬ 
tion  is  directly  seen  in  comparing  the  (v/c)  vs.  M  curves  of  Runs  5  and  6, 
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and  noting  that  the  difference  in  (p  c)  alone  brought  the  two  (v/c)  curves 
into  coincidence.* 

CALIBRATION  TESTS  AT  AVCO/LYCOMING  DECEMBER  1974 

In  this  section,  we  present  the  results  of  a  series  of  tests  run  at  Avco/ 
Lycoming  on  December  17-18,  1974.  The  tests  were  run  on  JP-4  fuel  at 
three  basic  temperatures.  These  were  50°F,  72°F,  and  101°F,  Two 
separate  and  independent  data  runs  were  made  at  each  basic  temperature 
in  order  to  obtain  numerical  estimates  of  the  precision  of  the  flowmeter. 
The  basic  data  were  taken  in  a  manner  identical  to  that  used  during  the 
June  1974  tests. 

The  numerical  data  taken  during  these  tests  are  presented  in  both  numer¬ 
ical  and  graphical  form  in  Appendix  A.  The  analysis  presented  in  this 
section  is  based  entirely  upon  these  numerical  data. 

As  mentioned  in  the  Conclusions,  the  (pc)  data  exhibited  drift  substantially 
in  excess  of  that  required  for  a  final  instrument.  It  was,  therefore,  de¬ 
cided  that  the  density  data  to  be  used  would  be  the  best  available  under  the 
circumstances.  This  was  obtained  with  the  hydrometer  which  is  perman¬ 
ently  installed  in  the  Cox  test  stand  at  Avco.  In  addition,  an  uncalibrated 
Simmonds  densitometer  was  obtained  on  loan  from  Simmonds  Precision 
Instrument  Corporation,  Vergennes,  Vermont,  and  data  was  taken  on  this 
during  all  tests  in  order  to  obtain  a  measure  of  its  stability  and  reproduc¬ 
ibility  as  a  potential  candidate  for  use  in  a  final  prototype  unit. 

As  was  noted  in  the  case  of  the  June  1974  tests,  two  of  the  six  runs  ex¬ 
hibited  a  slight  nonlinearity  at  very  low  flow  rates.  In  these  (December) 
tests,  a  Kenics  "Static  Mixer"  (see  Fig.  B-9)  was  incorporated  imme¬ 
diately  upstream  of  the  flow  cell.  As  can  be  seen  from  inspection  of  the 
data  plots  in  Appendix  A,  no  such  nonlinearities  were  observed  in  the 
December  test  results. 


*Based  on  subsequent  saturation  tests  and  analyses  of  the  densitometer, 
it  appears  that  the  excellent  stability  of  the  echo  comparator  pc  elec¬ 
tronics,  ~0.  1%,  observed  during  these  isothermal  calibration  tests  may 
have  been  obtained  as  a  tradeoff  with  linearity.  Such  a  tradeoff  could 
account  for  the  4.  5%  slope  difference  between  Run  4  and  Runs  5  and  6 
mentioned  on  the  previous  page.  Nonlinear  response  could  also  account 
for  the  observed  difference  of  a  few  percentage  points  between  calculated 
and  measured  p  c  change  in  individual  fuels  tested  at  Avco  during  fuel 
heating  and  cooling  cycles. 
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The  Avco  measurements  were  done  under  the  direct  supervision  of  Mr. 
Joseph  Balcher,  who  participated  in  all  phases  of  the  tests.  The  quantities 
measured  by  Avco  were  mass  flow  rate  M  (lb/hr),  density  p  (gm/cc),  and 
temperature  (deg  F). 

The  Panametrics  measurements  were  readings  of  the  digital  display  on  the 
front  panel  of  the  flowmeter  electronics  package.  The  units  of  these  mea¬ 
surements  are  machine  units,  as  was  the  case  in  the  June  1974  tests.  The 
quantities  measured  were  v/c^,  v/c,  v,  Tj  +  T£»  and  V  ,  where  v  =  fluid 
velocity,  c  =  sound  velocity,  Tj  and  T2  are  the  upstream  and  downstream 
transit  times,  and  Vs  =  Simmonds  densitometer  output  (volts).  As  in  the 
June  tests,  two  readings  were  taken  for  each  test  point,  one  of  these  cor¬ 
responding  to  reversed  transducer  leads.  As  already  shown,  the  average 
of  these  produces  a  result  which  eliminates  residual  instrumental  zero 
offset. 

ANALYSIS  OF  DATA  FROM  FIELD  TEST  #2 


As  mentioned  above,  two  separate  data  runs  were  made  for  each  of  the 
temperatures  (~50  F,  ~70  F,  and  ~100  F).  The  two  data  runs  corres¬ 
ponding  to  a  given  temperature  were  not  run  simultaneously.  If  two  such 
runs  are  taken  together,  included  in  the  scatter  of  data  will  be  such  effects 
as  mechanical  and  thermal  hysteresis  as  well  as  a  potential  long  term 
electronic  drift.  Thus,  the  data  can  be  analyzed  in  three  ways: 

(a)  Linearity.  This  is  expressed  as  the  average  or  standard  deviation, 
expressed  as  %  of  full  scale,  of  the  data  points  of  a  given  run,  from 
the  best  least  squares  straight-line  fit. 

(b)  Precision.  The  average  or  standard  deviation  from  the  least  squares 
linear  fit  of  the  data  points  of  two  separate  data  runs  gives  a  good 
measure  of  the  precision  of  the  flow  measuring  system  for  a  given  T. 

(c)  Temperature  effects.  The  average  or  standard  deviation  from  the 
least  squares  linear  fit  of  the  data  points  of  all  data  points  would  give 
a  measure  of  the  precision,  including  fluid  temperature  effects.  How¬ 
ever,  if  a  definite  trend  in  the  deviation  is  observed,  the  precision 
index  computed  in  this  manner  would  be  artificially  high.  Instead,  if 
such  a  trend  were  observed,  a  temperature-dependent  coefficient 
should  be  introduced  for  the  purpose  of  eliminating  the  effect  of  fluid 
temperature  on  precision. 

(d)  Accuracy.  Since  the  velocity  data  were  taken  as  a  function  of  quanti¬ 
tative  M  and  p  measurements,  the  potential  accuracy  of  the  system 
would  be  essentially  the  same  as  the  precision  discussed  in  (b)  above. 
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One  is  nearly  always  limited  by  the  amount  of  experimental  data  in  the 
specification  of  precision  indices.  In  this  series  of  tests,  sixty-five 
separate  values  of  M  were  run.  Although  we  would,  of  course,  like  to 
see  much  more  data,  we  believe  that  this  number  of  independent  mea¬ 
surements  constitutes  a  good  base  on  which  to  base  estimates  of  linearity, 
precision  and  temperature  dependence. 

ANALYSIS  OF  INDIVIDUAL  DATA  RUNS 


Each  of  the  six  data  runs  was  analyzed  by  means  of  computing  the  average 
deviation*  from  the  best  linear  least  squares  fit  to  the  data.  Since  fluid 
velocity  is  the  basic  quantity  measured  by  the  instrument,  we  have  nor¬ 
malized  the  Avco  measured  M  by  the  product  of  the  Avco  measured  p  and 
the  known  cross-sectional  area  A  =  1/2  in.  2  of  the  flow  cell,  to  give  units 
of  ft/sec.  These  are  given  in  Tables  9  through  14  and  are  the  standards 
against  which  v(exp)  is  measured  in  machine  units.  From  these  two  col¬ 
umns,  the  best  linear  fit  is  obtained.  Then,  the  value  of  velocity  in 
machine  units  is  computed  for  a  perfect  fit.  This  we  denote  by  v(calc). 

The  difference  Av  =  v(exp)  -  v(calc)  is  then  computed.  Next,  |Av|  is 
given  as  a  percent  of  full  scale.  These  results  appear  in  the  right-hand 
column.  Finally,  the  average  deviation  (a.  d.  )  is  computed  for  each  in¬ 
dividual  data  run.  These  values  appear  at  the  bottom  of  each  table.  Note 
that  the  density  (specific  gravity  units)  and  temperature  (  F)  are  given  at 
the  top  of  each  table. 

The  significant  result  of  these  analyses  is  that  the  average  deviations  from 
linearity  are  very  small:  0.15%,  0.10%,  0.10%,  0.14%,  0.  1  3%  and  0.  1 8% 
for  Runs  1  through  6,  respectively. 

ANALYSIS  OF  PAIRS  OF  RUNS 


As  mentioned  earlier,  two  separate  data  runs  were  made  at  each  tempera¬ 
ture.  The  time  sequence  of  these  is  given  below: 

Run  #  Temperature  (°F)  Time 


1 

2 

3 

4 

5 

6 


75 

101 

51 

72 

50 

102 


December  17,  a.  m. 
December  17,  p.  m. 
December  18,  a.  m. 
December  18,  a.  m. 
December  18,  p.  m. 
December  18,  p.  m. 


*Note  that  the  standard  deviation  is  equal  to  the  average  deviation,  mul¬ 
tiplied  by  the  factor  1.  25. 
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.  15% 


^  gj  hOhOhhNhOOO 


Data  taken  at  Avco/ Lycoming,  December  17,  1974  p.  m.  Average  deviation 

0.  10% 


Data  taken  at  A vco/ Lycoming,  December  18,  1974  a.  m.  Average  deviation  = 

0.  10% 
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Data  taken  at  Avco/ Lycoming,  December  18,  1974  p.  m.  Average  deviation 

0.  1 3% 
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Data  taken  at  Avco/ Lycoming,  December  18,  1974  p.  m.  Average  deviation 

0.18% 


Due  to  the  fact  that  no  two  runs  at  a  given  temperature  were  made  con¬ 
secutively,  if  we  analyze  the  composite  results  of  each  set  of  runs  at  the 
same  temperature,  we  can  obtain  a  good  measure  of  the  precision  of  the 
flowmeter.  This  has  been  done  and  the  results  are  given  in  Tables  15,  16 
and  17.  The  method  of  analysis  is  identical  to  that  used  for  the  individual 
runs,  except  that  the  combined  data  of  two  runs  is  used  in  each  analysis. 

The  significant  result  of  these  analyses  is  that  the  average  deviations 
from  the  best  linear  fit  are  0.  25%,  0.  14%  and  0.  21%  for  Runs  3  and  5, 

1  and  4  and  2  and  6,  respectively. 

It  should  be  noted  that,  since  the  Avco  density  measurement  is  subject  to 
a  possible  maximum  error  of  £  0.  2%,  it  is  best  to  take  the  measured 
density  at  the  median  temperature  (75°F),  at  which  p  =  0.  756  gm/cc, 
and  to  then  correct  the  readings  at  50  and  101  if  necessary.  This  can 
be  done  because  the  slope  of  the  JP-4  density  vs.  temperature  charac¬ 
teristic  curve  is  well  known  by  Avco  to  be  -0.  0004  gm^cc  per  deg  F. 

Thus,  when  analyzing  two  or  more  sets  of  data,  it  was  not  necessary  to 
include  the  potential  errors  of  multiple  measurements  of  density.  As  it 
turned  out,  it  was  necessary  to  correct  only  one  of  the  density  readings. 
This  was  the  density  measurement  of  Run  6,  in  which  the  measured  den¬ 
sity  was  0.  746  gm/cc  and  the  corrected  density  was  0.  744  gm/cc.  This 
correction  was  not  necessary  when  analyzing  the  individual  runs  for  lin¬ 
earity,  since  the  measured  density  was  a  constant  for  each  run.  It  was, 
however,  necessary  to  use  the  correction  in  the  composite  analysis  of 
all  six  rims,  which  we  next  describe. 

COMPOSITE  ANALYSIS  OF  THE  SIX  RUNS  . 


Next,  in  order  to  include  the  effects  of  variable  fluid  temperature,  we 
analyzed  the  combination  of  all  six  data  rims.  The  same  analysis  pro¬ 
cedures  were  used  as  in  the  foregoing  cases.  As  mentioned  above, 
density  was  corrected  on  the  basis  of  p  (75  F)  =  0.  756  gm/cc  and 
9  p  /3  T  =  0.  0004  gm/cc  per  deg  F. 

The  significant  results  of  this  analysis  were:  (a)  the  average  deviation 
was  0.  32%  and  (b)  a  definite  temperature -dependent  trend  is  observed. 
Notice  that  the  sign  of  Av,  given  in  Table  18,  is  very  predominantly 
positive  for  the  low  temperature  data,  is  mixed  for  the  room  temperature 
data,  and  xs  very  predominantly  negative  in  the  case  of  the  high  tempera¬ 
ture  data.  This  is,  of  course,  a  very  small  effect,  since  the  average 
deviation  for  the  six-run  composite  analysis  is  only  about  0.  1%  greater 
than  the  mean  of  the  average  deviations  of  the  runs  taken  two  at  a  time. 
However,  this  dependence  can  be  measured  in  more  extensive  tests  over 
a  considerably  wider  temperature  range,  and  the  final  instrument  can  be 
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TABLE  17.  RUNS  #2  AND  #6  -  JP-4  FUEL,  COMPOSITE  ANALYSIS 
p  fy  adjusted  to  0.  744 
T2  =  101  °F 
T  k  =  1 02°F 
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TABLE  18.  RUNS  #1  THROUGH  #6  -  JP-4  FUEL.  COMPOSITE  ANALYSIS 
All  p  's  adjusted  to  p  =  0.  756  for  T  =  75°F 
using  A p  =  0.  0004/deg  F 
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calibrated  so  that  any  such  temperature  dependence  need  not  affect  the 
measurement  accuracy. 

SUMMARY  OF  RESULTS 


The  analytical  results  of  the  December  Avco/Lycoming  tests  are  sum¬ 
marized  in  Table  19.  Although  these  have  been  individually  discussed 
above,  it  is  of  interest  to  inspect  them  in  one  table  and  to  note  how  the 
average  deviations  (a.  d. ' s )  vary  from  run  to  run  and  among  the  compos¬ 
ites.  The  mean  of  the  average  deviations  for  the  individual  run  analyses 
is  0.  13%;  the  mean  of  the  average  deviations  of  the  three  analyses  of  runs 
taken  in  pairs  is  0.  20%;  and  the  average  deviation  for  all  six  runs  taken 
together  is  0.  32%.  It  appears  as  though  the  precision  of  the  velocity 
measuring  instrument  can  be  given  by  an  average  deviation  of  about  +0.  2%, 
which  corresponds  to  a  standard  deviation  of  4^0.  25%.  The  observed  tem¬ 
perature  dependence  of  the  overall  composite  analysis  indicates  that  a 
final  instrument  should  be  calibrated  over  a  reasonably  wide  range  of 
temperatures. 

The  Simmonds  densitometer,  although  uncalibrated,  demonstrated  remark¬ 
able  stability  and  reproducibility.  This  can  be  seen  from  the  Vs  data  of 
Appendix  A.  Using  the  approximate  formula  (given  by  Simmonds) 

p(gm/cc)  =  0.71  (1  +  V/15)  (22) 

s 

where  Vg  =  output  of  densitometer  in  volts,  the  stability  of  the  densitom¬ 
eter  is  exceptionally  good.  Using  the  above  formula,  and  taking  the  ex¬ 
tremes  of  the  Vs  data  for  Runs  3  and  5  combined,  the  variation  in  indicated 
density  is  0.  1%  from  pmax  to  pm^ n.  A  similar  calculation  for  Runs  1  and 
4  gives  a  0.  2%  overall  indicated  density  variation,  while  Runs  2  and  6  give 
an  indicated  density  range  of  0.  1%.  Although  the  density  computed  from 
the  above  (rule  of  thumb)  equation  differs  from  the  corresponding  Avco 
value  by  about  one  percent,  a  calibrated  Simmonds  densitometer  should 
be  considered  as  the  prime  candidate  for  an  alternative  means  of  provid¬ 
ing  fluid  density  for  the  final  mass  flowmeter  system.  Simmonds  has 
made  extensive  measurements  on  JP-4  fuel  over  a  wide  range  of  tempera¬ 
tures  and  mixtures.  As  a  result  of  these  measurements,  Simmonds 
states  that  the  errors  associated  with  their  densitometer  are  conserva¬ 
tively  within  a  2a  value  of  0.  5%,  which  is  equivalent  to  an  average  devia¬ 
tion  of  0.  2%.  Although  the  Simmonds  instrument  might  not  be  used  in  the 
main  fuel  line,  it  could  be  used  in  a  bypass  line,  in  which  only  a  small 


’•‘Discussions  between  N.  Pedersen  and  D.  Stuart  of  Simmonds  Precision 
Instrument  Corporation. 
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fraction  of  the  total  flow  would  take  place*  the  temperature  of  the  fluid  in 
the  bypass  could  be  monitored  and,  if  necessary,  a  small  correction  could 
be  supplied  to  be  utilized  in  the  final  M  computation. 

VIBRATION  TESTS 


Referring  to  Fig.  4,  the  components  considered  most  likely  to  fail  in  a 
vibration  test  were  the  transducers,  the  pc  densitometer  probe  and  the 
soldered  screen. 

A  test  plate  was  designed  and  fabricated,  upon  which  the  above  components 
were  mounted.  Then,  on  July  5,  1974,  a  vibration  test  was  conducted  by 
Associated  Testing  Laboratories,  Inc.  ,  Burlington,  Massachusetts.  The 
test  generally  followed  MIL-E-5009D.  Test  levels  were: 

5  to  14  Hz  0.  5  in.  double  amplitude 
14  to  150  Hz  +_  5  g's 

150  to  500  Hz  +_  20  g's 

Also,  one  15-minute  logarithmic  sweep  (5-500-5  Hz)  was  performed  in 
each  mutually  perpendicular  plane.  No  resonances  were  noted. 

From  these  observations,  it  was  concluded  that  the  items  tested  can  with¬ 
stand  the  vibration  environment  specified. 


*Note  that  the  bypass  line  would  not  be  in  parallel  with  the  flow  cell,  but 
would  preferably  be  downstream  relative  to  the  flow  cell. 
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CONCLUSIONS 


A  mass  flowmeter  system  has  been  designed,  fabricated,  and  tested  for 
fuel  flows  in  excess  of  4000  lb/hr.  Two  separate  test  series  were  run. 

The  first  test  runs  were  made  in  June  1974,  while  the  second  were  made 
in  December  1974.  The  overall  results  of  the  June  Avco/Lycoming  tests 
show  that,  for  JP-4  fuel,  velocimeter  linearity  of  better  than  1%  and  re¬ 
producibility  from  one  day  to  the  next  on  the  order  of  +  1%  of  full  scale 
have  been  achieved.  During  the  tests,  it  was  observed  that  the  (Tj  +  T2) 
portion  of  the  measurement  was  operating  under  nonideal  conditions.  The 
problem  was  subsequently  corrected,  and  it  was  anticipated  that  the  lin¬ 
earity  and  reproducibility  of  the  v  and  v/c  results  would  be  consistently 
within  0.  5%  in  future  tests.  This  was  expected  to  be  confirmed  in  the 
second  series  of  calibration  tests  to  be  conducted. 

The  second  series  of  tests  at  Avco/Lycoming,  run  on  JP-4  in  December 
1974,  did  in  fact  confirm  the  above  predictions.  The  six  individual  veloc¬ 
ity  data  runs  showed  average  deviations  of  less  than  0.  2%  of  full  scale, 
while  the  mean  of  the  average  deviations  was  0.  13%.  The  average  devia¬ 
tion  associated  with  measurement  precision  was  0.  25%  or  less,  while  the 
mean  of  these  average  deviations  was  0.  2%.  The  data  were  taken  over  a 
temperature  range  of  ~50°F  to  ~100°F.  When  all  the  runs  were  analyzed 
together,  it  was  found  that  the  average  deviation  for  the  six  run  composite 
was  0.  32%.  Thus,  one  can  reasonably  conclude  that  a  slight  fuel  tempera¬ 
ture  effect  is  present  and  gives  a  deviation  from  linearity  of  ~0.  2%  per 
100°F.  This  can  be  taken  into  account  in  the  final  instrumental  calibration 
and  should  not  be  considered  to  be  a  measurement  error. 

The  instrumentation  required  to  perform  the  v  and  v/c  measurements  has 
been  demonstrated  to  be  stable  and  reliable  over  many  hundreds  of  hours 
of  operation.  There  is  no  doubt  that  a  reliable  and  accurate  flowmeter 
can  be  fabricated.  Also,  the  area-averaging  technique,  utilizing  a  1"  x 
1/2"  rectangular  duct  with  screened  ports  and  gradual  inlet  transition,  has 
been  demonstrated  to  provide  linear  phase  shift  as  a  function  of  flow  ve¬ 
locity  over  a  50:1  range  of  flow,  including  the  laminar,  transitional,  and 
turbulent  flow  regimes.  Departures  from  linearity  have  not  been  observed 
except  in  some  cases  at  unmixed  flows  less  than  10%  of  full  scale. 

Although  significant  improvements  have  been  made  in  reduction  of  the 
voltage  drift  of  the  pc  densitometer,  random  output  variations  on  the 
order  of  a  few  percent  are  still  observed.  Work  is  still  proceeding  to 
improve  the  stability  of  the  p  c  measurement  to  the  required  0.  5%  or 
better.  At  this  writing,  the  evidence  indicates  that  most  of  the  presently 
observed  output  variations  are  due  to  instabilities  in  the  probe.  Correc¬ 
tive  procedures  are  now  being  implemented,  and  include  a  potentially 
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improved  separated  probe  design.  The  construction  of  a  prototype  mass 
flowmeter,  however,  need  not  be  dependent  upon  successful  operation  of 
the  'c)  (v/c)  =  M  measurement  mode.  A  number  of  commercially  avail¬ 
able  densitometers  exist,  including  the  Simmonds  Precision  capacitance- 
*ype  instrument,  whose  stability  we  measured.  Such  a  densitometer  could 
be  installed  in  the  line  as  indicated  in  Fig.  11,  and  the  product  p  v  could 
be  directly  computed.  At  this  time,  it  appears  that  this  latter  approach 
is  to  be  preferred  over  the  (p  c)  (v/c)  approach. 

We  conclude  that  this  M  =  p  v  approach  could  determine  the  mass  flow  rate 
of  fuels  svch  as  JP-4,  JP-5  or  their  mixtures  with  a  precision  (standard 
deviation)  of  better  than  1%  of  full  scale,  over  the  full  flow  range,  tempera¬ 
ture,  vibration  and  contamination  parameters  of  interest  in  this  program. 
The  measuring  equipment  normally  operates  with  no  moving  parts,  would 
introduce  no  significant  pressure  drop,  and  would  respond  linearly  and 
rapidly  (response  time  ~1  s)  to  mass  flow  rate  in  diagnostic  engine  tests. 
Response  time  could  be  reduced  to  the  millisecond  range  (~10  to  100  ms) 
for  control  purposes,  by  using  high-speed  analog  or  digital  components  in 
the  electronic  portion  of  the  mass  flowmeter  system.  Thus,  not  only 
steady-state  mass  flow  rate,  but  also  dynamic  (transient  or  pulsating)  flow 
could  be  measured. 
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APPENDIX  A 

TEST  DATA  TAKEN  AT  A VCO/ LYCOMING 


In  this  appendix,  we  present  the  data  taken  at  Avco/ Lycoming  during 
two  series  of  Calibration  Tests.  These  were  run  in  June  1974  and  in 
December  1974. 

Since  the  December  results  are  the  more  significant,  these  are  presented 
first.  The  data  presented  in  the  tables  consist  of  points  such  as  2a,  2b, 
and  2av.  This  nomenclature  corresponds  to  normal  and  reversed  trans¬ 
ducer  connections,  and  the  average  of  the  two  readings,  respectively. 

The  following  quantities  were  read  from  the  digital  panel  display  of  the 
flowmeter:  v/c^,  Tj  +  T£,  v/c  and  v.  Numerically,  these  data  are  in 
machine  units.  The  output  of  the  Simmonds  densitometer,  VB,  was  also 
recorded.  The  temperature,  T,  and  mass  flow  rate,  M,  were  measured 
by  the  Avco  team.  The  flow  data,  v/c^,  v/c,  and  v,  are  also  presented 
in  graphical  form  in  the  accompanying  graphs.  In  the  flow  velocity  graphs, 
the  flow  rate  M  has  been  normalized  by  the  factor  p  A  to  give  flow  rate  in 
ft/sec,  where  p  =  Avco  measured  fluid  density  and  A  =  cross-sectional 
area  of  flow  cell  (=1/2  in.  ^). 

Since  the  data  were  taken  at  varying  temperatures,  this  normalization 
permits  easier  comparison  among  the  runs.  For  convenience,  another 
set  of  graphs,  namely,  pv  vs.  M  is  also  provided. 
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TABLE  A-l.  RUN  #1  -  JP-4  FUEL 
DENSITY  p  =  0.  756 


oo 

* 

CO 

O' 

H 

O 

00 

00 

rt 

H 

O' 

sO 

o 

O' 

in 

o 

■** 

M* 

If 

CO 

CM 

CM 

<N3 

m  co  i n  O'  cm  in  oo  (M 

00  O  O'  vO  H  ^  O'  O  00 

h  n  h  t'-  r*-  rr  ro  w 

1*1  rtffl  N  N  N  N  N  N 


00  00  00 
00  00  00 
(*1  1*1  (*1 
in  in  in 


N  (*l  00  U1  N  O'  00  rj*  (*1  N  N 

o  in  m  ^  co  cm  vO  O'  oo  cm  M*  co 

ooo  m  co  co  ooo 

(NJ  N  N  p— *  P-H  p— 4  p— <  i— <  p-H  H  H  H 


(M  O'  O  in  00  NO  O  00  O' 
in  M*  •'f  com  h  m  m  ©  m*  n-  co  •«**  O'  O'  ^  in  n  h  c 

^1,-irH  oovon»  O'  oo  rn  o  o  n  «  n  co  m  m«  incui 

r-  r'-  r-  ^  *f  n  n  n  rH  r-i  O'  O'  o  t'-  r-  h-  m  m  m 

i— i  i— h  >—*  o  o  d  ooo  ooo 


co  if  e 
oo  oo  oo 
m  co  co 
m  m  m 


ooo 
00  00  00 
co  co  co 
in  in  m 


N  CO  N 
00  00  00 
CO  cO  CO 

in  in  in 


cm  o  .M 
00  00  00 
CO  CO  CO 

in  in  in 


O  N  H 
00  00  00 
CO  CO  CO 

in  in  in 


CM  Tf  co 
00  00  00 
CO  CO  CO 

in  in  in 


cm  o  -h  so  ot'-M*  o  M*  n-  co  m  o'  h  m  c-  cm  oo  in  o  co  c 

O  N  00  1*1  O  CO  CM  CO  O  h  CMr-M*  ooom*  o  co  h 

^  cococo  O  O'  O  O  O'  O'  ooo  ooo  ooo  ooo 

>  o'  o'  O'  oo  r-  oo  r-  e  e  e  e  e  m  m  m  cococo 


U]  -*j|  oo  oo  oo 


m  m  m 


*f  c*l  *t 
r-  r'- 


CM  CM  CM 

h-  r-  r- 


i— i  «— i  i— i  o  o'  o  oo  n-  oo 
r-  r>-  n-  r-  MJ  h-  erne 


>  >  >  >  > 
(CtrQflt  (l)  42  flj  ;  l  42  rt  cc)  42  eel  ccj  42  rt 
cm  cm  cm  cococo  M*  M*  M*  in  m  in 


►  > 

(l)  42  Cl)  Cl)  a 

vo  \0  vo  t*-  f-  r- 


86 


O' 

00 

00 

• 

• 

• 

t'- 

rH 

rH 

vO 

rH 

CM 

in 

rH 

o 

in 

CM 

rH 

rH 

rH 

'O 

00 

n- 

CO 

m 

CO 

r~ 

O 

00 

00 

CO 

cO 

00 

vO 

• 

vO 

nO 

• 

CO 

oo’ 

rH 

'£> 

• 

in 

rH 

00 

rH 

O' 

00 

vO 

00 

00 

00 

rH 

o 

o 

n£> 

ro 

CO 

CO 

H 

rH 

rH 

H 

H 

n- 

vO 

(M 

O' 

in 

CM 

'0 

rH 

't 

cm 

CM 

r- 

sO 

lf> 

NO 

o 

p- 

r~ 

rH 

rH 

vO 

in 

CM 

CO 

s£> 

00 

p- 

00 

o 

O' 

O' 

o 

O' 

» t 

vO 

in 

CO 

CO 

CO 

rH 

CM 

rH 

o 

H 

o 

o 

O 

o 

o‘ 

o’ 

o’ 

• 

o 

O 

o 

• 

o 

• 

o 

o’ 

o’ 

o’ 

o 

o 

o 

o 

o 

M< 

r~ 

00 

00 

00 

CM 

CM 

pj 

00 

00 

00 

oo 

O' 

00 

O' 

O' 

O' 

O' 

rH 

ID 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CM 

o 

rH 

m 

m 

in 

in 

in 

in 

m 

in 

in 

in 

m 

in 

rf  PO  00  n-  O'  00  O  nOO 

o  cr  rf  o  ■'f  oo  m  -m 

ooo  o  rH  o  rf  m  in 

N  N  N  i— I  >— I  >— ( 


vO  O'  N 
•  •  • 
oo  n-  m 

T?  O' 
<M  m  M 


CO 

>  O 
> 


'£> 

in 

vO 

M1 

CM 

CO 

o 

O' 

o 

co 

vO 

vO 

vO 

vO 

sO 

'0 

m 

M3 

in 

m 

in 

>  > 
rt  ■  ■  ti  rt  43  rt 
00  00  00  O'  O'  O' 


nl  ^0  nj  (4  rO  n) 

OOO  H  rt  H 


TABLE  A- 2.  RUN  #2  -  JP-4  FUEL 
DENSITY  p  =  0.  744 
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TABLE  A-4.  RUN  #4  -  JP-4  FUEL 
DENSITY  p  =  0.  758 
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TABLE  A-5.  RUN  #5  -  JP-4  FUEL 
DENSITY  p  =  0.  767 
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Figure  A-14.  v  Versus  M/p  A. 
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Figure  A-15.  v  Versus  M/p  A. 
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Figure  A-l6.  v  Versus  M/p  A. 


E 


» 


^ 

\ 

\ 

\ 

\ 


\ 


\ 

\ 

\ 

\ 


O 


-t 

I 

0< 


►"5 


6 

v 


a 


© 

r-* 


o 


o 

in 


< 


a 


D 

C6 


t- 

^  i 

"  rX  <| 

<D 

u 

s> 

00 

•H 

u* 

_  O 

rJ 


r- 

-o 

r- 

o’ 


in 

o 

o 

lTi 


<1  li 


a  H 


SXIMH  AHVHIJH MV  *A 


o 

O 

O 

o 

o 

o 

o 

o 

o 

o 

o 

O 

o 

o 

o 

o 

o 

o 

o 

o 

cq 

X 

o 

N 

X 

"t 

rO 

rO 

rq 

rsl 

rg 

*—» 

127 


Kl/pA  (FT /SEC) 


RUN  6 


c 

oc 


e'¬ 

er 


a 


u 

V 

Q 


o 

NO 


< 

a 


\ 


o 

u 

2 

bo 


#  o 

<M 


nO 

r- 

o’ 

m 

a 


no 

o 


SXINA  AHVHXIGHV  '± 


\ 

\, 


o 

o 

O 

o 

O 

O 

o 

o 

o 

o 

o 

o 

O 

o 

O 

O 

o 

o 

o 

NO 

fO 

00 

o 

ro 

ao 

* 

m 

rO 

CO 

no 

r— ^ 

*•4 

128 


M/o  A  fFT/SEO 


RUN  ) 


o 

o 

o 

in 


f- 

O' 


I  o 

Pi  * 


o 

o 

o 


vO 

in 

Pm 

r-4 

0 

• 

in 

o 

r- 

ii 

ii 

a 

H 

I* 


T 


Slim  AHVHXI9HV 


o 

O 

o 

o 

o 

o 

ro 

M 

129 


& . 


Figure  A-19.  pv  Versus 
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June  1974  Test  Data 
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In  the  following  pages  we  present,  in  tabular  and  graphical  form,  the  com¬ 
plete  test  data  on  six  separate  runs,  which  were  taken  at  Avco/Lycoming 
in  June  1974,  The  tables  are  arranged  in  the  order  in  which  the  data  were 
taken.  The  graphs,  however,  are  ordered,  first,  according  to  the  mea¬ 
sured  quantity,  and  secondly,  according  to  the  run  number.  This  is  done 
for  easy  comparison  of  a  given  measured  quantity  (e.  g.  ,  v/c^,  v,  etc.  ) 
from  run  to  run. 
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Figure  A-34.  v/c  Versus 


Figure  A-35.  v  Versus  M. 


Figure  A-36.  v  Versus  M* 
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APPENDIX  B 

AREA-AVERAGING  ULTRASONIC  FLOWMETERS* 


BACKGROUND 


This  appendix  relates  in  general  to  ultrasonic  flowmeters  and  more  par¬ 
ticularly  to  area-average  measurement  of  fluid  flow  through  conduits. 

In  the  measurement  of  the  fluid  velocity  flow  rates  where  the  fluids  are 
confined  within  a  conduit,  one  nonintrusive  technique  which  has  proved 
to  be  successful  in  a  number  of  applications  is  the  use  of  electroacoustic 
transducers  to  generate  ultrasonic  waves  which  are  transmitted  diago¬ 
nally  across  the  conduit,  often  at  45  ,  such  that  a  significant  component 
of  the  ultrasonic  beam  is  transmitted  in  the  direction  parallel  to  the  di¬ 
rection  of  flow  of  the  fluid  through  the  conduit.  The  ultrasonic  modes 
generated  at  the  transducers  may  be  either  longitudinal  waves  or  shear 
waves,  with  the  latter  type  being  converted  by  refraction  at  a  solid/fluid 
interface  into  longitudinal  waves.  The  measurement  of  the  flow  velocity 
is  made,  usually,  by  electronic  circuits  which  respond  to  the  small  dif¬ 
ference  in  transit  time  between  a  longitudinal  wave  traveling  diagonally 
downstream  and  one  traveling  diagonally  upstream.  The  measurements 
may  be  made  with  waves  traveling  different  paths  or,  preferably,  over 
a  common  path.  Also,  the  transit  time  measurements  may  be  made 
sequentially  on  upstream  and  downstream  transits,  or,  by  employment 
of  pairs  of  transducers,  they  may  be  made  simultaneously.  The  deter¬ 
mination  of  transit  time  may  be  accomplished  by  clocking  the  time  dif¬ 
ference  between  the  leading  edges  of  two  video  pulses  or  alternatively  by 
measuring  the  phase  shift  between  emitted  continuous  waves  or  bursts  of 
continuous  waves.  Another  technique  for  flow  velocity  measurement 
using  ultrasonic  beams  involves  measurement  of  doppler  frequency  shifts 
of  scattered  waves.  Doppler  methods,  however,  are  presently  restricted 
mainly  to  blood  flow  velocity  measurements. 

In  situations  where  it  is  required  to  nonintrusively  measure  the  velocity 
of  flow  of  fluids  in  conduits,  the  specific  measurement  conditions  may 
vary  rather  widely.  For  example,  the  area  over  which  measurements 
are  to  be  made  may  vary  from  less  than  one  square  millimeter  in  capil¬ 
lary  tubes  to  several  square  meters  in  water  or  sewage  pipes.  The  time 
during  which  a  measurement  is  to  be  made  may  vary  from  a  few  milli¬ 
seconds  in  engine  control  applications  to  one  minute  in  chemical  process 
plants,  or  even  one  whole  day  for  power  plant  exhaust  stack  monitoring. 
The  precision  or  accuracy  requirements  are  usually  within  one  percent 
and  may  be  as  low  as  one  tenth  of  a  percent.  The  variation  in  flow  ve¬ 
locity  encountered  also  extends  over  a  wide  range  from  theoretically  zero 

v  The  content  and  language  of  this  appendix  are  based  on  a  U.  S.  patent 
application,  which  ha6  been  allowed  and  is  expected  to  issue  in  1975.18 
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velocity  at  the  walls  of  the  conduit  to  speeds  approaching  or  equal  to  the 
speed  of  sound,  or  even  faster,  in  supersonic  or  hypersonic  flow  as  en¬ 
countered  in  wind  tunnels.  Typical  liquid  flows  seldom  exceed  3  m/s, 
but  velocities  of  ~30  m/s  are  sometimes  encountered. 

In  the  flow  of  fluids  in  closed  conduits,  it  is  clear  that  the  flow  velocity 
will  vary  across  the  area  of  the  fluid,  usually  being  lowest  near  the  walls 
of  the  conduit  and  increasing  to  a  maximum  velocity  somewhere  near  the 
center.  The  profile  of  these  flows  will  vary  widely  from  case  to  case, 
and  may  include  asymmetries,  swirl,  and  reverse  flow.  Even  in  a  given 
case,  the  profile  is  generally  different  for  the  static  or  dc  component  of 
flow  as  compared  with  the  dynamic  or  ac  component.  For  example,  if 
one  analyzes  small  100-Hz  oscillations  superimposed  upon  a  larger  steady 
flow,  the  profile,  theoretically,  can  be  significantly  flatter  for  the  100-Hz 
portion  of  the  flow.  Interrogation  of  the  flowing  fluid  with  an  ultrasonic 
wave  that  averages  over  the  diameter  of  the  conduit  but  not  over  the  whole 
area  of  the  conduit,  provides  a  measurement  of  flow  velocity  v^  equal  to 
va/K,  where  v<j  is  the  diameter  averaged  flow  velocity,  va  is  the  area- 
averaged  flow  velocity  and  K  is  a  constant  which  is  less  than  unity,  the 
actual  value  of  which  depends  upon  the  flow  profile.  In  special  cases,  such 
as  laminar  flow  or  fully  developed  turbulent  pipe  flow,  K  can  be  calculated 
in  terms  of  the  Reynolds  number.  The  Reynolds  number  expresses  the 
ratio  of  inertial  to  viscous  forces.  Numerically,  this  dimensionless  ratio 
is  given  by  the  expression  Re  =  p  VD/n»  where  p  is  fluid  density,  V  is  free- 
stream  flow  velocity,  D  is  channel  width  or  diameter,  and  rj  is  viscosity 
coefficient.  For  laminar  flow  of  Newtonian  liquids,  where  the  Reynolds 
number  is  typically  less  than  two  thousand,  the  profile  across  the  diam¬ 
eter  is  essentially  parabolic  in  a  smooth,  regular,  circular  pipe.  Under 
these  conditions  an  ultrasonic  flow  measurement  over  the  diameter  pro¬ 
vides  a  value  for  flow  velocity-v^  "which  is  larger  than  the  true  area- 
average"  value  v&  by  33%  (K  =  0,  75,  1/K  =  1.  33).  Accordingly,  where 
the  specific  flow  profile  through  the  conduit  is  known,  the  values  of  the 
coefficient  K  may  be  utilized  to  convert  a  measurement  of  velocity  across 
the  diameter  to  the  area-averaged  velocity.  However,  in  a  number  of 
circumstances,  the  flow  profile  is  ?t  best  uncertain  due  to  non-Newtonian 
liquids,  or  due  to  the  proximity  of  oends,  inlets,  outlets,  side  ports  or 
other  possibly  unknown  perturbations  or  irregularities  within  the  conduit. 
Thus,  the  reliability  of  the  output  indication,  in  terms  of  average-area 
velocity,  can  deteriorate  rapidly.  Even  in  those  circumstances  where 
thj  fluid  type  is  known,  e.  g.  ,  JP-4,  diesel  fuel,  etc.  ,  the  value  of  K  may 
not  be  known,  particularly  when  temperature  gradients  are  present,  be¬ 
cause  it  depends  upon  the  value  of  Reynolds  number  Re,  which  can  itself 
vary  when  wide  temperature  ranges  are  encountered. 
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The  above  difficulties  have  been  recognized  for  some  twenty  years.  It 
is  not  surprising,  therefore,  to  find  that  in  order  to  obtain  sufficient  ac¬ 
curacy  of  ultrasonic  measurements  under  these  uncertain  profile  condi¬ 
tions,  several  approaches  have  been  proposed  by  earlier  investigators. 

Kritz  (1955)  published  nomograms  relating  va  to  v^  in  terms  of  Reynolds 
number,  Fig.  B-l.  A  graph  containing  this  information,  due  to  McShane 
(1974),  is  shown  in  Fig.  B-2.  Swengel  (1956)  originated  an  intrusive 
technique  using  a  number  of  distributed  points,  e.  g.  ,  using  slender  rods, 
or  using  an  array  of  transducers,  Fig.  B-3.  But  Swengel's  approach 
introduces  three  sources  of  error:  (1;  probes  perturb  the  flow;  (2)  re¬ 
flections  lead  to  multipath  ambiguities;  (3)  resonant  rods  or  separate 
transducer  elements  are  not  suitable  or  likely  to  provide  a  wave  of  uni¬ 
form  intensity,  as  required  for  uniform  weighting.  Petermann  (1959), 
although  recognizing  the  need  for  area-averaging,  did  not  describe  an 
adequate  area-weighting  procedure,  despite  his  use  of  multiple  bounces 
down  the  pipe  axis  (Fig.  B-4)  and  analysis  of  refractive  focussing. 

Fisher  and  Spink  (1972)  described  a  multichord  Gaussian  approach  which, 
for  pipe  diameters  of  about  0.  2  m  or  larger,  apparently  has  provided  the 
best  accuracy  to  date,  about  0.  1%  (Fig.  B-5)  when  4  chords  ace  used. 
Suzuki  et  al  (1972)  have  reported  on  measurements  in  orthogonal  planes, 
which  apparently  goes  part  way  toward  obtaining  a  better  area-average 
than  a  single-path  measurement  (Fig.  B-6),  at  least  for  large  pipes,  with 
diameters  from  ~0.  2  m  to  over  1  m.  See  Refs.  3-7. 

For  area-averaging  over  a  limited  range  of  flow  conditions,  Pfau  (1970) 
derived  a  preferred  location,  at  the  intersection  of  two  profiles,  where 
doppler  measurements  should  be  made.®  Range-gated  doppler  measure¬ 
ments,  to  map  out  the  profile,  or  at  least  segment  it  into,  say,  eight 
zo~es  (Figs.  B-7  and  B-8),  were  described  by  Haase,  Foletta  and  Meindl 
(1973). ^  Doppler  measurements,  however,  are  seldom  used  in  industrial 
flowmeters,  probably  due  to  the  variability  of  scatterers  and  sound  speed. 

In  USAAMRDL  TR  72-66,*  p.  20,  we  suggested  using  a  Kenics  static  mixer 
to  establish  an  essentially  flat  profile,  to  provide,  in  a  limited  region,  a 
situation  where  essentially  K  =  1  (Fig.  B-9).  (The  second  series  of  cali¬ 
bration  tests  in  the  present  program  used  this  type  mixer.  )  On  p.  64, 

Fig.  28  of  TR  72-66,  we  showed  an  offset  configuration  similar  to  designs 
recently  sold  by  at  least  one  ultrasonic  flowmeter  manufacturer.  The 
offset  provides  area-averaging  but  suffers  from  inlet  and  outlet  perturba¬ 
tions  and  undesirable  pressure  drops.  An  earlier  offset  design  was  used 
by  Noble  (1968).10 

Figure  B-10  illustrates  a  noviultra<sonic  approach  to  the  flow  profile 
problem,  the  "Annubar,  "  an  intrusive  area-averaging  pitot-tube  flowmeter 
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NUMBER 


Figure  B-l.  Percent  Deviation  From  Calculated  Calibration  Versus 
Reynolds  Number.  After  Kritz  (1955). 


Figure  B-2.  Ratio  of  Measured  Velocity  to  Average  Velocity  in  a  Round 
Pipe  Versus  Reynolds  Number.  After  McShane  (1974). 
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Figure  B-3.  Perspective  View  Shows  Use  of  Slender  Rods  Installed 
Intrusively  Within  Duct,  Intended  to  Provide  Area- 
Average.  End  View  Shows  Arrays  of  Piezoelectric 
Sources  (110  and  112)  Intrusively  Installed  Within 
Duct  (108),  Intended  for  Same  Area-Average  Purpose. 
After  Swengel  (1956). 


Figure  B-4. 


Multiple  Bounces  Down  Pipe  Sample  Fluid  Volume  Over 
Considerable  Axial  Distance,  But  Do  Not  Adequately 
Weight  the  Radial  Distribution  of  Flow  Velocity.  Multiple 
Bounce  and  Transducer  Configurations  After  Petermann 
(1959). 
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Figure  B- 


Figure  B- 


5.  Measurement  Plane  in  Round  Pipe,  Illustrating  Gaussian 
Quadrature  4-Chord  Technique.  After  Fisher  and  Spink 
(1972).  Chord  Locations  Corresponding  to  Numerical 
Integration  or  Quadrature  Methods  of  Gauss,  Chebycheff 
and  Lobatto  are  Given  by  Malone  and  Whirlow  (1971).  ^2 


6.  Two-Channel  System 

D:  Downstream  Direction 
U:  Upstream  Direction 

This  figure  shows  the  time  chart  of  an  installation  at  a  power 
station.  Two  pairs  of  probes  are  mounted  where  the  sound 
path  for  each  pair  crosses  at  right  angles.  The  flow  rate  is 
measured  from  the  sum  of  two  differential  sing-around  fre¬ 
quencies  on  two  paths.  The  effect  of  slant  or  circular  move¬ 
ment  of  flow  proved  to  be  largely  improved  by  such  a  method 
of  mounting.  After  Suzuki  et  al  (1972). 
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Figure  B-7.  Perspective  View  of  Velocity  Profile  During  One  Cardiac 
Cycle  in  a  Canine  Ascending  Aorta.  After  Haase,  Foletta 
and  Meindl  (1973). 


Figure  B-8.  Time-Average  Velocity  Profile  in  Human  Ascending  Aorta. 
After  Haase,  Foletta  and  Meindl  (1974). 
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Element  Number 


#  0 

4  8  16 

Number  of  Striations 


32 


Clockwise  Rotation 
Right-Hand  Element 


Counter  Clockwise  Rotation 
Interface  Left-Hand  Element 


The  fixed  geometric  design  of  the  unit  produces  unique 
patterns  of  flow  division  and  radial  mixing  simultaneously. 

In  laminar  flow,  a  processed  material  divides  at  the  leading 
edge  of  each  element  and  follows  the  channels  created  by 
the  element  shape.  At  each  succeeding  element  the  two 
channels  are  further  divided,  resulting  in  an  exponential 
increase  in  stratification.  The  number  of  striations  pro¬ 
duced  is  2n  where  n  is  the  number  of  elements.  In  either 
laminar  or  turbulent  flow,  rotational  circulation  of  a  pro¬ 
cessed  material  around  its  own  hydraulic  center  in  each 
channel  of  the  mixer  causes  radial  mixing  of  the  material. 

All  processed  material  is  continuously  and  completely 
intermixed,  resulting  in  virtual  elimination  of  radial 
gradients  in  temperature,  velocity  and  material  composition. 
Illustrations  courtesy  Kenics  Corporation,  Lawrence,  Mass. 

The  Kenics  Static  Mixer  used  in  the  December  1974  tests  is 
identified  as  follows: 

Serial  No.  1128  Model  No.  5-10-321-00 

Maximum  design  pressure.  Maximum  temperature, 

3800  psi  300°F 

Dimensions: 

.  835  O.  D.  x  .  625  I.  D.  x  6"  long 

Figure  h  9.  Kenics  "Static  Mixer".  This  Unit  is  a  Series  of  Fixed  Helical 
Elements  Enclosed  Within  a  Tubular  Housing. 
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Figure  B-10.  Example  of  Flowmeter  Responsive  to  Flow  Profile. 

Each  annular  segment  of  the  "Annubar"  primary  flow 
element  represents  an  area  of  equal  cross-sectional 
size.  Accurate  detection  of  flow  within  these  annular 
segments  is  achieved  by  large  sensing  ports.  The 
various  flow  rates  are  averaged  by  the  plenum  of  the 
upstream  element.  Source:  Ellison  Instrument 
Division,  Dieterich  Standard  Corporation,  Boulder, 
Colorado. 
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available  from  Ellison  Instrument  Division,  Dieterich  Standard  Corpora¬ 
tion,  Boulder,  Colorado  80302. 

To  further  emphasize  that  the  profile  problem  is  generic  to  all  fluid 
flowmeters,  we  may  refer  to  the  analyses  of  electromagnetic  flowmeters 
by  Bevir  (1970)  and  by  Haacke  (1974).  11»  12 

Turtle  (1974)  introduced  a  converging  throat  section,  Fig.  B-ll,  to 
"normalize"  the  velocity  profile,  but  at  the  expense  of  the  pressure  drop 
associated  with  the  reduced  cross  section.  13 

Despite  the  variety  of  the  foregoing  approaches,  accurate  measurement 
of  area-averaged  flow  velocity  in  small  conduits  (diameter  <  1  to  3  cm) 
of  the  type  of  concern  in  this  program,  over  the  full  range  of  flow  pro¬ 
files  to  be  encountered,  cannot  be  obtained  by  the  methods  of  the  other 
investigators  cited  above. 

On  the  other  hand,  our  area-averaging  flow  cell  and  transducer  design 
given  in  USAAMRDL  TR  72-66,  particularly  with  the  inclusion  of  im¬ 
provements  such  as  screens  over  the  45  transducer  portB  to  minimize 
eddies  in  the  measuring  path,  and  a  gradual  transition  at  the  inlet,  pro¬ 
vides  an  accurate  area-average.  Calibration  data  given  elsewhere  in 
this  report  show  that  the  device  is  reproducible  and  essentially  linear 
over  wide  ranges  of  Reynolds  number.  Since  this  desirable  performance 
has  now  been  demonstrated,  it  is  appropriate  to  elaborate  on  the  use  of 
a  rectilinear  flow  channel,  in  which  all  the  fluid  passing  through  is  mea¬ 
sured  by  an  obliquely  incident  rectilinear  beam,  with  equal  area  elements 
of  flow  being  given  equal  weight  by  the  uniform  beam. 

SUMMARY  OF  THE  AREA -AVERAGING  CONFIGURATIONS 


Broadly  speaking,  in  the  present  approach  the  correct  weighting  of  the 
flow  profile  in  the  measurement  of  flow  velocity  across  a  conduit  is  ac¬ 
complished  by  flowing  the  material  through  a  section  of  conduit  having  a 
rectilinear  cross  section  in  a  plane  transverse  to  the  direction  of  flow 
and  interrogating  the  fluid  with  a  substantially  plane  ultrasonic  beam, 
itself  having  a  rectilinear  cross  section,  with  the  beam  dimension  and 
the  conduit  dimension  along  one  coordinate  being  equal  and  with  the  beam 
dimension  in  the  other  direction  along  the  other  coordinate  being  equal  to 
or  less  than  the  conduit  dimension  along  the  same  coordinate.  The  mea¬ 
surements  are  preferably  made  in  a  pathlength  of  fluid  within  the  very 
near  field  of  the  transducer  generating  the  ultrasonic  beam  so  that  the 
curvature  of  the  interrogating  wave  by  the  flow  of  the  fluid  is  minimized. 
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Figure  B-ll.  Converging  Throat  Section  "Normalizes"  Flow.  After  Turtle  (1974) 


In  general,  the  pathlength  of  the  ultrasonic  beam  in  the  flowing  fluid 
should  not  exceed  d^/X.  ,  where  d  is  the  smaller  dimension  of  the  inter¬ 
rogating  beam  and  X  is  the  wavelength  in  the  fluid.  Thus,  in  a  typical 
situation  where  X  =  1  mm,  if  the  diameter  d  equals  10  mm,  then  the 
pathlength  should  not  exceed  approximately  100  mm.  Generally,  for 
preferred  operating  conditions,  the  pathlength  should  not  exceed  one 
quarter  of  this  distance,  i.  e.  ,  25  mm  long.  In  many  situations,  the 
rectilinear  conduit,  which  may  be  either  square  or  rectangular,  is  a 
flow  cell  introduced  into  the  flow  path.  The  flow  cell  assembly,  then, 
is  an  entire  apparatus  including  ultrasonic  transducers  mounted  to  pro¬ 
duce  the  appropriate  interrogating  ultrasonic  beams.  These  beams  which 
bidirectionally  interrogate  the  fluid  in  response  to  transducer  energiza¬ 
tion  may  be  given  their  rectilinear  shape  in  any  one  of  several  ways.  For 
example,  piezoelectric  generators  and  detectors  may  be  ground  to  the 
desired  shape,  or  a  circular  disk  may  be  employed  with  at  least  one  of 
its  surfaces  electroded  in  a  rectilinear  pattern.  Yet  another  approach 
is  to  collimate  generated  beams  of  any  shape  by  means  of  opposed  cav¬ 
ities  or  slots  having  rectilinear  cross  sections. 

By  employing  beams  of  rectilinear  shape  to  interrogate  the  fluid  in  con¬ 
duits  also  of  rectilinear  shape,  all  of  the  segments  between  the  ceiling 
and  the  floor  of  the  flow  channel  may  be  measured,  and  correct  area¬ 
averaging  may  be  provided. 

DESCRIPTION  OF  THE  DRAWINGS 


Figure  B-I2a  is  an  illustration  in  plan  view  of  a  flowmeter  cell  con¬ 
structed  in  accordance  with  the  principles  of  this  invention.  Figure 
B-12b  is  a  cross-sectional  view  along  the  line  bb  of  Fig.  B-l.  Figure 
B-12c  is  a  cross-sectional  view  taken  along  the  line  cc  of  Fig.  3-la. 
Figures  B-13-B-19  are  cross-sectional  and  end  views  of  portions  of 
alternative  embodiments  of  flowmeter  cells  constructed  in  accordance 
with  the  principles  of  this  invention.  Figure  B-20a  is  an  illustration  of 
the  probe  of  Fig.  B-20  taken  along  the  line  ;ia.  Figure  B-21  is  an  illus¬ 
tration  in  cross-sectional  view  of  another  configuration  of  a  flowmeter 
cell  constructed  in  accordance  with  the  principles  of  this  invention. 
Figures  B-22  and  B-23  are  illustrations  in  a  cross-sectional  view  of 
additional  embodiments  of  flowmeters  constructed  in  accordance  with 
the  principles  of  this  invention.  Figures  B-24  and  B-25  are  illustrations 
in  a  cross-sectional  view  of  portions  of  flowmeter  cells  constructed  in 
accordance  with  the  principles  of  this  invention,  illustrating  specific 
details  thereof.  Figure  B-26  illustrates  cross -polarized  shear  wave 
transducers. 
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Figure  B-12a.  Area-Averaging  Flow  Cell>  Top  View. 
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Figure  B-12c.  Cell  End 
View. 
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Figure  B-21.  Cavity  Reflection  Controls  Oblique  Interrogation. 


Figure  B-22.  Vee-Grooved  Walls,  Used  Nonrefractively. 


Figure  B-23.  Vee-Grooved  Walls,  Operating  Refractively. 
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Figure  B-24.  Screen  Over  Transducer  Cavity 
Reduces  Turbulence. 


o 


Figure  B-25.  Plate  Installed  in  Cavity  to 
Reduce  Turbulence. 


Figure  B-26.  Cross -Polarized  Shear 

Wave  Transducer  Assembly. 


DESCRIPTION  OF  PREFERRED  EMBODIMENTS 


With  reference  to  Figs.  B-l2a,  B-12b  and  B-I2c,  there  is  shown  an 
embodiment  of  an  ultrasonic  flowmeter  cell.  The  cell  body  4  is  formed 
o\it  of  a  suitable  material,  such  as  stainless  steel  or  aluminum.  Other 
suitable  materials  include  plastic,  glass  or  many  other  solids.  The  cell 
contains  an  axial  flow  channel  or  conduit  with  a  rectilinear  cross  sec¬ 
tion  having  a  width  dimension  W  and  a  height  dimension  H.  The  conduit  1 
typically  has  a  length  two  to  ten  times  longer  than  the  larger  dimension  of 
its  cross  section.  The  inlets  and  outlets  to  the  flow  channel  1  are  shown 
in  Fig.  B-12b  as  threaded  for  connection  to  the  line  through  which  the 
fluid  to  be  measured  is  flowing.  Tapered  fittings  or  flanged  connections 
may  be  substituted  for  the  threaded  inlets  and  outlets  shown,  gradual 
transitions  being  preferred.  Welded,  brazed,  soldered  or  epoxied  fittings 
may  also  be  used. 

A  pair  of  ultrasonic  transducers  2  are  positioned  within  a  pair  of  recesses 
3,  with  the  recesses  extending  in  a  direction  oblique  to  the  axis  of  flow  in 
the  conduit  1.  The  transducers  2  generate  a  beam  of  ultrasonic  waves 
which  is  characterized  by  a  rectilinear  cross  section  in  the  plane  normal 
to  the  axis  of  propagation  of  the  beam.  The  cross  section  of  the  beam  has 
a  dimension  W  by  kH,  where  k  has  a  value  typically  between  0.  5  and  1.  0. 
The  axes  of  the  recesses  3  are  parallel  to  one  another.  The  transducers 
2  are  arranged  to  launch  an  oblique  wave  in  the  cell  body  material  at  an 
angle  9  0  with  respect  to  the  normal  drawn  to  the  axis  of  the  conduit  1. 

Transducers  2  may  be  made  of  any  conventional  piezoelectric  material 
such  as  lead-zirconate-titanate,  lead  metaniobate,  etc.  The  oblique 
waves  transmitted  through  the  cell  body  at  angle  9  interrogate  the  fluid 
flowing  through  conduit  1  at  an  angle  9  j,  where  the  angles  9  Q  and  9  j  are 
related  to  the  sound  speeds  cQ  and  Cj  in  the  solid  and  fluid,  respectively, 
in  accordance  with  Snell's  Law:  cQ/sin  9  Q  =  c^/sin  9  j.  The  position  of 
transducer  2  with  respect  to  conduit  1  will  depend,  therefore,  on  the 
angles  9  Q  and  9  ^  and  also  on  the  number  of  times  the  ultrasonic  beam 
is  to  transverse  the  fluid  path.  Thus,  while  in  Fig.  B-12b  the  beam  is 
indicated  as  having  a  single  traverse  of  the  fluid  path,  by  arranging  for 
the  greater  separation  between  the  two  transducers  along  the  axis  of 
flow,  a  multiple  traverse  may  be  utilized,  to  increase  the  time  differ¬ 
ences  for  upstream  and  downstream  waves. 

The  conduit  within  the  cell  4  may  be  formed  by  being  broached  or  electric- 
discharge -machined  between  the  inlet  and  outlet  ports.  Alternatively,  the 
conduit  may  be  milled  out  and  the  coverplate  5  (as  shown  in  dotted  lines 
in  Fig.  B-I2c)  may  then  be  attached.  The  coverplate  5  may  be  either 
permanently  attached,  such  as  by  welding,  or  may  be  removably  attached, 


162 


*  «■  VIKMKl  W  *»*  »■ .  m<- 


^^'^*BIWW’IW'«T  •  ^t’W^SWWB'fn 


for  example,  by  strews  or  clamps.  The  coverplate  5  may  be  formed  of 
a  transparent  material  such  as  a  clear  plastic  so  that  the  flow  pattern  and 
clogging  or  buildup  of  residual  materials  may  be  visually  observed  with¬ 
out  the  distortion  of  a  curved  lens. 

The  orthogonal  walls  of  the  conduit  1  enable  other  nonintrusive  ultrasonic 
measurements  on  either  a  pulse-echo  or  through-transmission  mode  to  be 
made.  These  measurements  may  include  time  interval  measurements 
across  the  fluid  perpendicular  to  the  direction  of  flow  to  determine  the 
sound  speed,  cj,  of  the  fluid.  Additionally,  the  reflection  coefficient 
at  one  of  the  walls  may  be  measured  in  order  to  obtain  p  cj ,  which  then 
allows  for  determination  of  the  fluid  density,  p  .  The  attenuation  coeffi¬ 
cient,  a  ,  of  the  fluid  may  also  be  measured  to  provide  in  some  instances 
information  on  fluid  viscosity.  While  such  measurements  can  also  be 
made  in  curved-wall  conduits,  the  rectilinear  conduit  promotes  greater 
accuracy. 

In  Fig.  B-12b,  transducer  assemblies  6,  7  and  8  are  conventional  ultra¬ 
sonic  transducers  which  are  positioned  to  perform  density,  sound  speed 
and  attenuation  measurements.  In  Fig.  B-l2b,  the  transducers  are  illus¬ 
trated  as  being  in  the  same  plane  as  the  velocity  measuring  transducers  2. 
These  additional  transducers  may,  however,  be  located  in  an  orthogonal 
plane  or  planes,  so  that  the  beams  interrogate  the  same  volume  element 
of  fluid  at  the  same  time  that  the  flow  velocity  of  that  element  is  being 
interrogated.  The  axis  of  such  a  composite  assembly  would  ideally  in¬ 
tersect  the  flow  channel  1  axis  and  the  flow  interrogation  beam  axis  at 
a  point  9  illustrated  in  Fig.  B-l2b. 

Turning  now  to  Fig.  B-13,  an  alternative  arrangement  for  positioning  one 
of  the  transducer  assemblies  2  is  illustrated.  In  Fig.  B-13,  as  in  most 
of  Figs.  B-14  through  B-25,  only  one-half  of  the  measurement  cell  is 
shown,  since  in  those  instances  the  other  one  of  the  pair  of  transducers 
is  identical.  In  the  configuration  of  Fig.  B-13,  a  plug  21  is  inserted 
within  the  recess  3  between  the  transducer  2  and  the  fluid  within  the  con¬ 
duit  1.  The  plug  21  is  formed  of  materials  which  will  provide  preferred 
acoustic  properties  such  as  sound  speed,  density  and  attenuation  coeffi¬ 
cient  in  the  cell  measuring  region  where  it  is  undesirable  to  form  the 
entire  cell  body  4  of  a  material  having  these  specific  properties.  Addi¬ 
tional  materials  which  may  be  employed  for  this  purpose  include  graphite, 
various  alloys,  and  plastic  such  as  polytetrafluoroethylene.  In  a  specific 
example,  the  insert  21  is  formed  of  acrylic  plastic  in  which  the  shear 
wave  sound  speed  is  about  1  mm/p.  s,  and  if  the  fluid  is  water,  then  cj 
equals  0.  5  mm/ps.  If  transducer  2  generates  a  vertically  polarized  (SV) 
shear  wave  at  an  incident  angle  where  6  0  =  30  ,  the  longitudinal  wave  in 
the  water  will  propagate  at  a  refracted  angle  9  ^  of  ~45  . 
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In  Fig.  B-14,  the  transducer  2  is  mounted  on  a  wedge  31  which  is  coupled 
to  the  outside  surface  of  the  cell  body  4.  This  provides  a  technique  for 
generating  the  oblique  wave  without  the  requirement  of  recessing  the  body 
of  the  cavity. 

In  Fig.  B-15,  a  transducer  assembly  is  illustrated  in  which  transducer  2 
is  mounted  within  a  transducer  cavity  41  which  extends  at  an  oblique  angle 
entirely  through  the  wall  of  the  cell  4.  In  this  configuration  the  angle  of 
incidence  of  a  produced  beam  does  not  depend  upon  the  ratio  cQ  and  Cj 
since  the  transducer  is  coupled  to  the  fluid  itself.  In  a  specific  example 
of  an  arrangement  for  a  transducer  in  accordance  with  Fig.  B-15,  the 
transducer  element  was  a  piezoelectric  material  of  circular  shape,  19  mm 
in  diameter,  electroded  over  a  rectangular  area  of  12.  7  mm  by  6.  35  mm. 
The  transducer  assembly  was  sealed  into  the  cavity  41  by  means  of  an 
O-ring.  In  the  specific  example  above,  the  rectilinear  dimension  of  the 
,  conduit  1  was  W  =  H  =  12.  7  mm. 

In  Figs.  B-16  and  B-17,  there  is  illustrated  a  slightly  modified  version 
of  Fig.  B-15,  in  which  the  W  by  H  flow  channel  1  is  defined  by  a  first 
square  or  rectangular  pipe  42,  and  the  transducer  channels  are  defined 
by  a  pair  of  square  or  rectangular  pipes  43  and  44  which  may  be  welded, 
brazed  or  epoxied  at  the  required  positions  and  angles,  e.  g.  ,  45  .  For 
illustrative  purposes,  Fig.  B-16  shows  both  transducer  channels  or  ports 
on  the  same  side  of  the  flow  channel  pipe,  for  a  double  traverse.  Obvi¬ 
ously,  they  may  be  on  opposite  sides,  as  is  generally  intended  in  the  other 
Figs.  B-13-B-15.  A  cell  assembled  as  in  Figs.  B-l6  and  B-17  may  be 
us  -  1  in  some  instances  as  a  probe  which  accurately  senses,  with  little 

pe  rturbation  of  the  total  initial  profile  of  a  large  pipe  or  duct,  all  the  fluid 
entering  the  cell's  inlet.  By  mounting  this  probe  such  that  its  flow  axis  is 
at  a  suitable  point  such  as  on-axis,  or  off-axis  at  radius  r,  where  v(r)  ~ 
va  for  the  Reynolds  number  range  of  interest,  a  useful  indication  of  flow 
may  be  achieved.  On  the  other  hand,  either  by  scanning  radially  with  one 
probe  or  by  averaging  the  readings  obtained  with  a  number  of  small  probes, 
a  measure  of  flow  velocity  may  be  obtained  in  a  duct  or  pipe  which  is  so 
large  that  it  cannot  be  easily  interrogated  nonintrusively  over  100%  of  the 
total  flow  channel  area,  especially  where  the  large  channel  cannot  be 
square  or  rectangular. 

In  the  type  of  application  shown  in  Figs.  B-l6  and  B-17,  while  the  trans¬ 
ducers  are  nonintrusive  with  respect  to  the  small  flow  channel,  the  ex¬ 
terior  shape  and  dimensions  of  the  probe  must  be  designed  to  minimize 
drag  and  more  generally  to  minimize  the  pressure  drop  and  other  per¬ 
turbations  introduced  by  a  small  probe  in  a  large  conduit.  In  view  of  the 
fact  that  the  obstructive  cttiss  section  of  the  inlet  portion  of  the  probe  can 
be  made  quite  small  and  streamlined,  and  because  disturbances  on  the 
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outside  of  the  probe,  downstream  of  the  inlet,  do  not  perturb  the  flow 
measurement,  and  because  the  orthogonal-walled  W  by  H  flow  channel 
is  itself  unobstructed,  this  type  of  probe  may  be  particularly  useful  in 
profiling  and  survey  measurements  in  fluids  of  much  larger  cross  sec¬ 
tion  than  the  probe. 

Another  form  of  a  probe-type  cell  is  shown  in  Figs.  B-18  and  B-19.  The 
major  portion  of  the  cell  body  45  may  be  plastic,  e.  g.  ,  acrylic,  for  use 
with  water.  To  increase  the  reflectivity  of  the  inner  surface  at  mega¬ 
hertz  ultrasonic  frequencies,  a  thin  layer  of  higher  density  material  46 
such  as  0.  1  mm  thick  stainless  steel  may  be  bonded  to  the  wall  ogposite 
the  transducers.  The  two  transducers  interrogate  the  fluid  at  45  ,  the 
angle  of  refraction  of  longitudinal  waves  mode-converted  from  SV  shear 
waves  incident  from  the  plastic  at  30  .  The  transducers  may  be  elec¬ 
trically  shielded  by  surrounding  the  outer  surface  of  the  probe  with  addi¬ 
tional  0.  1  mm  thick  stainless  steel,  48,  grounded  electrically  to  the 
electromechanical  connector  shield  47.  This  shield  47  may  be  a  solid 
(inflexible)  metal  tube,  which  contains  the  two  transducer  "high"  leads, 
and  which  can  be  used  to  orient  the  probe  by  rotation  about  its  axis,  to 
measure,  for  example,  swirl  or  cross  flow  components,  a6  well  as  axial 
flow  velocity  in  a  large  conduit.  Plastic  probe  bodies  may  also  be 
shielded  by  conductive  coatings,  such  as  conductive  epoxy,  applied  to 
their  exterior  surfaces. 

A  functionally  equivalent  construction  may  be  built  by  making  the  major 
portion  of  the  cell  in  Fig.  B-18  out  of  metal  or  ceramic,  and  by  inserting 
acrylic  or  other  low-sound  speed  buffers  as  in  the  embodiment  21  of 
Fig.  B-13. 

A  clamp-on  flowmeter  probe  is  shown  in  Figs.  B-20  and  B-20a  appropriate 
to  some  pipes  which  are  generally  of  circular  cross  section,  but  which 
are  flexible  or  malleable  enough  so  that  they  can  be  mechanically  de¬ 
formed  by  orthogonal  pairs  of  parallel  clamps  49.  In  the  region  of  mea¬ 
surement,  the  originally  circular  pipe  is  constrained  externally  so  that 
an  essentially  square  or  rectangular  flow  channel  is  provided.  Trans¬ 
ducers  may  be  coupled  directly  to  an  opposing  pair  of  flattened  pipe  walls, 
or  they  may  be  coupled  or  attached  to  one  of  the  clamps,  such  that  the 
ultrasonic  beams  are  subsequently  pressure-coupled  or  otherwise  con¬ 
ventionally  coupled  to  the  opposite  walls.  If  a  circular  pipe  of  inside 
radius  R  is  deformed  to  a  square,  the  length  of  the  interior  side  S  of  the 
resulting  square  is  readily  calculated  by  equating  interior  "perimeters.  " 
Thus,  setting  2ttR  =  4  s,  then  S  =  irR/2.  This  determines  the  required 
transducer  or  beam  width,  previously  denoted  W,  as  in  Fig.  B-14,  in 
terms  of  the  round  pipe's  interior  dimension.  As  a  numerical  example, 
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if  R  =  10  mm,  S  =  W  =  15.  9  mm.  In  other  words,  S  is  approximately  60% 
greater  than  R,  or  in  terms  of  the  original  inside  diameter  2R,  S  is  ap¬ 
proximately  20%  less  than  the  original  inside  diameter. 

The  combination  of  clamps  and  transducers  shown  in  Figs.  B-20  and 
B-20a  may  also  be  applied  to  originally  square  or  rectangular  thin-walled 
pipe,  to  help  maintain  orthogonal  flat  walls  despite  high  internal  fluid 
pressures,  which  tend  to  bow  the  sides  of  unrestrained  pipe  that  is  not 
sufficiently  stiff  by  itself. 

Comparing  the  areas  AQ  and  AgQ  of  the  circular  and  square  regions, 
respectively,  it  is  seen  that 

Aq  =  ttR2;  AgQ  =  S2  =  (rrR/2)2;  Aq/AsQ  =  4/tt  *  1.27  (B-l) 

Therefore,  the  average  flow  velocity  vQ  in  the  circular  region  will  in¬ 
crease  by  27%  to  1. 27  vQ  in  the  square  region,  if  the  fluid  is  incompres¬ 
sible.  The  mass  flow  rate,  however,  is  the  same  in  both  regions,  for 
any  fluid  confined  to  the  conduit,  whether  or  not  the  fluid  is  compressible, 
for  steady-state  flow. 

In  Fig.  B-21,  a  cavity  51  inclined  at  22.  5°  to  the  normal  to  the  flow  chan¬ 
nel  axis  reflects  the  incident  longitudinal  wave  beam  to  interrogate  the 
fluid  at  45°. 

In  Fig.  B-22,  the  internal  surface  of  the  conduit  1  is  subdivided  into  vee- 
grooves  6l  in  the  regions  between  the  transducers  2  and  the  fluid.  This 
maintains  propagation  along  a  specific  oblique  path,  and  yet  reduces  the 
amount  of  material  that  must  be  removed  near  the  intersection  of  the 
transducer  cavity  and  the  flow  channel.  Another  purpose  of  the  vee* 
grooves  is  to  reduce  the  scale  of  turbulence,  and  to  provide  supports  for 
membranes,  screens  or  acoustic  windows  used  as  flow  straighteners ,  as 
described  in  connection  with  Figs.  B-24  and  B-25. 

Figure  B-23  differs  from  Fig.  B-22  in  that  the  transducer  launches  a  wave 
normal  to  the  flow  axis.  But  this  wave,  after  impinging  on  the  vee-grooves 
6l,  is  refracted  obliquely  across  the  fluid.  When  the  ratio  cq/cj  is  large, 
as  it  is  for  metal/air  interfaces,  the  oblique  angle  is  relatively  insensitive 
to  changes  in  cQ  and  Cj,  propagation  being  in  a  direction  in  air  nearly  nor¬ 
mal  to  the  vee  surfaces,  that  is,  at  nearly  45  in  the  air.  However,  in 
Figs.  B-22  and  B-23,  care  must  be  exercised  so  that  diffraction  effects 
do  not  confuse  the  measurement. 
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Figures  B-24  and  B-25  illustrate  two  ways  of  reducing  the  undesirable 
effects  of  eddies,  which  are  generated  at  high  flow  velocities.  Since 
these  eddies  are  not  symmetrically  confined  to  the  region  interrogated 
by  the  beam,  their  effects  are  not  perfectly  canceled  by  upstream  and 
downstream  interrogations,  and  so  they  produce  nonlinear  response. 

The  degree  of  nonlinearity  depends  on  the  size  and  nature  of  the  eddy 
compared  to  H.  For  example,  small  eddies  will  not  produce  departures 
from  linearity  greater  than  1%,  when  H  is  sufficiently  large.  However, 
when  circumstances  require  that  interrogation  be  at  a  particular  angle, 
and  that  the  transducer  cavity  diameter  be  comparable  to  H,  and  if  Re 
ranges  from  '"'100  to  '"<100,000,  then  it  is  preferable  to  provide  a  means 
of  damping  or  eliminating  these  eddies.  This  means  should  also  provide 
a  high  transmission  coefficient  for  ultrasonic  energy,  e.  g.  ,  at  least  50%, 
and  often  >  90%.  In  Fig.  B-24,  a  wire  cloth,  or  a  screen  mesh,  81,  is 
used  to  dampen  the  eddies,  and  yet  provide  a  high  transmission  coefficient. 
The  wire  cloth  81  is  positioned  essentially  in  the  planes  of  the  sidewalls  of 
the  W  by  H  flow  channel  1,  i.  e.  ,  a  distance  H  apart.  The  wire  cloth,  typ¬ 
ically  made  of  stainless  steel,  may  be  silver  brazed,  soldered  or  epoxied 
in  a  small  recess  machined  in  the  sidewall,  or  it  may  be  joined  to  an  in- 
sertable  transducer  assembly.  A  suitable  mesh  size  for  use  from  ~1  MHz 
to  5  MHz  with  common  liquids,  c  ~  1000  m/s,  n  <  100  centipoise,  would 
have  a  mesh  opening  of  about  1  mm,  defined  by  wires  of  diameter  about 
0.  10  to  0.  25  mm.  Finer  meshes  are  preferred,  provided  that  they  do  not 
clog,  that  they  are  strong  enough,  and  that  they  do  not  introduce  serious 
diffraction  effects.  The  spaces  in  the  wire  cloth  encourage  '.hermal  and 
compositional  equilibrium  to  be  maintained  between  the  cavity  and  the  flow 
channel.  This  assures  that  the  sound  beam  propagates  without  refractive 
bending  on  passing  through  the  screen. 

The  arrangement  shown  in  Fig.  B-25  illustrates  an  alternative  to  the 
embodiment  of  Fig.  B-24,  in  that,  instead  of  the  screen,  a  membrane 
or  plate  91  is  used  as  an  acoustic  window.  For  operation  essentially  at 
a  single  frequency,  and  at  a  single  temperature,  the  member  91  may  be 
fabricated  or  tuned  to  resonance,  even  if  its  characteristic  impedance  is 
substantially  different  from  that  of  the  fluid  (e.  g.  ,  steel/water,  impedance 
ratio  *  30:1).  However,  for  interrogation  of  most  liquids  near  45  ,  it  can 
be  shown  that  plate-like  window  materials  like  plastics  or  graphites,  be¬ 
cause  of  inherent  attenuation,  are  not  particularly  efficient  as  resonated 
devices.  Thus,  the  thickness  of  member  91  is  generally  chosen  based  on 
mechanical  and  transmission  coefficient  requirements,  not  resonance  re¬ 
quirements.  A  typical  thickness  would  be  <,  1  mm.  For  some  cases  the 
use  of  plastics  or  graphites  is  undesirable,  but  metal  is  allowable,  some¬ 
times  as  thin  as  0.01  mm.  As  a  representative  example,  stainless  steel 
shim  stock  0.  05  mm  thick  is  appropriate  with  common  liquids  and  a  fre¬ 
quency  of  about  1  MHz  or  less. 
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When  the  fluid  is  air  or  other  low-density  gas,  a  membrane  can  be  used, 
provided  the  frequency  is  low  enough.  In  Fig.  B-25,  passageways  may 
need  to  be  provided  to  assure  that  the  fluid  (cQ)  is  the  same  on  both  sides 
of  member  91.  Otherwise  the  interrogating  beam  will  not  emerge  in  the 
fluid  parallel  to  the  beam  incident  in  the  transducer  cavity. 

The  difference  in  transit  times  between  the  upstream  traveling  beam  and 
the  downstream  traveling  beam  provides,  as  discussed  earlier,  an  indi¬ 
cation  of  the  flow  velocity  of  the  fluid  within  flow  conduit  1.  The  upstream 
and  downstream  beams  can  be  separated  either  by  initiating  the  beams  at 
different  times  so  that  the  measurements  are  essentially  made  in  rapid 
sequence  or,  preferably,  by  simultaneous  measurements. 

While  electronic  measuring  methods  can  be  employed  using  only  a  single 
transducer  on  either  end  of  the  measuring  path,  there  are  advantages  to 
be  obtained  if  two  transducers  are  used  at  each  end  of  the  path.  As  illus¬ 
trated  in  Fig.  B-26,  a  pair  of  shear  wave  transducers  may  be  orthogo¬ 
nally  polarized  with  respect  to  one  another,  that  is,  +_  45  relative  to  the 
plane  of  incidence.  One  transducer  118  is  the  transmitting  transducer, 
while  the  other  transducer  119  is  the  receiver  transducer.  The  transducers 
118  and  119  are  separated  from  one  another  by  a  layer  of  solid  material 
120  having  a  similar,  or  preferably  the  same,  characteristic  acoustic 
impedance  as  the  transducers  themselves.  For  example,  this  interme¬ 
diate  may  be  formed  of  unpolarized  or  depolarized  piezoelectric  material. 

In  some  situations  the  small  nonuniformity  of  the  intensity  of  the  plane 
wave  in  the  fluid  may  be  compensated  for  by  "shading"  the  transducer 
electrodes.  For  example,  a  small  portion  of  the  central  part  of  the 
electrode  may  be  removed,  so  that  the  corresponding  part  of  the  trans¬ 
ducer  is  not  activated.  This  would  compensate  for  a  weakened  acoustic 
field  at  the  edges  of  the  beam. 

To  illustrate  cells  built  in  accordance  with  the  above  principles,  we  show 
exploded  and  assembly  view  photographs  of  the  aluminum  cell  used  in  the 
previous  contract,  Figs.  B-27  and  B-28.  In  contrast  to  this  earlier  de¬ 
sign,  the  more  recent  design,  Figs.  2,  3  and  4,  incorporates  a  screen 
soldered  over  the  transducer  ports  (and  sight  ports,  if  present),  that  is, 
along  the  sidewalls  of  the  flow  channel.  Another  improvement  aimed  at 
reducing  eddies  is  the  welding  of  smooth  transitions  from  circular  pipe 
to  the  rectangular  flow  channel.  To  simplify  the  cell,  we  have  taken  ad¬ 
vantage  of  the  (v/c)  (p  c)  formulation  and  have  eliminated  the  separate 
transd\icer  previously  used  with  the  separate  time  intervalometer.  We 
have  also  built  an  experimental  cell  of  stainless  steel,  leaving  out  the  pc 
densitometer  probe,  since  this  probe  may  not  always  be  required.  When 
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Figure  B-28.  Assembly  View  of  Flow  Cell  Used  in  Reference  1 
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it  is  required,  it  can  be  located  just  downstream  of  the  v/c  cell,  in  a 
standard  tee  fitting. 

A  transparent  acrylic  cell  was  also  built,  so  that  bubbles,  contaminants, 
and  flow  patterns  could  be  seen  under  test  conditions.  The  cell  was  ma¬ 
chined  in  two  halves  and  then  cemented  together.  The  rectilinear  slots 
were  thus  produced  routinely,  by  milling  the  necessary  channels.  Circu¬ 
lar  transducer  ports  and  inlet  ports  were  bored  after  the  assembly  was 
cemented  together. 

To  illustrate  that  the  above  concepts  are  applicable  to  fluids  at  high  tem¬ 
perature, .  we  show  in  Fig.  B-29  a  stainless  steel  experimental  cell  tested 
up  to  300  C.  Figures  B-30  and  B-31  show  internal  details  of  this  cell, 
particularly  the  square  flow  channel,  and  the  Dowtherm  passageways  to 
maintain  the  high  temperature. 

In  cases  where  zigzag  paths  or  area-averaging  are  not  required,  the  ex¬ 
perimental  cell  of  Figs.  B-29  to  B-31  may  be  replaced  by  a  simpler  cell 
such  as  that  shown  in  Fig.  B-32.  In  Fig.  B-32,  the  flow  channel  princi¬ 
pally  is  comprised  of  standard  1-in.  stainless  steel  pipe.  The  Dowtherm 
jacket  principally  is  concentric  3-in.  pipe.  Transducer  mounts  are  also 
1-in.  pipes,  welded  at  45  ,  and  flanged  to  receive  transducer  holders.  In 
one  design  the  transducer  holders  are  of  the  top  hat  design,  fabricated  of 
titanium  because  of  its  relatively  low  acoustic  impedance  (see  Fig.  E-l). 
The  cell 8  of  Figs.  B-29  and  B-32  have  the  same  flange-to-flange  dimen¬ 
sions  for  the  process  fluid  and  the  Dowtherm  lines.  But  being  substantially 
hollow,  the  concentric  design  weighs  only  about  half  as  much  as  its  solid 
counterpart. 

A  lightweight,  relatively  inexpensive  configuration  but  retaining  the  area¬ 
averaging  feature,  consists  of  nominally  1-in.  (2.  54  cm)  stainless  steel 
square  tubes  welded  together,  Fig.  B-33.  It  contains  a  stainless  steel 
screen  (42  mesh,  0.  0055  in.  wire  diameter)  flush  with  the  flow  tube  side- 
walls.  Note  that  the  transducer  holder,  fabricated  from  a  standard  SS  304 
1-in.  screwed  cap  fitting,  contains  a  square  hole  into  which  a  mating  low- 
impedance  window  is  installed.  A  square  piezoelectric  transducer,  not 
shown,  is  coupled  to  the  window.  Inlet  and  outlet  transitions  are  also 
fabricated  of  the  standard  screwed  caps.  Overall  length  of  this  flow  cell 
is  about  1  ft  (~30  cm).  Weight  is  about  3  lb  (1.  3  kg),  considerably  less 
than  the  flow  cell  of  Fig.  2,  which  weighs  some  10  lb  (4.  3  kg). 

Readers  wishing  to  compare  the  above  flow  cell  designs  with  alternatives 
developed  by  other  investigators,  are  referred  to  an  invited  ultrasonic 
review  article  recently  published  by  one  of  the  authors,  4  and  to  the  ref¬ 
erences  cited  therein. 


Figure  B-29.  Experimental  Stainless  Steel  304  Flowmeter  Cell,  Dowtherm- Jacketed, 
Used  in  Tests  to  Optimize  Ultrasonic  Frequency  and  Pathlength,  for 
Two-Phase  Prepolymer  Flow  Up  to  ~300°C. 


Figure  B - 33.  Welded  Area-Averaging  Flow  Cell  Fabricated  of  Square  Stainless 
Steel  Tubes  and  Standard  Fittings. 
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COMPARISON  OF  ULTRASONIC  PROFILING  AND 
AREA -AVERAGING  METHODS 


Alternatives.  To  conclude  this  area-averaging  appendix,  and  to  provide 
additional  perspective,  we  compare  the  profiling  potential  of  the  trans¬ 
mission  method  (particularly  where  rectilinear  flow  is  interrogated 
obliquely  by  a  rectilinear  beam  transmitted  over  a  full  cross  section  of 
the  flow  channel,  as  used  in  the  previous  and  present  programs^®)  with 
the  profiling  potential  of  the  six  other  ultrasonic  flowmetering  alternatives 
listed  in  Table  B-l. 
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TABLE  B-l.  ULTRASONIC  FLOWMETERING  METHODS 

No. 

Method 

Measured  Parameter 

Measurand 

1 

Transmission 

Times  of  Flight 

v,  v/c  or  v/c^ 

2 

Reflection  (Doppler  Scatter) 

F  requency  Shift 

v/c 

3 

Beam  Drift 

Drift  Distance  or  Angle 

v/c 

4 

Vortex  Shedding 

Shedding  F  requency 

v 

5 

Correlation 

Correlation  Time 

V 

6 

Noise 

Acoustic  Pressure 

vn,  n  >  1 

7 

Hot  Wire 

Transit  Time  in  Wire 

Wire  Temp. 
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The  profiling  comparison  may  be  effected  by  commenting  on  or  presenting 
cell  configurations  which  sample  or  weight  the  flow  profile  so  that  the 
area-averaged  value  may  be  determined.  Method  6,  noise,  may  be  omitted 
from  further  profile  consideration  here  as  it  is  presently  limited  mainly 
to  discrete  flow  switch  applications  (on/ off,  threshold  or  leak  detection)  as 
opposed  to  wide-range  flow  measurement.  This  limit  is  due  to  the  gen¬ 
erated  noise  amplitude  and  spectrum  responding  to  several  factors  other 
than  v,  e.  g.  ,  geometry,  inlet  conditions,  scale  of  turbulence,  pipe  rough¬ 
ness,  etc. 1^.  Profiling  configurations  for  the  remaining  methods,  1  to  5, 
however,  are  compared  in  Fig.  B-34.  See  also,  Figs.  B-35andB-36, 

Transmission  Method.  Let  us  summarize  first,  the  principal  character¬ 
istics  of  the  rectilinear  cell  (Figs.  2-4,  B-34a).  Representative  curves  of 
constant  v  were  obtained  in  1926  by  Nikuradse^  and  are  reproduced  in 
Fig.  B-34a.  In  the  present  program,  reference  to  Appendix  A  shows  that 
in  the  JP-4  tests  of  December  1974,  essentially  linear  response  (O  <,0.  2%) 
was  obtained  for  M  from  ~100  to  5000  lb/hr,  or  v  from  ~0.  2  to  8  ft/ s 
(6  to  250  cm/s).  To  show  that  this  range  comprises  flow  in  the  laminar, 
transitional  and  turbulent  regimes,  we  now  calculate  the  Reynolds  numbers 


177 


1 

I 

i 

i 

i 

{ 


s 

f 

I 


(e)  Correlation  Method 

Figure  B-34.  Comparison  of  Ultrasonic  Profiling  and  Area- Averaging  Methods:  Trans¬ 
mission,  Reflection,  Beam  Drift,  Vortex  Shedding  and  Correlation. 
Concepts  Illustrated  in  (a)  and  (b)  Have  Been  Generally  Demonstrated 
in  the  References  Cited. 
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Figure  B-35.  Flow  Profile  Probe  Concepts  for  Hot  Wire  Ultrasonic 
Anemometer.  Electrical  Self-Heating  Connections  to 
Terminals  Such  As  A  and  B  Are  Not  Shown. 29 


Above:  Ultrasonic  flowmeter  duct  illustrates  "square  meter"  insert  concept  of 
area-averaging.  It  combines  a  flowmetering  region  of  square  cross  section,  an 
ultrasonic  interrogating  beam  of  rectangular  envelope,  and  inlet  and  outlet  tran¬ 
sitions  to  pipe  of  circular  cross  section.  Sheet-metal  duct  can  be  installed  to 
operate  at  hydrostatic  pressure  within  pipes  of  normal  or  heavy  schedules. 


Below:  Area-averaging  transducer  plate  provides  for  a  multiplicity  of  parallel 
pairs  of  transducers.  This  insertable  plate  may  be  welded  in  place  at  an  angle 
(e.  g.  ,  45°)  to  the  axis  of  a  large  circular  pipe.  Proximity  of  adjacent  beams, 
each  of  which  yields  a  flow  velocity  value,  vj,  is  such  that  nearly  all  of  the 
cross  section  is  interrogated  along  i  parallel  chords  in  the  generally  elliptical 
measurement  zone.  Weighting  of  each  Vj  is  proportional  to  the  W.  product. 


Two  Methods  for  Interrogating, Essentially  the  Full 
Cross  Section  of  Fluid  Flowing  in  Large  Circular 
Conduits. 
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Re  =  v  Dfr/v.  The  hydraulic  diameter  Dj^  for  a  rectangular  pipe,  accord¬ 
ing  to  Schlichting,  ^0  page  414,  equals  4  times  the  area,  divided  by  the 
wetted  perimeter.  Thus,  Dj^  =  (4)(1  )(l/2)/3  in.  =  2/3  in.  =1.7  cm  for  our 
specific  cell,  Fig.  2.  Referring  to  Ref.  1,  Fig.  2,  v  ~  0.0125  stokes  for 
JP-4  at  50  F,  and  0.  009  stokes  at  100  F.  Approximating  these  values  by 
v  =  0.  01,  the  above  area-average  flow  velocities  v  correspond  to  a  Re 
range  of  about  800  up  to  about  33,  000.  In  the  lower  range  (Re  <  2000) 
flow  is  clearly  laminar,  in  mid-range  (Re  >  2000)  transitional,  and  in 
the  upper  range,  turbulent.  The  graphs  and  calculations  of  deviations 
from  linearity  in  Appendix  A,  prove  that  with  the  cell  of  Fig.  2,  for  the 
40:1  flow  range  now  seen  to  encompass  all  three  flow  regimes,  the  recti¬ 
linear  technique  inherently  provides  for  an  area -averaging  measurement 
without  any  need  for  subsequent  profile  correction  by  a  K-factor.  Again, 
the  deviation  from  perfect  linearity  is  expressed  by  CT  <  1/4%  of  range. 

The  only  other  ultrasonic  method  presently  known  to  provide  comparably 
accurate  area-averaging,  is  a  multichord  method.^'  ”  (But  see  also  Refs. 

5  and  7.  )  So  far,  this  multichord  method  has  been  applied  to  relatively 
large  pipes,  typically  of  diameter  >  1  ft  (>  0.  3  m).  Referring  to  Fig.  B-5, 
we  note  that  to  obtain  data  over  chords  which  are  sufficiently  far  off  the 
pipe  axis,  the  usual  practice  has  been  to  drill  ports  at  prescribed  loca¬ 
tions.  In  some  instances  it  would  be  desirable  to  avoid  penetrating  the 
pipe,  e.  g.  ,  for  reasons  of  safety.  However,  for  the  usual  combinations 
of  pipe  material  and  fluid,  and  using  externally -mounted  transducers, 
the  elastic  properties,  refraction  and  mode  conversion  phenomena  pre¬ 
vent  one  from  transmitting  over  a  chord  far  off  the  pipe  axis.  A  typical 
practical  limit  is  on  the  order  of  1/4  of  the  pipe  radius.  To  avoid  this 
limit,  and  at  the  same  time  obtain  propagation  over  a  path  independent  of 
refraction,  Fig.  B-34a  suggests  the  use  of  pa.  allel-faced  internal  and 
external  wedges,  or  the  use  of  internal  reflectors.  Where  their  installa¬ 
tion  and  presence  are  allowable,  such  obstructions  would  permit  one  to 
achieve  area-averaging  by  transmission  over  preferred  chords  without 
jeopardizing  the  pipe  integrity  as  a  penetration  might.  If  obstructions 
are  intolerablei  internal  cavities  may  provide  an  acceptable  alternative.^ 

| 

Reflection  Method.  Doppler  methods  for  profiling  are  illustrated  in  Fig. 

B-34b,  and  include  range -gating, 9  traversing^  (see  also  Fig.  B-18),  and 
utilization  of  the  intersection  of  fixed  or  swept  beams. *  Potential  error 
sources,  particularly  for  the  intersecting  beam  methods,  include  non- 
uniform  beam  intensity,  nonani^orm  distribution  of  scatterers,  and  atten¬ 
uation.  In  all  cases  one  assumes  that  the  scatterers  move  at  the  local  v. 

This  assumption  is  known  to  be  net  generally  valid  for  particulate  scat¬ 
terers  in  gases,  and  may  sometimes  be  invalid  for  scatterers  in  liquids. 
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For  the  reflection  method,  as  well  as  for  Methods  3,  4  and  5  listed  in 
Table  B-l,  the  rectilinear  flow  cell  concept  sometimes  may  be  incor¬ 
porated  to  take  advantage  of  its  simpler  geometry  (at  least  regarding  the 
ultrasonic  measurements)  compared  to  the  circular  pipe. 

Beam  Drift.  This  method  reportedly  has  provided  accuracy  on  the  order 
of  5%  for  measurements  of  air  flow  in  ducts.  22  it  is  not  considered  prac¬ 
tical  for  liquid  flow  since  the  Mach  number  is  so  low.  A  schematic  is 
given  in  Fig.  B-34c. 

Vortex  Shedding.  This  method  is  usually  based  on  the  measurement  of 
the  v-dependent  frequency  at  which  vortices  are  shed  off  a  single  strut 
placed  across  the  turbulent  flow  stream.^’  ^  To  obtain  profile  data, 
particularly  for  cases  motivated  by  uncertainties  in  meter  calibration 
factor,  the  basic  concept  might  be  expanded  to  include:  a  compound  strut 
having  several  zones  of  differing  shape  or  diameter;  multiple  struts; 
transducer  arrays  to  interrogate  vortices  shed  at  different  locations  such 
as  at  different  radii  in  the  fluid,  the  response  to  specific  locations  being 
enhanced  in  some  cases  by  focusing  on  the  wakes  of  struts  at  said  specific 
locations.  Reflectors  might  be  used  in  conjunction  with  pulse  echo  inter¬ 
rogation  to  limit  each  transducer  to  respond  to  only  one  specific  wake. 

See  Fig.  B-34d,  also,  White,et  al, and  Yamasaki  and  Rubin  (1974). ^3 

Correlation.  Coulthard  (1973)  and  Ong  and  Beck  (1974)  hav^developed 
ultrasonic  correlation  flowmeters  in  the  past  few  years. Precisions 
are  reported  to  be  in  the  2  to  5%  range.  Again,  to  achieve  area-averaging, 
the  rectilinear  cell  may  be  considered  in  view  of  its  profile -compensating 
advantages  compared  to  the  usual  circular  pipe  (Fig.  B-34e). 

Hot  Wire  Profile  Probe.  In  contrast  to  Methods  1-6,  wherein  flow  ve¬ 
locity  is  determined  by  measuring  the  speed  or  amplitude  of  sound  which 
is  generated  or  transmitted  in  the  fluid  itself,  this  last  method,  7,  is 
based  on  measuring  sound  speed  in  a  foreign  sensor,  the  hot  wire 
anemometer  probe.  In  effect,  this  transforms  the  flow  profile  problem 
to  a  temperature  profile  problem. 

27 

Hot  wire  anemometers  are  well  established  as  flow  measuring  devices. 
Means  for  ultras onically  measuring  average  temperature  or  temperature 
profile  in  wires  which  are  self-heated  electrically  have  also  been  estab¬ 
lished.^®’  ^  Average  temperature  in  a  wire  can  be  determined  by  measur¬ 
ing  the  transit  time  over  the  wire,  by  reflection  or  transmission  techniques, 
or  their  combination.  Temperature  profile  in  a  wire  can  be  determined  by 
reflection  techniques.  Figure  B-35  illustrates  these  techniques. 
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Temperature  may  be  averaged  along  wires  located  in  a  plane  which  is 
perpendicular  or  inclined  to  the  flow  axis.  In  the  latter  case  the  paths 
could  coincide  with  the  Gaussian  quadrature  paths  of  Ref.  6.  Alternatively, 
temperature  could  be  profiled  in  one  or  more  segmented  wires  located  on 
a  diameter  or  on  chords. 

In  summary,  ultrasonic  flow  profiling  or  area-averaging  has  been  dem¬ 
onstrated  for  Methods  1  and  2,  with  accuracies  better  than  1%  so  far 
limited  to  the  two  transmission  methods  of  Fig.  B-34a.  The  use  of  key 
elements  used  successfully  in  the  two  transmission  methods  (rectilinear 
cell,  multiple  paths)  might  be  adapted  to  Methods  2  to  5,  but  the  actual 
potential  and  limitations  of  such  adaptations  are  not  yet  known. 


I 
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APPENDIX  C 
0  c  DENSITOMETER 


In  USAAMRDL  TR  72-66,  we  described  the  general  principle  of  ultra¬ 
sonic  reflectance,  or  impediometry  as  it  is  sometimes  called,  as  applied 
to  determining  the  characteristic  acoustic  impedance  p  c  of  aviation  fuels.  1 
We  showed  that  the  ideal  probe  material  would  have  its  own  characteristic 
impedance  Zj  close  to  that  of  the  fuel;  that  Zj  would  not  be  sensitive  to 
temperature;  and  that  the  probe  material  should  be  readily  available,  ma- 
chineable,  rugged,  impervious,  and  corrosion-resistant  or  inert  relative 
to  the  fuel,  including  contaminants.  The  geometry  of  main  interest  was  a 
cylindrical  rod  of  stepped  diameter,  to  provide  the  required  echoes  A  and 
B,  as  in  Fig.  C-l . 


Figure  C-l.  pc  Probes  With  Stepped  Diameter,  and  With  Uniform 
Diameter. 

The  present  program  utilized  a  nearly-symmetric  probe,  also  shown 
in  Fig.  C-l.  Materials  for  the  newer  probe  included  fused  silica  and 
X-cut  quartz.  The  newer  probes  generally  were  slightly  larger  than 
their  predecessors,  to  minimize  second-order  interferences  due  to 
mode  converted  pulses  generated  spuriously  at  the  sidewalls  and  at  the 
ends. 

Although  not  initially  obvious,  it  turned  out,  due  to  practical  considera¬ 
tions,  that  in  order  to  determine  the  temperature-immunity  of  the  probe, 
a  temperature-immune  echo  comparator  was  almost  mandatory.  Con¬ 
versely,  to  verify  the  echo  comparator,  a  temperature-immune  probe 
was  virtually  essential.  In  the  program,  an  iterative  procedure  was  fol¬ 
lowed  in  an  attempt  to  yield  ape  Bystem  which  was  sufficiently  immune 
to  temperature,  without  sacrificing  rangeability,  linearity  or  other  im¬ 
portant  characteristics. 

As  test  data  became  available  for  analysis  during  the  program,  it  grad¬ 
ually  became  apparent  that  the  usual  transducer  bonding  and  probe 
mounting  procedures  were  inadequate.  Tests  of  coupling  or  bonding 
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agents  such  as  propylene  glycol  (liquid),  Dow  resin  276 V9  (viscous  fluid), 
epoxy,  electrostatic  and  a  gold/tin  eutectic  solder  led  us  to  select  an 
epoxy  which  is  usable  from  -200  C  to  +100  C.  Tests  of  mounts  including 
oppositely -tapered  rods,  silicone-rubber-potted  assemblies,  epoxied-in- 
sleeve  assemblies,  and  compression  fittings  (Conax,  Swagelok)  led  us  to 
select  the  epoxied  sleeve  combination.  However,  further  test  results 
may  indicate  that  this  selection  should  be  revised.  Accordingly,  it  is 
considered  appropriate  here  to  introduce  several  alternative  p  c  probe 
designs  (Fig.  C-2),  several  of  which  have  already  been  fabricated  (Fig. 
C-3).  Some  of  these  alternatives  merely  relieve  thermal  stresses  at  the 
transducer  bond  or  simplify  the  mechanical  mounting  procedure.  Others 
magnify  the  p  c  sensitivity  by  causing  the  B  echo  to  be  influenced  by  not 
just  1 ,  but  2  or  3  solid/liquid  reflections  before  being  measured.  Also, 
the  possibility  of  substituting  a  low-Z^  rod,  such  as  certain  isopaustic 
glasses,  graphites  or  possibly  other  materials,  would  offer  another 
factor  of  2  improvement  in  sensitivity. 

The  desirability  of  improving  probe  sensitivity  and/or  stability  will  be 
understood  from  the  mathematical  analysis  given  below.  The  analysis 
also  implies  that  p  c  system  performance  may  be  improved  significantly 
by  using  a  narrowband  echo  comparator  approach,  rather  than  the 
broader  band  approach  embarked  upon  at  the  outset  of  the  present  pro¬ 
gram.  The  narrowband  approach,  as  presently  envisioned,  would 
utilize  the  regeneration  of  continuous  waves  from  the  A  and  B  echoes, 
somewhat  analogous  to  the  regeneration  of  continuous  upstream  and 
downstream  waves  which  has  now  been  successfully  demonstrated  in 
the  flow  velocimeter  portion  of  this  program. 

Next,  we  present  a  mathematical  analysis  of\the  accuracy  required  in 
the  measurement  of  A  and  B  to  obtain  p  c  within  a  specified  fraction  of 
1  %>  assuming  one  uses  the  nearly-symmetrical  probe  of  fused  silica. 

This  will  be  followed  by  a  brief  description  of  a  proposed  narrowband 
echo  comparator,  which  appears  capable  of  measuring  A  and  B  accurately 
enough  so  that  p  c  could  be  determined  in  the  range  +^  0.  1  to  0.  3%.  Com¬ 
bining  this  p  c  result  with  a  v/c  measurement  of  comparable  accuracy, 
the  mass  flow  rate  iJl  could  be  determined  with  a  probable  error  of  less 
than  +_  0.  5%. 

Now  we  shall  provide  an  error  analysis  of  the  p  c  measurement,  as  a 
function  of  the  measured  parameters.  The  purpose  of  the  analysis  is 
to  determine  the  accuracy  to  which  the  amplitudes  of  the  two  reflected 
pulses  must  be  measured. 

The  p  c  measurement  consists  of  the  independent  measurements  of  the 
amplitudes  of  essentially  identically-shaped  10-MHz  rf  pulse  bursts 
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reflected  from  a  fused  silica/ air  interface  and  from  a  fused  silica/liquid 
interface,  where  the  liquid  may  be  JP-4  fuel.  The  first-mentioned  pulse, 
A,  is  100%  reflected  and  serves  as  the  normalizing  factor  in  the  measure¬ 
ment.  The  particle  displacement  reflection  coefficient  from  the  fused 
silica  lod/liquid  interface  is  given  by 


R 


P1  Ci  -p2C2 
Pl  Cl  +  p2°2 


(C-l) 


where  Zj  =  p  j  =  acoustic  impedance  of  fused  silica,  and  Z2  =  P  2  c2  = 
acoustic  impedance  of  the  liquid,  p  =  density,  and  c  =  sound  speed. 


Equation  (C-l)  can  be  rewritten  as 


Z, 


Z 


/  A  -  B 
1  \  A  +  B 


) 


(C-2) 


where  R  =  B/A,  B  =  amplitude  of  pulse  reflected  from  the  fused  silica/ 
liquid  interface,  and  A  =  amplitude  of  pulse  reflected  from  the  fused 
silica/air  interface.  Note  that  A  and  B  vise  the  measured  quantities, 
and  for  now  we  take  Zj  =  constant.  Taking  the  logarithmic  derivative 
of  (C-2), 


2BfiA  -  2  A  5  B 

A2  _  2 

A  -  B 


In  terms  of  the  reflection  coefficient, 


_2 

A 


RCA  -  5  B 
(1  -  R2) 


If  the  A  and  B  pulses  have  mutually  random  noise, 


_2 

A 


R  |  6  A |  +  |  SB 
(1  -  R2) 


and  if  the  relative  noise  amplitudes  are  the  same  for  A  and  B, 


6  A 
A 


( C  -  3 ) 


(C-4) 


(C-5) 


(C-6) 
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In  the  case  of  fused  silica  and  JP-4,  it  can  be  shown  that  R  =  0.  86.  Thus, 


5Z2 

6  A 

Z2 

=  9.  3 

A 

(C-7) 


In  the  actual  measurement,  A  =  1000  millivolts.  Equation  (C-7)  there¬ 
fore  predicts  that  the  fractional  error  in  measuring  2^  -  p  2  c2  *s  ab°ut 
0.  93%  per  millivolt  of  noise  on  the  pulse  amplitude.  Laboratory  mea¬ 
surements  of  relative  p  c  error  vs.  fi  (A  -  B)  have  shown  about  0.  57c  per 
millivolt,  which  is  in  reasonable  agreement  with  the  above  analysis. 

From  the  above  analysis,  it  can  be  concluded  that,  to  insure  an  overall 
p  c  measurement  error  of  less  than  0.  57o,  the  "A"  and  "B"  measure¬ 
ments  should  each  be  more  precise  than  ~0.  5%.  This  can  be  accom¬ 
plished  if  narrowband  cw  techniques  are  used  and  if  standard  precision 
AC  voltage  measuring  techniques  (+_  0.  027o)  are  used.  This  approach, 
which  is  similar  in  many  respects  to  the  high-frequency  portion  of  the 
(01  -  0  2>  measurement,  is  presently  being  investigated  as  an  alternate 
and  more  precise  means  of  measuring  the  amplitudes  of  the  reflected 
waves  from  the  fused  silica/air  and  the  fused  silica/liquid  interfaces. 

Along  with  the  narrowband,  reconstituted  cw  approach  to  improving  the 
p  c  determination,  some  improvement  in  probe  design  is  anticipated. 
Referring  again  to  Fig.  C -2,  we  point  out  that  the  stepped-diameter 
probe  (a)  used  in  the  previous  contract  tended  to  crack  at  the  stress  - 
raising  corner.  The  nearly-symmetrical  improvement  (b)  has  proven 
to  be  more  rugged  and  much  better  damped  when  fully  potted,  but  it 
suffers  somewhat  from  nonidentical  crystal  bond  stresses  when  the 
sleeve  in  which  it  is  potted  expands  on  heating.  The  bond  stresses  due 
to  differential  thermal  expansion  could  be  eliminated  by  potting  only  one 
rod,  but  this  would  introduce  other  problems,  particularly  under  vibra¬ 
tory  loads.  Therefore,  we  considered  alternatives  such  as  (c),  (d), 
where  both  A  and  B  rods  cou^d  be  fully  potted,  yet  a  small  gap,  ~l/8  in. 
(~3  mm),  could  be  left  between  the  pair  of  neighboring,  matched  crystals. 
Separated  rods  also  offer  more  design  freedom  in  the  electronics,  rela¬ 
tive  to  gating  and  frequency  selection.  This  separated-rod  design  appears 
relatively  promising  for  future  work  with  the  pc  concept. 


*Commercially-available  instrument  alternatives  may  also  be  considered. 
Princeton  Applied  Research  Corporation,  for  example,  recently  intro¬ 
duced  their  Model  1 88  digital  ratiometer  (which  utilizes  coherent  detection 
principles)  for  measuring  A/B  ratios. 
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Another  group  of  probe  alternatives  is  illustrated  in  (e)  to  (j),  where 
only  one  bond  is  used,  and  the  reference  reflectors  are  readily  fab¬ 
ricated  in  probe  materials  such  as  aluminum  (subsequently  anodized) 
or  titanium.  (Eliminating  the  reference  reflector  or  probe,  not  shown, 
is  also  a  possibility.  ) 

The  probes  (k)  and  (1)  provide  double  and  triple  reflections  of  the  A  and 
B  echoes,  for  corresponding  increases  in  p  c  sensitivity. 

Trapezoids  (m)  are  included  principally  for  academic  interest.  Poten¬ 
tially,  these  could  measure  viscosity  (using  shear  wave  transducers 
and  Sg)  as  well  as  pc,  leading  to  a  determination  of  kinematic  viscosity 
and  Reynolds  number. 

At  this  writing,  the  drift  in  the  p  c  densitometer  is  principally  due  to  the 
probe.  Probe  drift,  in  turn,  is  thought  to  be  attributable  mainly  to  the 
transducer/ rod  bond(s).  Secondary  sources  of  probe  drift  are  attributed 
to  minor  changes  in  the  reflectivity  at  the  periphery  of  the  rods,  which 
influence  the  contribution  of  rays  which  strike  the  periphery  at  grazing 
incidence.  Such  secondary  sources  of  error  can  be  eliminated  by  making 
the  rod  diameter  large  enough  to  delay  both  mode  converted  and  noncon- 
verted  peripherally-reflected  waves  beyond  the  time  of  flight  of  the  direct 
waves  whose  end-wall  reflection  corresponds  to  the  fuel  pc. 

One  may  summarize  the  problems  covered  in  the  above  discussion  by 
stating  that  much  of  the  difficulty  in  determining  p  c  by  the  reflection  co¬ 
efficient  method  stems  from  the  fundamental  difficulty  of  measuring  the 
relatively  small  difference  KA-LB  between  two  relatively  large  numbers 
KA  and  LB,  when  |  R |  is  nearly  equal  to  unity. 

Instead  of  measuring  the  reflection  coefficient,  if  one  measures  the 
transmission  coefficient,  that  fundamental  difficulty  is  avoided.  However, 
one  must  avoid  two  new  errors  due  to  (1)  attenuation  in  the  fluid  by  ab¬ 
sorption  and  scattering  (turbulence)  and  (2)  diffraction  (beam  spread). 

Up  to  now,  in  this  program,  we  have  not  pursued  the  pc  transmission 
approach  precisely  because  of  the  computational  complexity  associated 
with  eliminating  the  contributions  of  absorption,  scattering  and  diffrac¬ 
tion  losses.  However,  given  the  present  state  of  the  art  of  micro¬ 
processors,  plus  the  experimental  difficulties  observed  with  the  pc 
reflection  approach,  the  transmission  approach  is  worth  considering. 

Perhaps  the  easiest  way  to  understand  the  basis  for  this  approach  is  to 
examine  the  equation  (C-8)  for  energy  transmission  Tg  between  media 


190 


of  impedances  Zj  and  Z 2  for  plane  waves  at  normal  incidence  and  the 
corresponding  graph,  Fig.  C-4. 

Tg  =  4  ZJZ2/(Z1  +  z/  (C-8) 

The  graph,  Fig.  C-4,  is  based  on  the  values  listed  in  Table  C-l.  Since 
we  are  concerned  with  transmission  (see  Fig.  C-5)  across  two  interfaces 
per  path,  we  plot  the  term  corresponding  to  Tg  ,  namely,  10  log  Tg  or 
20  log  Te  in  dB  versus  the  fluid  impedance  Z  ;>  for  four  typical  materials 
that  could  be  used  either  as  buffer  rods  or  as  the  flow  cell  structure. 
From  Table  C-l,  or  Fig.  C-4,  it  will  be  seen  that  the  sensitivity  of  20 
log  Tg  to  Z£  is  relatively  independent  of  probe  impedance  Zj,  but  in¬ 
creases  for  smaller  Z2.  Quantitatively,  the  sensitivity  (AZ2)/Z2(AdB) 
of  the  four  quasiparallel  curves  in  Fig.  C-4  is  about  1%  per  0.  1  dB  for 
Z2  ~  0.  1  g/cm^-p  s.  This  means  that  to  determine  p  c  of  the  fluid  to  1%, 
one  needs  to  measure  the  transmitted  wave  amplitude  to  1%.  To  the 
extent  that  the  effect  of  attenuation  and  diffraction  losses  can  be  elim¬ 
inated,  this  means  that  amplitude  measurement  requirements  are  about 
one  order  of  magnitude  less  stringent  for  the  transmission  approach  than 
for  the  reflection  approach. 

The  diffraction  loss  is  calculable  in  terms  of  the  Seki  parameter  S  = 
kz/a^  where  k  =  C2/f,  z  =  pathlength,  and  a  =  radius  of  buffer  rod.  ^ 

The  attenuation  loss  is  measurable  using  the  differential  path  inherent 
in  the  rectangular  cell  of  Fig.  C-5.  For  this  cell  of  dimensions  1  cm  x 
2  cm,  each  of  these  losses  may  be  expected  to  be  approximately  1  dB 
through  suitable  choice  of  frequency  (~1  to  5  MHz)  and  rod  diameter 
(~1  cm). 

The  pc  transmission  approach,  in  summary,  would  consist  of  measuring 
two  transmissions  A  and  B,  normalizing  them  to  reference  echoes  Arj 
and  Ar2  in  the  transmitter  rods,  correcting  for  diffraction^  and  atten¬ 
uation,  and  finally  computing  Tg  and  then  Z2  =  p  c  in  the  fluid. 

While  this  transmission  coefficient  procedure  appears  much  more  com¬ 
plicated  than  using  an  expression  like  (C-2),  it  will  be  recognized  that 
in  the  transmission  procedure  the  accurate  determination  of  Z2  no  longer 
depends  on  the  extremely  precise  determination  of  (A-B).  In  effect,  the 
transmission  coefficient  approach  exchanges  computational  complexity 
for  measurement  complexity.  Given  the  capabilities  of  modern  micro¬ 
processors,  this  exchange  may  significantly  favor  the  transmission 
approach  over  the  reflection  approach  within  the  framework  of  a  p  c  densi¬ 
tometer.  Furtner,  when  S<<1,  computation  is  simplified  -  see  Fig.  C-5. 
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Figure  C-4.  Transmission  Loss  for  Two  Interfaces  (Solid-Fluid-Solid)  Versus  Fluid 
Impedance  Z^,  for  Four  Buffer  Rod  Materials  of  Different  Impedances  5 
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TABLE  C-l.  LOSS  IN  dB  FOR  SOLID-FLUID-SOLID  TRANSMISSION 
PATH,  ATTRIBUTABLE  TO  IMPEDANCE  MISMATCH 
BETWEEN  Z  AND  Z 

1  Ct 

Solid 

zi 

,  2 
g/cm  -jjl  s 

Z2 

g/cm  s 

Loss  =  20  log  T 

e 

dB 

Fused  Silica 

1.  30 

0.  03 

21.  09 

0.  04 

18.  72 

0.  05 

16.  91 

0.  07 
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0.  10 

11.  53 

0.  15 

8.  6l 

0.  20 

6.  70 

Titanium 

2.  77 

0.  03 
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0.  04 
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0.  05 
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0.  07 
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0.  10 
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0.  15 

14.  20 
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0.  04 

35.  97 

0.  05 

34.  04 

0.  07 

31.  16 

0.  10 

28.  12 

0.  15 

24.  68 

0.  20 

22.  26 

f 

? 

" 


1 


i 

) 

i 

1 


( 

1 


DIFFRACTION  LOSS,  dB 


j 


f 

i 


Hypothetical  Data  Obtainable 
in  Transmission  Test,  and 

Extrapolation  to  Obtain  A0. 
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SEKI  PARAMETER16  S  =  X»z/a2  PATHLENGTH,  CM 

Figure  C-5.  Differential  Path  pc  Transmission  Cell.  Pathlengths  zA 
and  zB  in  the  fluid  may  be  in  the  proportion  1:2,  e.  g.  , 

1  cm  x  2  cm.  Received  amplitudes  A  and  B  are  attenuated 
by  absorption,  diffraction  and  mismatch  losses.  At  5  MHz 
the  absorption  coefficient  may  be  ~0.  1  to  ~2  dB/cm.  The 
diffraction  loss  for  S  <<  1,  is  essentially  proportional  to 
pathlength.  In  this  special  case,  A  and  B  may  be  easily 
corrected  (as  illustrated  graphically  in  semilog  plot)  for 
absorption  plus  diffraction  losses ,  the  extrapolated  value 
Aq  being  p  c- dependent:  Aq  =  A  /B. 
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We  conclude  this  pc  appendix  by  noting  that  in  some  instances,  such  as 
operation  only  over  a  narrow  temperature  range,  the  differential  path 
configuration  of  Fig.  C-5  may  not  be  necessary.  The  v/c  transducers 
may  suffice.  That  is  to  say,  when  the  attenuation  and  beam  spread  are 
substantially  constant,  and  if  the  transducer  piezoelectric  activity  and 
sensitivity  coefficient  are  also  substantially  constant,  then  the  amplitude 
of  the  wave  transmitted  across  the  v/c  measuring  region  may  be  inter¬ 
preted  in  terms  of  the  fuel  impedance  Z^. 


i 


i 
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APPENDIX  D 
ECHO  COMPARATOR 


The  echo  comparator,  Fig.  D-l,  consists  internally  of  a  power  Bupply, 
a  10-MHz  rf  burst  pulsed  oscillator  and  gating  circuitry,  an  echo  com¬ 
parator  receiver,  and  output  circuitry.  A  block  diagram  of  the  echo 
comparator  is  shown  in  Fig.  D-2.  Tracking  and  linearity  test  data  are 
given  in  Figs.  D-3  and  D-4. 

To  understand  the  operation  of  the  echo  comparator,  consider  Fig.  D-5, 
showing  the  nearly-symmetrical  p  c  probe,  and  broadband  or  video  echoes 
A  and  B  obtained  under  spike  excitation.  Figure  D-6  shows  that  the 
spectra  of  A  and  B  peak  broadly  below  10  MHz,  and  continue  out  at  least 
to  10  MHz.  This  means  that  if  the  p  c  crystal  is  driven  with  an  rf  burst 
containing,  say,  10  cycles  of  a  10-MHz  carrier,  i.  e.  ,  a  lp.s  wide  burst, 
the  echoes  will  closely  follow  the  drive  waveform. 

The  pulsed  oscillator,  then,  drives  the  pc  crystal  at  about  10  MHz,  at  a 
prf  of  about  200  Hz.  Probe  buffers  are  cut  to  lengths  L^  and  LB  such 
that  end  echoes  A  and  B  arrive  in  lps  time  slots  beginning  at  10|is  and 
12jjls  after  the  initiation  pulse.  These  echoes  are  selected  by  adjustable 
echo  selectors  A  and  B,  amplified,  peak  detected,  and  finally  compared 
in  sum  and  difference  amplifiers.  Both  echoes  are  amplified  equally, 
but  by  an  amount,  controlled  by  the  AGC  amplifier,  such  that  A®  2  V. 
Output  dc  levels  KA  +  LB  and  KA  -  LB  are  fed  to  the  flow  calculator  via 
coaxial  cable  connections,  where  their  quotient  is  formed,  proportional 
to  p  c  of  the  fuel.  K  and  L  are  manually  adjusted  potentiometers  set  such 
that  KA  =  LB  when  probe  ends  A  and  B  radiate  into  equal  impedances, 
namely,  dry  air. 

The  liquid  density  sensed  by  the  p  c  probe  may  be  computed  from  the  echo 
comparator  outputs,  the  probe  impedance  Z^,  and  the  sound  speed  c ^  in 
the  liquid: 

p  =  (Z^/c^)  (KA  -  LB)/(KA  +  LB)  (D-l) 

An  example  of  test  data  and  the  resulting  calculation  are  given  in  Table 
D-l  and  in  the  corresponding  graph,  Fig.  D-7. 
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pc  PROBE 


O  KA  +  LB 


Figure  D-2.  Echo  Comparator  Block  Diagram.  Legend: 

(1)  Pulser  and  Gating  Circuitry;  (2)  RF  Amplifier; 
(3)  Envelope  Detector:  (4)  Peak  Detector;  (5)  Peak 
Detector;  (6)  Summing  Amplifier;  (7)  Difference 
Amplifier. 
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RTV 


AIR 


1  p  s  /  div 
0.  1  V/div 


Figure  D-5.  Epoxy-Bonded  p  c  Densitometer  Probe,  Using  Sandwiched 
10  MHz  X-Cut  Crystal  Between  Gold-Plated  Fused  Silica 
Buffer  Rods.  Oscillogram  Showing  Reference  Echo  A  From 
Shorter  Rod,  and  Echo  B  From  Longer  Rod.  Echo  B  De¬ 
creases  From  the  Value  Shown,  When  It  Contacts  a  Liquid 
Whose  pc  More  Nearly  Matches  the  Acoustic  Impedance  of 
Fused  Silica. 
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A  SPECTRUM 


TABLE  D-l.  ECHO  COMPARATOR  TEST  DATA 


Liquid 

P 

,  3 
g/cm 

c 

/  * 
m/ s 

KA-LB 

volts 

KA+LB 

volts 

Calc,  p 
3 

c/cm 

** 

p  norm 

g/cm _ 

Acetone 

.  79 

1172 

.  318 

4.  085 

0.  86 

Benzene 

.  87 

1290 

.  380 

4.  015 

0.  95 

Water 

.  998 

1500 

.  490 

3.  905 

1.  087 

Nitrobenzene 

1.  20 

1459 

.  559 

3.  84 

1.  297 

1.  1908 

Carbon  Disulfide 

1. 26 

1149 

.  472 

3.  925 

1.  36 

1.  249 

Chloroiorm 

1.49 

986 

.  487 

3.  915 

1.  637 

1.  502 

Carbon  Tetrachloride 

1.  595 

925 

.  479 

3.  92 

1.715 

1.  575 

^Corrected  for  test  temperature,  25  to  26°C. 

3 

Calculated  value,  normalized  to  value  fo^  water,  0.  998  g/cm  .  This 


column  may  be  compared  with  first  p  column. 


DENSITY  COMPUTED  FROM 
ULTRASONIC  MEASUREMENTS,  G/CM 


Null  readjusted 


(29  April  1 974) 


NOMINAL  DENSITY,  G/CM 


0 

APPENDIX  E 

ANALYSIS  OF  THERMAL  EFFECTS 


ACOUSTIC  PROPERTIES  OF  FUELS  - 
TRANSMISSION.  REFRACTION _ 

The  sound  speed  and  density  of  several  fuels  are  given  in  Ref.  1,  Figs.  1, 

10  and  22.  The  attenuation  coefficient  may  be  estimated  from  data  in  Ref. 

11  cited  in  Ref.  1.  The  temperature-dependence  of  these  parameters  may 
be  computed  or  estimated  from  the  data  cited.  From  this  data  it  will  be 
evident  that  p  ,  c  and  p  c  of  aviation  fuels  exhibit  much  larger  temperature- 
dependences  than  elastic  solids.  Therefore,  the  reflection,  refraction  and 
transmission  across  such  solid/liquid  interfaces  will  exhibit  temperature- 
dependences  dominated  by  the  changes  in  the  acoustic  properties  of  the 
fuel. 


The  application  of  the  reflection  coefficient  to  measuring  fuel  p  c  has  been 
dealt  with  in  Appendix  C.  Let  us  deal  here  with  transmission  and  refrac¬ 
tion  phenomena.  The  energy  transmission  coefficient  T  for  plane  waves 
at  normal  incidence,  for  propagation  from  medium  1  to  2,  or  from  2  to  1 , 
is  given  by 


T  =  4  2^/12,  +  z/ 


(E-l) 


or 

T  =  (4Zi/Z2)/(Z1/Z2+  l)2 


(E-2) 


where  the  Z's  are  characteristic  acoustic  impedances:  Z^  =  Pj  ci»  ^2  = 
p  2c2’  Equivalently,  the  magnitude  of  the  transmission  loss  is  expressed 
in  decibel  notation  as 

dB  loss  =  lOlog  (1/T)  (E-3) 

T  is  determined  by  the  impedance  ratio  Z^/Z^t  as  in  (E-2).  The  pc 
nomogram  introduced  in  1965  by  one  of  the  authors  can  be  used  in  con¬ 
junction  with  a  second  graph  of  loss  vs.  Z^/Z^to  aid  in  visualizing  or 
estimating  mismatch  losses.  This  is  shown  in  the  two  graphs  comprising 
Figs.  E-l  and  E-2.  In  the  log -log  p  c  nomogram,  Fig.  E-l,  note  that 
lines  of  constant  Z  are  drawn  with  slopes  of  -1,  and  further,  that  the  ratio 
of  two  impedances  is  proportional  to  the  linear  distance  between  said  im¬ 
pedance  values,  measured  along  the  45  impedance  axis,  which  is  the 
line  of  slope  equal  to  +1.  Accordingly,  the  mismatch  from,  say,  steel 
to  water,  is  determined  by  the  ratio  of  impedance  values  between  the 
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Figure  E-l.  Impedance  Nomogram. 


intersections  of  the  45°  impedance  axis  and  the  two  dashed  iso-Z  lines 
passing  through  stainless  steel  (SS)  and  water,  namely,  Zj/Z2  ~  4.5/ 
0.15=  30.  The  transmission  loss,  Fig.  E-2,  is  9  dB.  Over  the  tempera¬ 
ture  range  and  fuel  types  of  interest  in  this  contract,  fuel  impedance  Z2 
ranges  from  approximately  0.  05  to  0.  15  g/cm^-ps.  The  corresponding 
losses  can  be  readily  estimated  for  various  Z^’s.  Values  in  the  loss  vs. 
mismatch  graph  will  be  seen  to  correspond  to  the  calculated  values  con¬ 
tained  in  Table  E-l  . 


TABLE  E-l.  ENERGY  TRANSMISSION  COEFFICIENT  T  AND 
TRANSMISSION  LOSS  CALCULATED  FOR  PLANE 
WAVES  AT  NORMAL  INCIDENCE  FOR  VARIOUS 

IMPEDANCE  RATIOS 

zi 

Z2 

T 

Loss 

dB 

1 

1 

0 

2 

.  89 

.  51 

3 

.  75 

1.  25 

6 

.49 

3.  10 

10 

.  33 

4.  81 

14 

.  25 

6.  04 

30 

.  12 

9.  04 

38 

.  10 

10.  00 

100 

.  039 

14.  07 

200 

.  020 

17.  03 

300 

.  013 

18.  78 

400 

.  010 

20.  02 

500 

.  0080 

20.  99 

1000 

.  0040 

23.  99 

2000 

.  0020 

26.  99 

4000 

.  0010 

30.  00 

Another  way  of  depicting  the  relationship  between  T,  Zj  and  Z2  was  given 
recently  by  Golis  (1974).  33 

A  useful  approximation,  valid  to  the  extent  that  Z£  <<  Zj,  expresses  T 
as  follows: 


T  *  4Z2/Zj 


(E-4) 
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It  can  be  shown  that  for  normal  incidence  the  sound  pressure  transmission 
coefficient  expressing  the  ratio  of  pressure  received  at  a  perfectly  termi¬ 
nated  (matched  backing)  transducer  of  acoustic  impedance  Z j,  compared 
to  that  transmitted  from  an  identical  transducer,  is  also  approximated  by 
4  Z2/ZJ,  if  attenuation  and  beam  spread  are  either  negligible  or  compen¬ 
sated  for.  In  this  simplified  case,  the  fuel  characteristic  impedance  is 
directly  proportional  to  the  received  amplitude,  or  received  voltage  A: 

P2c2  =  KAZl  (E-5) 

where  K  is  a  constant  that  includes  the  transducer  activity  and  sensitivity. 
This  equation  can  be  shown  to  apply  also  to  transmission  between  appro¬ 
priate  buffer  rods  of  impedance  Zj.  Thus,  we  have  a  transmission  coef¬ 
ficient  impediometer  or  p  c  densitometer. 

In  passing,  it  is  interesting  to  note  that  for  fluids  where  Z2  <<  Zj,  if  one 
utilized  a  parallel  branch  stilling  chamber,  as  in  Fig.  11,  p.  41,  in  which 
p  c  is  determined  by  Eq.  (E-5),  and  provided  that  attenuation  in  the  fluid 
is  either  negligible  or  compensated  for,  a  more  sensitive  transmission 
coefficient  measurement  of  fuel  p  c  would  be  obtained  than  if  the  reflection 
coefficient  method  of  Appendix  C  were  attempted.  Further,  for  gaseous 
fuels,  by  measuring  c  both  in  the  freestream  as  well  as  in  the  parallel 
branch,  the  average  freestream  gas  temperature  T£  would  be  obtainable 
from  the  square  of  the  sound  speed  ratio  and  the  measured  or  thermo¬ 
statically  controlled  branch  temperature  T^: 

Tf  =  Tb  <  w2  (E-61 

The  freestream  gas  impedance  (p  c)£  is  then  given  by 

(p c)(  =  (POh^Tr;  =  KAZj  cb/c£  IE-7) 

Furthermore,  the  masf  flow  rate  in  the  freestream  can  be  determined  as 
the  product  of  freestream  gas  impedance  and  Mach  number: 

Mf  =  (pc)f  (v/c)f  =  KAZ1cbv/cf2  (E-8) 

Finally,  at  least  for  a  limited  range  of  gas  compositions,  either  the  gas 
average  molecular  weight  M,  or  the  M/7  ratio,  could  be  determined  from 
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2 


(E-9) 


the  measured  and  temperature  T^: 

“  =  TR„  VCb2  °r  "  RoTb/Cb 

where  7  =  average  ratio  of  specific  heats  and  RQ  is  the  universal  gas  con¬ 
stant.  From  these  M  or  M/7  values,  the  energy  content  of  the  gaseous 
fuel  could  be  computed  in  Btu/lb.  Consequently,  the  rate  of  energy  trans¬ 
fer  in  Btu/s  could  be  metered  for  diagnostic  or  control  purposes. 

Refraction  phenomena  are  of  interest  in  this  contract  for  two  reasons. 
First,  a  flush-mounted  wedge  comprises  one  potential  approach  to  avoid¬ 
ing  the  generation  of  eddies  in  the  v  measuring  region  of  the  flow  cell. 
Second,  dating  back  to  the  inception  of  the  previous  program,  the  possi¬ 
bility  of  measuring  flow  with  clamp-on  or  externally -mounted  transducers 
has  been  recognized  as  desirable,  if  it  could  be  accomplished  without 
sacrificing  performance. 

At  oblique  incidence,  the  solid/liquid  transmission  and  mode  conversion 
effects  are  more  complex  than  at  normal  incidence.  Having  dealt  with 
these  energy  partition  and  mode  conversion  problems  in  a  recent  NASA 
report,  let  us  limit  the  present  discussion  simply  to  phenomena  that 
are  readily  quantified  using  Snell's  Law: 


c^ /sin  9  j 


c2/sin02 


(E-l  0) 


where  9  j  =  angle  of  incidence  (in  the  solid),  and  9  2  =  refracted  angle  (in 
the  fuel).  The  change  in  6  2  caused  by  temperature -dependent  or 
composition-dependent  changes  in  C2  are  readily  understood  from  a  graph 
such  as  Fig.  E-3,  plotted  for  9  j  =  45°  and  for  one  value  of  Cj,  namely, 

Cj  =  1500  m/s.  (This  is  approximately  the  value  of  the  shear  wave  ve¬ 
locity  in  ATJ  graphite.  ) 

A  generally  more  useful  graph  of  Snell's  Law  is  given  in  Fig.  E-4  for 
three  representative  values  of0j,  namely,  30  ,  45  ,  60  .  (The  inter¬ 
ested  reader  will  find  a  linearized  version  of  this  graph  in  Golis  (1974), 
where  the  sinusoidal  "distortion"  of  9  j  and  9  2  axes  enables  the  velocity 
ratio  parameters  to  be  plotted  as  straight  lines.  ) 

Of  the  several  consequences  of  changing  refracted  angles,  two  of  the  most 
obvious  and  most  important  are:  (1)  the  axially -projected  path  changes, 
and  (2)  the  beam  could  miss  the  receiver  entirely.  Referring  to  Fig.  E-5, 
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500  I  00  000  I  200  1000  1  soo 

<  2 >  m/s 

Figure  E-3.  Refraction  Graph  for  9  j  =  45°,  =  1500  m/s 


Figure  E-5.  Wedges  Beveled  for  9  j  =  30,  45,  60°.  Top,  Wedges  Positioned  for  Com¬ 
munication  Over  Path  at  45°  to  Flow  Axis.  Bottom,  Refraction  Changes 
for  cj  =  1500  m/s,  C2  =  980  and  1020  m/s  at  30°C  and  20°C,  Re¬ 
spectively.  Transducer  Diameter  =  2a. 
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the  axial  pathlength  L  is  given  in  terms  of  the  pipe  or  duct  inside  dimen¬ 
sion  D  as 


L  =  D  tan  9  ^  (E-ll) 

Note  that  in  Fig.  E-5,  the  centers  of  the  wetted  surfaces  of  the  wedges 
are  along  a  45°  angle  to  the  pipe  axis.  However,  the  transducers  them¬ 
selves  of  diameter  2a  may  be  at  different  angles.  Under  the  constraint 
that  the  receiver  transducer  wedge  intercepts  at  least  half  the  beam,  and 
under  the  assumptions  that  tne  beam  does  not  spread  at  all,  the  limits  on 
9  £  were  computed  in  terms  of  D/a  for  three  orientations  of  transducers 
(9  i  =  30°,  45°,  60°).  By  trigonometrically  analyzing  Fig.  E-5,  the 
limits  on  0  2  are  readily  shown  to  be 


limiting  9  2  = 

1  +  2  &/\IT  D 

for  9 

=  30°  'j 

limiting  6  2  = 

1  +  a  \TZ/ D 

for  9 

.  r  °  I 

=  45 

>  (E-l  2) 

limiting  9  ^  = 

1  +  2  a/D 

for  9 

) 

The  calculated  values  in  Table  E-2  show  that  for  small  D/a,  more  latitude 
in  9  2  can  be  tolerated.  The  influence  of  9  ±  is  more  pronounced  for  small 
D/a.  See  also  Fig.  E-6.  By  combining  this  information  with  the  Snell's 
Law  graph,  Fig.  E-4,  the  required  D/ a  can  be  determined  in  terms  of  the 
range  of  C2/cj  to  be  accommodated.  However,  data  reduction  must  take 
into  account  the  effective  L  -  see  Eq.  (E-ll).  Except  for  the  unique 
wedge/fuel  case  where  the  fractional  change  in  c^  per  C  is  identically 
equal  to  the  fractional  change  in  C2  per  °C,  L  would  change  with  tempera¬ 
ture,  and  in  any  event,  would  change  with  composition,  requiring  compen¬ 
sation.  For  this  reason,  the  screen  is  preferred  over  the  wedge. 

Referring  again  to  Fig.  E-5,  the  bottom  part  of  this  illustration  shows 
explicitl>  the  change  in  path  due  to  a  change  in  C2.  For  purposes  of  illus¬ 
tration  let  us  analyze  the  case  where  9  j  =  60°;  Cj  =  constant  =  1500  m/s; 

D  =  constant  =  1  in.  =  2.  54  cm;  and  c2  =  1000  m/s  at  25°C,  with  a  tem¬ 
perature  coefficient  of  -4  m/s  per  C.  The  refracted  angles  are,  there¬ 
fore,  6  2  (20  C)  =  36.  08°,  9  7  (30°C)  =  34.46°.  The  axial  paths  are  com¬ 
puted  by  (E-ll)  to  be:  L  (20bC)  =  0.  729  in.  ;  L  (30°C)  =  0.  686  in.  The 
temperature-dependence  of  L  over  this  range  is  AL/AT  =  -0.43 %  per  °C, 
where  AL  =  0.  043  in.  =1.1  mm  corresponds  to  the  10°C  temperature 
difference. 
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TABLE  E-2.  MINIMUM  AKD  MAXIMUM  BOUNDS  ON  REFRACTED 
ANGLE  9  2,  FOR  NONSPREADING  BEAM,  SUBJECT 
TO  CONSTRAINT  THAT  AT  LEAST  HALF  THE  BEAM 
IS  INTERCEPTED,  AS  A  FUNCTION  OF  D/a  AND 
6  j  =  30,  45  AND  60°. 


9  x  =  30 


2  min 


2  max 


0  =45 


2  min 


2  max 


•l- 

60° 

9  _  . 

2  mm 

9  , 

2  max 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 


22.  9 
31.  6 
35.  4 

37.  6 

38.  9 

39.  9 

40.  6 
41.1 

41.  5 


65.  1 
57.  6 
54.  2 
52.  2 
50.  9 
50.  0 
49.4 
48.  9 
48.  5 
48.1 


- 

67.  5 

16.  3 

59.  6 

27.  9 

55.  8 

32.  9 

53.  5 

35.  7 

52.  1 

37.  4 

51.  0 

38.  6 

50.  2 

39.  5 

49.  6 

40.  1 

49.  2 

40.  7 

48.  8 

- 

71.  6 

0 

63.  4 

18.4 

59.  0 

26.  6 

56.  3 

31.  0 

54.  5 

33.  7 

53.  1 

35.  5 

52.  1 

36.  9 

51.  3 

37.  9 

50.  7 

38.  7 

50.  2 
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Figure  E-6.  Angular  Bounds  on  Refracted  6  to  Assure  Reception  of  At  Least  Half  the 
Beam,  As  a  Function  of  D/a  and  Q  1  .  Notation:  See  Fig.  E-5,  Top. 


A  second  consequence  of  refraction-induced  path  changes  is  that  the  path- 
length  in  the  wedges  changes,  too.  Let  us  trace  the  central  ray  from  the 
upper  (transmitter)  wedge  to  the  lower  (receiver)  wedge.  Denoting  the 
wedge  paths  W£q  =  w  (20°C)  and  w^q  =  w  (30°C),  it  can  be  shown  that  the 
difference  in  wedge  paths  is  W30-w20  =  ^L/2  =  0.  866  AL,  or  0.  866  6  L 

for  a  1°C  change,  where  5l  =  0.  1  AL. 


If  not  compensated,  the  above  refraction-induced  changes  in  fluid  path- 
length  and  wedge  pathlength  would  lead  to  errors  in  v,  c,  c2,  v/c  and 
v/c2  approximately  as  follows  for  an  assumed  1°C  temperature  increase 
from,  say,  25°C  to  26°C. 


Error  in  v:  The  phase  difference  A <f>  decreases  by  0.43%. 


V?T 


Error  in  c:  The  fluid  path  decreases  from  V  +  (L  +  6  L)  to 
Numerically,  the  fluid  path  decrement,  L  •  6l,  is  approximately  0.  0025 
in.  The  wedge  path  increment  is  0.  0037  in.  If  we  calculate  the  corres¬ 
ponding  change  in  transit  time  At  =  t25-t26  =  (wec*8e  Pa*h  increment)/ cj  - 
(fluid  path  decrement)/c2«  the  result  is  At  <  1  ns,  which  is  a  negligible 
fraction  of  the  total  transit  time.  Therefore,  in  the  present  example  the 
errors  in  v,  v/c  and  v/c^  are  dominated  by  the  change  in  fluid  axial  path 
L.  Percentage  errors  for  the  three  terms  are  substantially  identical, 
about  0.  4%  per  °C. 


To  avoid  thermal  refraction  errors  associated  with  flush-mounted  wetted 
wedges  (which  errors  would  obviously  be  unacceptable  for  temperature 
excursions  much  larger  than  1°C)  one  would  generally  need  to  continuously 
monitor  c^  and  C£  and  then  compute  v  or  related  terms,  taking  into  ac¬ 
count  the  consequences  of  the  changing  c^  and  c^.  Alternatively,  one  can 
choose  to  avoid  refraction  changes  by  keeping  Q  £  constant  by  means  of  an 
oblique  cavity.  This  was  our  choice  in  the  present  program. 


ACOUSTIC  SHORT  CIRCUIT 


As  seen  from  Fig.  E-l,  c  in  metals  is  typically  about  4  to  5  times  faster 
than  in  the  fuels  of  interest.  Obviously,  there  arises  the  possibility  of  a 
metal-borne  wave  arriving  at  the  receiver  before  the  fuel-borne  signal. 
Even  if  the  metal  path  were  somewhat  longer  than  the  fuel  path,  the  travel 
time  in  the  metal  can  still  be  shorter  than  in  the  fuel,  hence,  the  term, 
acoustic  short  circuit.  Less  obvious  are  the  consequences  of  the  many 
vibration  modes  that  can  be  supported  by  a  structure  as  complex  as  a 
typical  flow  cell.  For  example,  suppose  that  a  short-duration  transmitted 
pulse  is  inadvertently  coupled  to  the  metal  cell  or  pipe  wall.  A  long  pulse 
train  will  be  generated  because  the  different  modes  propagate  at  different 
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velocities  in  the  metal  structure.  Thus,  the  short  circuit  rings  on, 
interfering  with  the  fuel-borne  wave.  In  the  case  of  an  external-mount, 
reference  to  Fig.  E-l  will  show  that  the  transducer -metal  pipe-transducer 
path  unfortunately  tends  to  be  well  matched  compared  to  the  transducer- 
fuel-transducer  path. 

Still  less  obvious  is  the  fact  that  even  when  the  signal  to  short  circuit 
noise  ratio  is  high  enough  so  that  the  signal  is  readily  identified  and 
distinguished  from  the  short  circuit  noise  on  an  oscilloscope  (e.  g.  ,  S/N  = 
10), the  error  in  measuring  flow  velocity  will  generally  be  unacceptably 
large  and  temperature-dependent.  We  will  show  below,  by  means  of  a 
simple  vector  analysis,  that  a  S/N  ratio  approaching  100:1  is  required 
if  short  circuit  errors  are  to  be  suppressed  below  1%. 

Referring  now  to  the  vector  (or  phasor)  diagram,  Fig.  E-7,  we  identify  a 
reference  vector  pointed  "east",  pairs  of  flow  velocity  signal  vectors  Su 
and  S^  representing  upstream  and  downstream  waves,  and  noise  vectors 
N.  Relative  to  the  phase -measuring  system  described  in  the  body  of  this 
report,  full-scale  flow  corresponds  to  a  phase  difference  of  180  between 
Su  and  S^. 

Let  us  assume  the  S/N  ratio  =  10.  We  will  consider  the  potential  phase 
errors  caused  by  collinear  vectors  N  adding  vectorially  to  Su  and  at 
10%  of  full-scale  flow,  and  also  near  full-scale  flow. 

At  10%  of  full-scale  flow,  Su  and  S<j  are  symmetrically  rotated  +  9°  from 
the  reference  vector.  The  small  circular  orbits  centered  on  the  S  tips 
represent  the  allowable  contributions  of  N.  It  appears  reasonable,  from 
a  reciprocity  argument,  that  N  will  be  the  same  for  both  upstream  and 
downstream  paths.  For  simplicity,  at  the  10%  of  full-scale  flow,  we  show 
the  situation  where  the  N  vectors  both  point  "north".  The  resultant  wave 
appearing  at  the  upstream  transducer  is  denoted  Ru,  and  at  the  down¬ 
stream  transducer,  R<j.  Intuitively,  N  seems  to  cancel.  But,  unfortun¬ 
ately,  calculation  of  the  angles  involved  shows  that  while  the  true  angle 
or  phase  difference  between  Su  and  is  18.  00  ,  the  angle  between  Ru 
and  Rj  is  17.  82°,  the  error  being  1%.  At  lower  flow  velocities  the  error 
diminishes,  but  does  not  vanish  except  at  zero  flow. 

At  the  higher  flow,  say,  at  about  94%  of  full  scale,  and  with  N's  taken 
"westward",  a  S/N  of  10  leads  to  an  error  of  7%  A  S/i1  ratio  of  100 
diminishes  this  error  source  to  0.  6%. 

At  a  given  moment,  the  angular  relationship  between  the  S's  and  the  N's 
will  depend  on  v  and  the  sound  speeds  in  the  liquid  and  solid  media.  But 
these  speeds  depend  on  temperature  and  composition.  Clearly,  a  rather 
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Figure  E-7.  Vector  Diagram  Showing  Errors  Due  to  Acoustic  Short 
Circuit  Noise  N  Adding  Vectorially  to  Flow  Velocity 
Signals  Su  and  Sj.  For  S/N  =  10,  Vectors  are  Drawn 
for  Flow  Velocities  of  10%  and  About  94%  of  Full-Scale 
Flow.  The  Magnitude  of  the  Resultant  Percentage  Error 
is  ( |  A  0  g  -  A  0  ^  |  /A  0  g )  (100%).  Three  Triangular 
Waveforms  at  Bottom  Show  Notation  for  Analysis  of 
Pulse  Signals  in  the  Presence  of  Interfering  Acoustic 
Short  Circuit  Noise.  See  also:  NASA  CR-112313  (R), 
Ref.  15,  pp.  C-l  to  C-6. 
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high  S/N  ratio  is  required,  preferably  100:1  or  higher,  to  assure  accu¬ 
racy  of  1%  or  better  over  a  wide  range  of  conditions. 

At  the  present  time,  ultrasonic  flowmeters  are  also  available  from  other 
manufacturers,  wherein  pulse  techniques  are  used,  rather  than  continuous 
waves.  For  the  reader  who  may  question  whether  the  foregoing  vector  or 
phasor  argument  applies  to  the  case  where  short  circuit  interferes  with 
measuring  arrival  times  of  pulses  transmitted  across  the  fluid,  we  offer 
the  following  brief  analysis. 

Let  us  assume  that  "arrival  time"  is  taken  to  be  the  midpoint  of  the  lead¬ 
ing  edge.  This  means  that  attenuation  effects  are  being  eliminated  by 
using  methods  such  as  automatic  gain  control,  with  fixed  threshold,  or 
by  using  a  variable  threshold  which  is  continually  reset  to  50%  of  the  peak 
received  pulse  amplitude,  or  by  other  means. 

The  received  pulse  waveform  may  have  various  shapes  such  as  exponen¬ 
tially  enveloped  sine  wave  burst,  video  triplett,  etc.  For  simplicity  we 
treat  the  case  of  a  triangular  pulse.  This  makes  it  easier  to  explain  the 
interfering  effect  of  noise  N  upon  the  desired  signal  S. 

Referring  to  the  lower  dr:’v.  ings  in  Fig.  E-7,  we  see  that  corresponding 
to  a  triangular  pulse  of  period  T,  amplitude  S  and  noise  of  amplitude  N, 
the  threshold,  set  at  S/2,  will  be  reached  too  early,  or  too  late,  by  the 
interval  6  ,  where  6  is  the  timing  error.  Algebraically,  6  =  +_  (T/4)/ 

(S/N)  =  +_  tr/(0.  8  S/N),  where  tr  =  rise  time  of  pulse.  The  maximum 
error  in  determining  time  differences  upstream  minus  downstream  is 
25.  In  general,  timing  errors  will  vary  between  zero  and  26,  depending 
principally  on  the  period,  geometry,  S/N,  the  materials  used  for  the 
conduit  and  fluid,  their  temperature  coefficients  of  sound  speed,  and  the 
temperature  range.  Numerically,  for  T  =  1  ps,  6  =  25  ns  if  S/N  =  10, 
and  6  =  2.  5  ns  if  S/N  =  100.  For  T  =  0.  lps,  6's  would  be  proportionately 
reduced  to  2.  5  ns  and  0.  25  ns  for  S/N's  of  10  and  100,  respectively. 

It  is  sometimes  argued  that  the  short  circuit  N  is  a  constant  and  therefore 
can  be  accounted  for  by  calibration  of  the  device.  This  view  would  appear 
to  be  supported  by  the  definitions  of  repeatability  and  reproducibility.  Ac¬ 
cording  to  Ref.  31,  these  terms  express  the  closeness  of  agreement  among 
repeated  measurements,  approaching  from  the  same  or  opposite  directions, 
respectively,  but  under  the  same  operating  conditions. 

But  in  actual  practice,  minor  changes  in  temperature,  composition,  or 
dampening  on  the  conduit  interior  or  exterior  (e.  g.  ,  bubbles,  scale, 
deposits)  will  influence  N  and/or  S.  Therefore,  departures  from  cali¬ 
brated  performance  will  be  observed,  unless  both  the  time  (or  phase) 
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and  amplitude  relationship  of  N  and  S  remains  constant,  or  unless  S/N  is 
so  large  that  variations  in  S/N  are  negligibly  small. 

CELL  DIMENSIONS  AND  INTEGRITY 

Cell  dimensions  x  increase  as  temperature  increases  by  the  amount  AT 
according  to  the  well-known  relation 

x  =  x  (1  +  a- AT)  (E-13) 

where  a  is  the  coefficient  of  thermal  expansion.  See  Fig.  E-8,  especially 
curve  2.  For  AE5I  type  304  stainless  steel,  a  =  10.  4  x  10"  per  °F  =  18.  7 
x  10“^  per  °C.  For  operation  from  -65°F  to  +l60°F  (-54°C  to  +7l°C)  the 
ffAT  term  amounts  to  (10.  4  x  10”  )  ( 1 6 0  +  65)  =  2.  34  x  10”  .  From  mid¬ 
temperature,  47.  5°F,  dimensions  will  therefore  change  by  no  more  than 
+_  1.  17  mils  per  inch.  The  flow  channel  area,  nominally  1.  000  in.  x  0.  500 
in.  =  0.  500  sq  in.  ,  therefore  changes  maximally  by  +_  [1.  0017  x  .  5005- 
.  500]  or  +_  1.  35  x  10“^  sq  in.  This  is  a  fractional  change  of  +_  0.  27%. 

The  distance  between  velocimeter  transducers  is  nominally  3  in.  This 
distance  changes  by  +3  x  1.  17  =  +  3.  51  mils,  or  about  0.  1%.  The  effect¬ 
ive  axial  pathlength  L,  nominally  1  in.  ,  changes  by  j^l.  17  mils,  or  again 
about  0.  1%. 

From  Eqs.  (11)  and  (12)  on  page  26,  one  sees  that  the  v/c  formulation 
suppresses  the  effect  of  changing  cell  dimensions.  But,  even  in  the  case 
where  v  is  measured  (to  be  multiplied  by  p  ),  the  dimensional  charges  are 
all  well  below  1%.  If  required,  a  simple  temperature  measurement  could 
be  used  to  generate  the  small  correction  term  to  reduce  such  thermal 
effects  well  below  0.  1%. 

Regarding  cell  integrity,  different  models  of  our  flow  cells  made  of  AISI 
304  have  been  subjected  to  temperature  from  cryogenic  levels,  -200  C, 
up  to  +300  C,  at  pressures  up  to  several  hundred  psi,  for  periods  ranging 
from  hours  to  months.  No  structural  failures  have  occurred.  The  parti¬ 
cular  flow  cell  most  used  in  this  program.  Fig.  1,  includes  several  welded 
joints,  and  soldered  joints  between  the  flow  channel  walls  and  the  screens. 
The  solder  alloy  composition  was  97  Sn,  3  Ag,  fluxed  with  a  mixture  of 
zinc  and  ammonium  chlorides.  The  joint  is  reported  to  comply  with  Mil 
Standard  454B.  During  soldering  the  screens  were  held  in  position  using 
opposed  alumina  wedges.  Similar  solder  joints  have  been  exposed  to 
liquid  nitrogen  and  boiling  water  thermal  shocks  with  no  ill  effects. 


I 

Figure  E-8.  Expansion  Characteristics  of  Metals,  Alloys 

and  Ceramics.  From  "Materials  and  Techniques 
for  Electron  Tubes",  by  Walter  H.  Kohl,  ©  1959 
by  Litton  Educational  Pub.  ,  Inc.  Reprinted  by 
Permission  of  Van  Nostrand  Reinhold  Company. 


The  last  temperature-dependent  effect  to  be  analyzed  in  detail  in  this 
appendix  concerns  the  theoretical  transmission  coefficient  from  the 
X-cut  crystal  to  the  Si02  buffer  rod,  across  the  thin  epoxy  bond  layer. 
The  basis  for  this  analysis  is  contained  in  texts  such  as  Kinsler  and 
Fry  (1962)  and  a  similar  analysis  recently  has  been  reported  by  Canella 
(1974)  for  liquid  coupling  layers.^,  35 


The  epoxy  used  in  this  analysis,  Tra-Con's  Type  2115,  was  cast  into  a 
block  of  approximate  dimensions  1  Li.  x  1  in.  x  3/8  in.  Sound  speed  was 
measured  at  room  temperature,  liquid  nitrogen  temperature,  and  at  boil¬ 
ing  water  temperature.  For  comparison  purposes,  calculations  were 
also  carried  out  for  a  common  plastic,  acrylic.  The  acrylic  sound  speed 
data  of  Arnold  and  Guenther  (1966)  covering  30  C  to  1  00  C  were  supple¬ 
mented  by  our  measurements  at  25  C  and  -200°C.  Specimen  thicknesses 
were  measured  at  the  temperature  extremes,  but  the  ~1%  change  in  thick¬ 
ness  was  neglected  in  computing  transmission  coefficients,  as  the  expan¬ 
sion  effect  is  dwarfed  by  the  changes  in  sound  speed. 

2 

The  impedance  of  X-cut  quartz  is  1.  5  g/cm  -ps;  that  of  fused  silica  is 
1.  3  g/cm2-(i  s. 

Consider  a  sandwich  model.  The  outer  members  have  impedances  Zj  and 
Zy  The  middle  layer  impedance  is  Z2,  its  thickness  is  I,  and  its  wave 
number  is  k.  The  energy  transmission  coefficient  at  normal  incidence  is 
of  the  form 


T  = 


4Z1Z3 


(Z1  +  Z3)2  COS2  kl  +  (Z2  +  Z1  Z3/Z2)2  Sin2  k| 


(E-14) 


The  analysis  of  the  temperature-dependence  of  T  at  10  MHz  proceeds  in 
two  steps.  First,  we  calculate  ki  as  a  function  of  temperature:  ki  = 

2  it  10^i/c.  Second,  we  calculate  T,  and  the  corresponding  loss  in  dB, 
as  a  function  of  ki.  Results  are  listed  in  Tables  E-3  and  E-4,  and  are 
plotted  in  Figs.  E-9  and  E-10.  Also,  from  Fig.  E-l,  zacryiic  «  0.  3 
g/cm^-ps.  We  will  use  this  value  for  ZepOXy>  too. 
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We  proceed  to  analyze  thermal  effects  in  the  bond  layer  as  follows.  To 
estimate  i  ,  we  note  that  the  crystal  and  the  rod  ends  are  fine  ground  to  an 
rms  (root  mean  square)surface  finish  between  about  8  and  l6p  in.  They 
are  epoxied  together  under  pressure,  so  the  average  bond  thickness  prob¬ 
ably  does  not  exceed  twice  the  roughness.  Therefore,  i  is  almost  cer¬ 
tainly  less  than  40p  in.  ,  or  1  pm.  (We  will  neglect  the  thermal  expansion 
of  i .  )  Next,  from  Fig.  E-9,  the  ordinates  of  the  bottom  dashed  curve, 
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TABLE  E- 3.  PARAMETER  k I  AT  10  MHz  AS  A  FUNCTION 
OF  TEMPERATURE  AND  THICKNESS  OF 
ACRYLIC  AND  EPOXY 


kf  at 

kl  at 

■gH 

Temp. 

c 

1=1  pm 

3  pm 

Material 

°C 

m/s 

(0.  04  mil) 

(0.  1  mil) 

Acrylic 

-200 

2900 

.  0217 

.  0650 

.  217 

II 

+  30 

2700 

.  0233 

.  0698 

.  233 

It 

+  100 

2470 

.  0254 

.  0763 

.  254 

Epoxy 

-200 

2890 

.  0217 

.  0652 

.  217 

II 

+  25 

2422 

.  0259 

.  0778 

.  259 

II 

+100 

1614 

.  0389 

.  1168 

.  389 

Figure  E-9.  Graph  of  Parameter  ki  at  10  MHz  as  a  Function  of 
Temperature  and  Thickness  of  Acrylic  and  Epoxy. 
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Figure  E-10.  Graph  of  Transmission  Loss  Versus  ki 
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TABLE  E-4.  ENERGY  TRANSMISSION  COEFFICIENT  AND 
TRANSMISSION  LOSS  FOR  SANDWICH  MODEL 
OF  TRANSDUCER  BONDED  TO  ROD 


z  z  z 

1  3  2 

(g  rams  /cm^ -microsecond) 

k  1 

(radians) 

T 

(energy  ratio) 

Loss 

(dB) 

1 

.  5  1 

.  5  0 

1 

0.  01 

0.  99446 

0.  02413 

0.  03 

0.  95223 

0.  21258 

> 

f 

0.  10 

0.  64282 

1.  91911 

0.  3 

0.  01 

0.  99942 

0.  00252 

0.  03 

0.  99484 

0.  02247 

' 

t 

0.  10 

0.  94571 

0.  24242 

1 

0 

0.  01 

0.  99998 

0. 00009 

0.  03 

0.  99984 

0.  00069 

> 

'  \ 

f 

0.  10 

0. 99827 

0.  00752 

1.3  0.1 

0.  01 

0.  99015 

0.  04299 

0.  03 

0.  95371 

0.  20584 

\ 

t 

0.  10 

0.  67295 

1.  72017 

0.  3 

0.  01 

0. 99441 

0.  02435 

0.  03 

0.  99053 

0.  04132 

> 

f 

0.  10 

0. 94859 

0.  22922 

1 

0 

0.  01 

0.  99489 

0.  02225 

0.  03 

0.  99480 

0.  02264 

\ 

' 

i  \ 

1 

0.  10 

0.  99381 

0.  02697 

from  -54  to  +71°C,  increase  from  ki  =  0.  028  to  0.  030.  Now,  using  Fig. 
E-10,  the  corresponding  increase  in  transmission  loss  is  found  to  be 
0.  004  dB  for  the  solid  middle  line  (X-cut  quartz  rod)  and  0.  003  dB  for  the 
dashed  middle  curve  (fused  silica  rod).  These  computed  temperature- 
induced  increases  are  negligibly  small,  corresponding  to  echo  amplitude 
changes  of  0.  04  and  0.  03%,  respectively,  for  temperature  excursions 
from  -54  to  +71°C. 

Observed  thermal  effects  on  the  p  c  probe  in  air  are  on  the  order  of  0.  1% 
per  °C.  We,  therefore,  conclude  that,  to  the  extent  that  the  above  model 
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is  valid,  the  observed  temperature -dependence  of  echo  amplitude  changes 
in  the  present  p  c  probes,  is  not  attributable  to  changes  in  kl  or  sound 
speed  in  the  bond  layers.  Mechanisms  not  included  in  the  above  model, 
such  as  stress  or  strain,  possibly  contribute  to  the  observed  thermal 
effects.  Independent  of  the  cause,  if  the  effects  are  reproducible,  com¬ 
pensation  based  on  a  temperature  measurement  of  the  p  c  probe  would  be 
straightforward. 

VELOCIMETER  TRANSDUCERS 


The  velocimeter  transducers  are  made  of  a  modified  lead  metaniobate 
formulation,  with  a  Curie  point  above  400  C.  Diameter  is  0.  75  in. 
Resonant  frequency  is  5  MHz.  We  have  operated  such  transducers,  when 
epoxied  to  an  ATJ  graphite  backing  member,  before,  during  and  after 
cyclic  immersion  in  liquid  nitrogen  and  boiling  water.  We  have  also 
evaluated  the  transducer  material  alone  in  activity  and  sensitivity  tests 
up  to  320°C.  One  of  the  high-temperature  test  fixtures  is  shown  in  the 
photograph  of  Fig.  E-ll.  Operation  generally  has  been  reliable  despite 
these  thermal  dwell,  cycling  and  shock  exposures.  Therefore,  the 
velocimeter  transducer  assemblies  are  considered  more  than  adequate 
for  the  present  application. 

p  c  DENSITOMETER  PROBE 

The  p  c  densitometer  probe  consists  of  a  mounted  assembly  of  two  cylin¬ 
drical  rods  of  fused  silica,  one  wetted,  one  dry.  Each  rod  is  epoxied  at 
its  chrome -gold-metallized  end  to  an  X-cut  10-MHz  crystalline  quartz 
transducer.  The  quartz  transducer  operates  piezoelectrically  in  its  alpha 
phase  up  to  573°C.  (In  some  experiments  the  rods  were  made  of  X-cut 
quartz,  oriented  transversely  the  same  as  the  X-cut  transducer  crystal 
to  minimize  differential  expansion  between  fused  and  crystalline  quartz.  ) 

The  temperature  coefficient  relating  to  the  X-cut  quartz's  frequency 
stability  is  -20  x  10"^  per  °C.  However,  the  coefficients  relating  to 
activity  and  sensitivity,  djj  and  g33>  have  not  been  measured,  to  our 
knowledge,  for  this  material.  One  purpose  of  the  reference  echo  A  in 
the  p  c  probe  is  to  compensate  for  changes  in  d^  and  that  are  pre¬ 
sumed  to  occur  due  to  thermal  effects  or  aging. 

The  effect  of  temperature  on  the  difference  in  attenuation  in  the  dry  and 
wetted  rods  may  be  computed  based  on  published  data  of  the  attenuation 
coefficient  or  vs.  temperature  (Krause  et  al,  1970).  See  Fig.  E-12.  For 
the  temperature  range  of  interest,  -54  C  to  +71  (L  a  decrease^  some¬ 
what,  as  temperature  increases,  by  about  3  x  10”  dB/cm  per  C,  at  20 
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MHz.  The  round  trip  difference  in  rod  lengths  in  our  p  c  probe  is  typi¬ 
cally  about  1.  2  cm.  Thus,  for  excursions  of  +_  62.  5°C  from  the  mid¬ 
temperature,  8.  5°C,  the  attenuation  difference  at  20  MHz  would  not  ex¬ 
ceed  +  0.  02  dB.  At  10  MHz,  the  frequency  used  in  the  p  c  probe,  atten¬ 
uation  changes  would  not  exceed  0.  01  dB,  or  +_  0.  1%,  a  negligibly  small 
value.  However,  in  the  unlikely  event  where  significant  temperature 
differences  exist  between  the  two  rods,  the  attenuation  difference  would 
increase  somewhat.  It  is  concluded  that  temperature-dependent  attenua¬ 
tion  effects  in  the  fused  silica  rods  do  not  contribute  significantly  to  p  c 
errors. 
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Figure  E-12.  20-MHz  Longitudinal  Acoustic  Attenuation  Versus 
Temperature  for  Some  Ohalcogenide  and  Oxide 
Glasses,  After  Krause  vt  al.*  ' 
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LIST  OF  SYMBOLS 


B 

BCD 


c 

D 

°1 

f 

H 

K 


K,  L 


k 

L 

1 


M,M  ,MZ 


M 


M 

M 

N 

R 


Echo  amplitude 

Zeroth  order  Fourier  coefficient 

DC  term  in  Fourier  series 

Transducer  radius 

Echo  amplitude 

Binary  coded  decimal 

Sound  speed 

Diameter 

Hydraulic  diameter 

Frequency 

Height  dimension 

Flow  profile  meter  coefficient 

Echo  amplitude  adjustment  factors 

Wavenumber,  2  tr  /X. 

Axial  pathlength 
Thickness 

Measurements  of  phase  shift 
Mass  flow  rate 

•  # 
dM/ dt  =  time  rate  of  change  of  M 

Average  molecular  weight 

Noise  amplitude;  Term  in  Fourier  series 

Radius;  Reflection  coefficient 
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LIST  OF  SYMBOLS  (cont'd) 


Re 

R 

o 

S 

T 

t 


T,  T 

e 


V 


v 

v 


W 


X 

Z 

a 


7 

0 


A 

0 

\ 


Reynolds  number 
Universal  gas  constant 
Signal;  Perimeter 
Temperature;  Period 
Time 

Energy  transmission  coefficient 

Time  intervals  downstream,  upstream 

Voltage  signal  output,  Simmonds  Precision  densitometer 

Flow  velocity 

Area-averaged  flow  velocity 
Width  dimension 
Cell  dimension 

Acoustic  characteristic  impedance 

Attenuation  coefficient;  Expansion  coefficient 

N-th  order  Fourier  phase  angle 

Phase  shifts  in  fluid 

Phase  shifts  in  Receivers  1  and  2 

Specific  heat  ratio 

Timing  error 

Difference 

Angle  of  transmission 
Wavelength 

231 


jw.aurw.-**  < 


LIST  OF  SYMBOLS  (cont'd) 


Kinematic  viscosity 

Density 

Time  interval 

Phase  shifts  of  carrier  frequency 

Phase  shifts  of  the  fundamental  modulation  frequency 


Angular  frequency 


